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« SBUBCIBf 

BxtMMione  of  Frofesaor  Crocoo*o  fimdanent«X  thaoi^otloiil 
of  tibt  ^ttieaoaoa  of  eoatbustioa  iaotubllitjr  in  liquid  pf’qpollant  rool'iit 
motors  bars  boon  mds«  Coaplsto  snslysss  of  the  specific  Q$ams  of 
eentrstod  ooeibustieu  «t  tstIous  ststlons  and  distributed  combustion 
aloOK  the  chamber  <^0  are  includod  in  the  present  report#  Additlcsitl 
analjsos  including  m the  characteristics  of  tangential  modes  of 
pressure  osoiXlatioas  are  nov  under  uagr«  and  uLll  be  repcnrtad  upon  idien 
eoaplete# 

During  the  three*aontn  period  cohered  bgr  this  r sporty  wide  hat  ccn> 
tinned  ma  the  construction  and  inatallatlon  of  a nonopropellant  rocket 
notcr  and  constanWata  feed  <qrstem#  Nearly  aU>  sarrlcing  and  auxiliary 
aqoipmant  has  been  completed  and  installed#  with  aereral  exeaptionsy 
all  instrumentation  components  hare  been  receiredy  and  calibration  and 
installmtlon  of  the  inEtruments  is  in  progress#  Dslays  in  delirsry  of 
eritieal  system  ccaponents  hai«  necessitated  minor  resehedalingy  and  it 
is  doubtful  that  operation  of  the  sonstanWate  syatam  er  the  mompee^ 
pellant  rocket  motor  sill  be  aceespliahed  before  the  end  of  the  calandar 
year. 

The  Frlneeton«IfIT  pressure  transducer  has  undergone  enrlrcnmental 
testing  on  a rocket  motor  at  the  NAd  Lesis  Laborstoryy  and  appeeucs 
satisfactory  for  monopropellaat  rocket  tests  from  this  standp€»lnt# 
Failure  of  a pielnq>  ocoinred  at  a oban^r  tanperature  of  ap|!^^imtely 
6200^F  under  "screaming*'  eonditions  sufficiently  s«svsre  to  esnise 
ture  of  the  sater-oooled  test  rocket  and  some  design  i?ciMXicutlcns  have 
been  initiated  in  order  to  inorease  the  operational,  temperatiu!^  range 
of  the  transducer# 

Static  calibration  tests  of  the  differential  lrinc£ton<^I^f  wit 


rniimimmmoii  . 


SOMMARY  (Coat»d.) 

ifttra  bagua  »ft«r  soma  dol«y  dae  t®  l«alcas«  flroa  the  ts*an»due€ir«e  back- 
I»r«siure  chamb<nr»  uid  the  remits  oS  the  first  series  of  teats  iadicate 
satlsfaoterjr  llnearlt^r  and  reproducibility.  Opcnrational  testing  in 

I 

MCA  rocket  chaiabere  and  static  ca3.ibratio»  at  Priucetm  idU.  continue 
daring  the  next  report  period. 
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II.  HfTHODOgriON 
A.  Ob^leot 

BuA«r  Ceatract  HOas  52*713“8  has  been  undertaken  as  & pajE^t  of  the 
Jet  propulsion  reaeareh  prcgrase  of  the  D^rtment  of  Aeronautical  ^«> 
iaoerlng  at  Princeton  to  *>conduet  an  investigation  of  the  generral  prob- 
len  of  oonbustioQ  instability  in  liquid  propellant  rocket  engines. 

This  program  shall  consist  of  theoretloal  analyses  and  experimental 
verification  of  theory.  The  ultimate  objective  shall  be  the  oollection 
of  sufficient  data  that  shall  permit  the  rocket  engine  designer  to  pro- 
duce povsr  plants  idiloh  are  relatively  free  of  the  phenomena  of  insta- 
bility. Interest  shall  center  in  that  form  of  \mstable  operation  which 
is  characterised  \ssr  high  frequency  vibrations  and  is  commonly  known  as 
* screaming  * <*9 

B«  History 

Interest  at  Princeton  in  the  problem  of  combustion  instability  in 
liquid  propellant  rocket  motors  waa!  given  itgpetus  by  a Bureau  of  Aero- 
ikutics  <tyupo3lu»  held  at  the  Naval  Research  Laboratory  on  the  7th  and 
8th  of  December  1950.  This  Interest  resulted  in  theoretical  analyses 
by  Professors  N.  Sunmerfield  and  L»  Crocco  of  this  Center. 

Professor  Sumnerfield’s  wrark,  "Theory  of  Bnotable  Combustion  in 
Liquid  Propellant  Rocket  Systems"  (JARS,  Sept*  1951),  considers  the 
effects  of  both  inertia  in  the  liquid  propellant  feed  lines  and  combus- 
tion chamber  capacJ.tim.ee  with  a constant  ©cabustioa  ttae  lag,  and 
applies  to  the  case  of  low  (up  to  about  200  cycles  per  second)  fre- 
quency oscillations  s<»ietims3  called  "chugging". 

Professor  Crocco  advanced  tho  conc^t  of  the  pressure  depondenca 
of  tho  time  lag  in  uid-1951|  his  papar,  "Aspessts  of  CembuafcJ.oa  Stabil'^ 

ity  in  Liquid  Propellant.  Rocket  Motors"  (JABS,  Nov*  19$I  asd  Jae-Fob 
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B«  History  (Cont'd.) 

19^2),  presents  the  f undos  ontals  resulting  from  this  concept 4,  end  «®u\3yaie8 
the  cases  of  low  f*^equency  instability  with  noaopropeliants,  low  frequency  | 

instability  with  btpropellantc  and  high  frequency  instability,  with  co«-  | 

bustion  ecmcentrated  at  the  end  of  the  combusti«as  chamber* 

Desiring  to  submit  the  concept  of  a pressure  dependent  time  lag  to 
e3q>erimental  test  a preliminary  propoisAl  was  made  by  this  Dhiirersity  to 

i' 

the  Bureau  of  Aeronautics  in  the  summer  of  1951  and,  followiag  a formal 
request,  a revised  proposal  was  sabmitted  idiieh  resulted  in  the  present 
contract  dated  30  April  1952* 

Analytical  studies  of  dlstiributed  combustion  had  been  carried  on  in 
the  meantime  under  Professor  Croeco  and  within  the  sponsorship  of  the 
Guggenheim  Jet  Propulsion  Center  by  S.I*  Cheng  and  were  issued  as  his 
Fh.D*  Thesis,  "Intrinsic  High  Fi‘equeney  Combustion  Instability  in  a ! 

liquid  Propellant  Rocket  Motor",  dated  April  1952* 

Time  was  devoted,  in  anticipation  of  the  contract,  during  the  first  ! 

third  of  1952  to  coistructing  facilities,  securing  personnel,  and  plan-  [ 

ning  the  e:3q)erlmental  approach* 

During  the  first  thz'se  month  period  of  the  contract  year,  personnel  ' 

and  facilities  at  the  new  James  Forrestal  Research  Center  were  assigned, 
and  the  initial  phases  of  the  experimental  program  were  planned  in  sors 
detail* 

A constant  rate  mcnopropellant  feed  system  was  completely  designed  I 

1 

and  preliminary  designs  of  the  ethylene  oxide  rocket  motor  and  the  in- 
strumentation systems  were  worked  outc  features  of  the  projec- 

ted systems  included  a pulsi33g  unit  to  cause  oscillations  in  propellant 

flow  rate,  a water-cooled  strain-gage  pressure  pickup  designed  for 
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IHTRCDUCriOH  (Cont«<i.) 

B.  ELRtory  (Cont'd.) 

flU0h  asouz!t.ing  in  th«  rocket  chaiBber^  and  sev-araX  possible  methods  for 
dynamlo  measorement  of  an  oscillating  propellant  flow»rato« 

Searches  were  made  of  the  literature  for  sources  of  information 
on  combustion  instability  and  ethylene  oxides  and  visits  to  a number 
of  activities  working  on  liquid  propellant  rocket  combustion  instabil- 
ity problems  were  made  for  puxposes  of  familiarisation  with  equipment 
and  results* 

Subsequent  efforts  are  presented  In  detail  in  this  report# 

G*  personnel 

The  personnel  assignments  for  this  project  have  been  changed 

since  the  first  quarterly  report  was  issued*  Personnel  at  the  present 

time  ara  assigned  as  followst 

1*  Professor-in-charge«  Dr*  Luigi  Crooeo  (One- third  time) 

2*  Asaiatant  Fsc'ofescor  (Theoretical  Studies)^  Dr*  S#I*  Cheng 
(Half-time) 

3#  Research  Snginews  Dr.  Jwrry  Orey 

li*  A8s*t*  Research  Sogineert  (Position  vacant) 

armduate  Assistants  (2)t  Q.B.  Hatthewa 

D.T*  Harrja  (Half-time) 

6.  A Technician  and  two  meohanioa 

A critical  gap  exlata  >in  the  position  of  Assistant  Research  Engin- 
eer# l(^h  effort  has  been  devoted  to  the  procurement  of  a qualified 
perscaa  with  some  background  in  electronic  instrumentation,  but  without 
success  to  date* 

Supporting  services  are  furnished  as  required  bf  the  project# 
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n.  IMTRQMCTION  (Coitfc'd.) 

C*  Schedule 

^le  Bajor  problem  encountered  to  date  In  the  investigation  hae  been 
the  di.ffioultj  in  obtaining  materials  and  equipment  idthout  prohibitively 
ICBig  delivery  delays*  Nearly  every  major  piece  of  equipment  received  has 
been  delayed  from  one  to  six  months*  and  many  ..ritioal  items*  long  over- 
due at  this  mrlting*  are  not  expected  for  several  veeks  or  months  due  to 
Bumufacturing  or  scheduling  difficulties  on  the  part  of  the  several  sub- 
oontractcrs*  In  scane  instances*  it  has  bean  necessary  to  use  makeshift 
equipment  in  order  to  avoid  excessive  delays  in  the  program*  but  it  is 
hoped  that  most  of  the  critical  items  will,  be  available  by  the  end  of 
the  year* 

Poor  deliveries  have  necessitated  a change  in  schedule*  Feed  systcsi 
operatl<mal  tests  will  probably  begin  in  January  after  completion  of  In- 
dividual component  tests*  and  the  first  monopropellant  rocket  motor  runs 
are  e]q>ected  early  in  February*  Part  of  the  individual  instrument  tests 
and  calibBfations  can  be  run  concurrently  'with  the  rocket  motor  and  feed- 
system  by  utilisation  of  the  facilities  offered  by  the  coig;oDeut  teat 
panels* 

Dssign  of  the  bipropellant  feed  .'iiystem  is  being  held  iq>  until  pre- 
liminary results  of  monopropellant  system  tests  are  available*  in  order 
that  the  use  of  any  unsatisfactory  components  < can  be  avoided  in  the  new 
design*  It  is  felt  that  the  small  amount  of  time  idiich  may  be  lost  as  a 
result  of  this  procedure  will  bo  regained  vilth  interest  by  using  the 
monopropellant  system,  in  effect*  as  a test-bed  for  the  bipropellant  ln-» 
stability  studiaa,  and  basing  design  of  the  latter  on  63q)erifflace  gained 
by  actual  operation* 
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in.  THECREl’ICAL  SIUDIES 

k bri©f  V(fT±m  of  the  fundsiaental  theory  vf  combustioa  l?sstabllJ.ty 
in  liquid  propellant  rocket  aotops  as  set  forth  hsr  Br.  Luigi  Oroeeo  is 
prssentsd  here  for  general  information  purposes«  Also  included  as 
pendiees  are  complete  presentations  of  the  tuo  anaiyses  tdiose  results 
were  given  very  briefly  in  the  first  quarterly  report.  These  appendi- 
ces are  referred  to  occasionally  in  the  ansulng  qualitative  discu8sicffii«» 

Ibe  theory  of  homogeneous  combustion  as  it  exists  today  is  an  ex- 
tremely con^lioated  onei  ho«ever>  this  is  a moderately  simple  problem 
ahen  conqpared  to  the  heterogeneity  of  the  combustion  prooese  vhich  takes 
place  in  the  chamber  of  a rocket  motor.  In  normal  rocket  combustion^ 
processes  involving  atomisation,  liquid  mixing,  diffusion,  heating, 
vaporisation,  and  chemical  reactions,  all  highly  sensitive  to  injector 
and  chamber  geometry,  produce  such  an  Intricate  picture  that  no  attseq}t 
can  be  made  to  predict  the  behavior  of  the  propellants  betueen  the  in- 
jector and  the  rocket  nozsle.  It  can  be  said,  however,  that  these  pro- 
cesses are  dependent  upon  physical  parameters  such  as  injection  rate 
and  pattern,  temperature,  pressure,  velocity,  mixture  ratio,  etc.,  and 
thus  the  processes  will  take  place  at  a constant  rate  only  if  these 
physical  factors  are  constant  in  time.  Now,  it  seems  that  with  the  con- 
ditions prevailing  in  ordinary  rockets,  such  a possibility  can  never 
occur,  so  that  the  local  physical  parameters  always  fluctuate  more  or 
less  with  time,  and  the  process  rates  idth  them.  As  a result,  no  per- 
fectly smooth  combustion  can  exist,  and  fluctuations  of  larger  ss^ler 
amplitude  are  always  present.  If  fluctuations  are  large,  we  arbitrarily 
call  the  ocmbustion  trough**,  if  net,  it  is  "smooth".  Now,  normal  com- 
bustion may  be  defined  as  a combustion  in  wlilch  procooses  going  on  ot  a 
Certain  location  have  no  influence  on  the  presosoeo  at  a sufficiently 
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THEOREnCAL  STUDIES  (Ooai»d,) 

distent  locaticm}  i»e.,  fluctuations  in  local  physical  parajnet©**8  occur 
at  randoBS  with  respect  to  location  in  the  changer*  Since  the  local 
process  rates  are  thus  varying  a random  manner*  the  a'trerage  effect 
of  the  variaticos  over  a sufficiently  large  time  or  space  will  be  sero* 
This  holds  true  regardless  of  the  amplitudes  of  the  local  fluctuations* 
and  hence  it  follows  that  stresses  on  the  walls*  heat  transfer  eoeffi* 
dents*  and  nechanical  vibrations  of  the  system  will  be  practically  un- 
affected by  roughness  of  the  combustion« 

NOW  suppose  that  some  factor*  external  to  the  processes  themselves* 
causes  the  fluctijatims  of  one  or  more  of  the  physical  parametars  to  be 
coordinated  so  that  the  local  process  rates  are  no  longer  distributed  at 
randon*  and  a correlation  now  exists  between  fluctuations  at  tro  arbl» 
trarily  distant  points  in  the  chamber.  If  appropriate  conditions  exist* 
these  coordinated  fluctuations  can  feed  back  to  the  physical  parameters 
themselves*  creating  a situation  of  self- amplification  or  "instability"* 
Even  though  the  amplitude  of  these  fluctfiations  may  be  no  larger  than 
that  of  the  randcmi  "rough  burning"  variations*  the  effects  of  the  organ- 
ised oscillation  on  stresses*  heat  transfer*  and  mechanical  vibrations 
are  far  more  severe*  since  biui  averaged  effects  of  such  fluctuations* 
contrary  to  the  case  of  normal  combustion*  are  not  aero*  These  organ- 
ized oscillations  can  be  sufficiently  severe  to  caitse  failure  or  burnout 
of  the  rocket  and  instrumentation  in  time  intervals  of  the  order  of  a 
fraction  of  a second*  and  this  is  v . reason  for  the  difficulty  in  ob- 
taining a sufficient  mass  of  reliable  experimental  informatlone  In,  view 
of  this  lack  of  an  experimental  basis*  formulation  of  a theory  which 
■ describes  these  organized  fluctuatimc*  or  Instabilities*  must  be  iiuide 

by  deductive  rather  than  inductive  reasoaing.  This  is  done  by  seni’ch- 
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III,  THgaREriCAL  STTOISS  {Cottfe«d,) 

lag  for  wAchanlms  \<hioh  Edght  h&  r«0potis5.bl4r  t<^  the  coordination  of 
ptgrsioal  parasaetera,  fomuliating  tkrocgh  theoretical  davelopmenta  the 
eonditioRS  for  appearance  of  inetahilltys  predicting  the  efi^cta  of 
Tariatl^a  in  the  geonetry  of  the  ^aten  and  other  qtiantltieff^  and 
iinally  collaring  with  the  available  ej^arimsntal  inforeatlon* 

Several  tjnpes  of  instability  e«n  be  defined,  depending  m the 
canee  for  coordlrievticn  of  fluctuations  and  on  the  affected  parasteters* 
The  first  type  Investigated  ms  first  observed  by  von  Karaan  and  his 
gronp,  and  ia  caused  by  pressure  oscillations  idiieh  are  set  \m  by  the 
offset  of  balance  between  burning  and  ^chauat  rates,  resulting  in  var> 
latiens  of  the  injection  rates*  This  type  of  instability  usually 
occurs  at  frequencies  lower  than  100  cycles  per  second  and  ie  called 
"low  frequency  instability"  or  "chugging"* 

Another  type  of  instability  idiioh  can  occur  Independently  of 
variations  in  injection  rate  depends  on  mve  notion  In  the  gases  and 
Its  effect  an  the  combustion  prooess*  This  phenonenon  is  ooaq>arable 
to  acoustic  resonance,  and  the  characteristic  frequencies  are  close 
to  the  natural  frequencies  of  the  chanbar.  Since  these  natural  fre* 
qpencies  are  generally  of  the  order  of  hundreds  or  thousands  of 
cycles  per  second,  this  typo  of  instability  is  usually  designated  as 
"high  ihrequency  instability",  or  "screaming"* 

WKT9  propagation  through  the  propellants  in  the  feed  system  end 
the  coordination  of  injection  rate  fluctuatl«is  produce  a third  type 
of  instability  known  as  hydraulic  resonance,  idtlch  usually  occurs  at 
Arequemoy  ranges  someidiere  between  chugging  and  scrcamlngu 

Other  possible  coordinating  cames  can  be  found,  aud  t&  each  of 
them  oorrespcsids  a different  type  of  potsai.ial  combustion  instability* 
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in,  THEORETICAL  STUDIBS  (Cont»d,) 

The  first  two  tjipes  here  been  subjected  to  theoretlosl  trestmenty  the 
basic  assuMptlons  and  fundanental  results  of  ndxloh  will  now  be  discussed 
without  going  into  the  rather  coi^lioated  details  of  the  analysis. 

As  has  been  done  in  many  other  classes  of  phenomena,  a theoretical 
criterion  for  instability  can  be  obtained  by  considering  ideal  systeias 
with  the  processes  taking  place  as  they  would  during  normal  combustion, 
and  then  applying  a disturbance.  If  the  disturbance  effects  fall  to 
sero  with  time,  we  conclude  that  the  syst«a  is  stable;  if  not,  the  system 
is  either  neutral  or  unstable. 

Thus  a prerequisite  for  the  application  of  this  criterion  is  a 
oerbaln  knowledge  of  the  normal  ceaibustlon  process,  about  which,  as  pro- 
▼lously  discussed,  so  little  is  known.  Fortunately,  however,  it  turns 
out  that  In  most  cases  the  processes  need  not  be  known  in  detail,  and 
scHoe  kind  of  synthetic  picture  is  sufficient.  If  we  consider,  for  in- 
stance, the  history  of  given  small  portions  of  two  propellants  which 
are  destined  to  react  together,  we  realise  that  the  processes  idilch  take 
place  from  the  instant  of  injection  to  the  moment  idien  the  evolved  hot 
gases  begin  contributing  to  the  chamber  pressure  cannot  have  any  direct 
effect  on  instability,  From  this  moment  on,  however,  the  process  of 
development  of  the  products  of  reaction  may  create  the  organised  os- 
cillations ^Ich  are  generally  recognized  as  instability.  In  a rocket 
mobor,  the  reaction  proceeds  more  or  less  gradually  until  it  is  approx- 
imately complete.  Now,  if  for  the  given  portions  of  prc^ellants  this 
mathematically  complicated  continuous  process  is  replaced  by  a hypothet- 
ical discontinuous  process,  we  obtain  the  considerable  simplification 
that  only  two  quantities,  the  time  delay  and  the  space  delay,  are  si;if- 
ficient  to  schematically  represent  the  combustion  process  of  the  given 

portions  of  propellantsj  Of  course,  the  same  simpllficatijMi  can  be 
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Made  fop  evesTy  ofchar  portion,  and,  finally  a satisfactosry  pictnr«  of  nor~ 

^ mal  eombnatico  can  be  obtainfad  1:^^  attributing  to  every  portion  of  each 

propellant  a tine  lag  and  a spaoe  delay,  the  latter  usually  being  defined 
vectorially.  In  other  words,  we  can  make  the  approximation  that  the  gas 
generated  at  each  point  in  the  combustion  chanber  was  Introduced  at  a 
well-defined  point  in  the  injeotioa  system  at  a well-defined  tine  before, 
for  each  propellant o 

The  closeness  with  \diich  the  approximation  of  a discontinuous  com- 
bustion can  fit  the  actual  continuous  process  depends  on  the  particular 
ease  concerned*  It  is  the  opinion  of  this  contractor  that  for  mai^  pro- 
pellant combinations,  the  rat»>dete(mining  processes  are  of  a physical 
^ nature  (e*g*,  atomisation,  heating  and  nixing)*  For  these  cases,  the 

actual  production  of  burned  gases  takes  place  in  a tine  interval  suf- 
ficiently short  so  that  the  assumption  of  discootlnuous  reaction  is  not 

* 

too  far  from  the  truth*  However,  if  the  controlling  rate  is  the  rate  of 
a chemical  process  which  is  slow  even  at  high  temperatures  (e.g*,  devel- 
opment of  free  oxygen  from  dissoolatiaa  of  NO)  then  the  approximation 
will  become  worse*  Unfortunately,  since  precise  information  on  the  act- 
ual combustion  mechanisms  and  their  rates  is  not  available,  the  pcsto> 
lation  of  a well-defined  time  delay  seems  to  be,  even  in  these  eases, 
the  only  reasonable  working  assunq>tlon* 

The  approximation  that  the  propellants  give  no  contribution  to  the 
pressure  before  they  are  suddenly  transformed  into  the  final  products 
has  the  interesting  consequence  that  from  a fluid-dynamical  point  of 
view  the  flow  can  bo  considered  to  consist  only  of  burned  gases,  generated 
by  appropriately  distributed  soTa*ces,  and  therefore  having  teuperature 
C'jtnditions  close  to  those  of  the  ccnsta^it  prossuro  combustiono  Phyaieally, 
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this  result  can  perhaps  be  accepted  as  sufficiently  s.pprariiBate  only 
if  the  Wned  gases  are  In  a state  of  Intense  circulation  (or  good 
'''recirculation")  back  to  the  InJeotoTf  tdiere  otherwise  the  heat  ex- 
changed with  the  incoming  propellants  would  create  a ten^erature  de- 
fect. (On  the  other  hand;  if  the  recirculation  is  poor  the  afore- 
said assumption  about  the  nearly  unifom  gas  temperature  needs  a 
correction.)  In  this  respect  the  schematic  picture  of  coabtustlon  in 
a rocket  differs  ooag>letely  from  that  of  an  air  ccsabustor.  Substan- 
tlally«  the  flow  in  a rocket  can  be  considered  as  the  flow  of  a gas 
with  constant  stagnation  teap^'ature  and  variable  flow  rate»  while  the 
flow  in  the  air  combustor  is  a flow  with  constant  rate  and  variable 
stagnation  teBg>erature« 

Considering  now  the  effect  of  disturbances  on  this  assumed  picture 
of  c<mbastion«  we  see  that  during  the  time  lag  between  injection  of 
our  given  portions  of  propellants  and  the  evolution  of  hot  gases»  the 
effects  of  the  physical  factors  on  the  process  rates  art*  changed  by  the 
disturbances  applled«  and  hence  the  tine  and  space  delays  for  the  given 
portions  of  propellants  will  be  altered.  A complete  description  of 
this  situation  is  again  hai^tered  by  the  total  lack  of  any  quantitative 
knowledge  of  the  processes  concerned^  so  that  it  becomes  necessary  to 
postulate  some  plausible  relationship  between  the  delays  and  the  his- 
tory of  the  physical  factors  involved.  For  exai^le^  if  the  effects  of 
only  pressure  and  temperature  of  the  gases  surrounding  the  propellant 
are  considered;  m can  make  the  schematised  assumption  that  during  a 
certain  initial  time  interval  (before  mixing  of  the  propellants), 
pressure  and  tespcrature  have  no  effect  on  the  process,  and  that  dur- 
ing the  regaining  portion  of  our  overall  time  lag,  the  rato  of  the 
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pro<H>8«  is  oniforwly  affected  pressure  and  tewperature  mriations, 
S«^ose  that  the  manner  in  which  the  process  rate  is  affected  bgr 
pressiare  and  teaparature  during  this  "pressure-sensitire”  portion  of 
the  orerall  time  lag  can  be  expressed  as  a certain  function  f(P»T)* 

We  assume  that  for  the  giren  portl^s  of  propellants^  the  pressure- 
sensitive  portion  % of  the  time  lag  is  ended  (and  the  burned  gases 
are  evolved)  «hen  the  processes  have  reached  a certain  state  of  com- 
pletion, fills  can  be  expressed  more  concisely  by  the  relation 

\f(p,T)dt’  » constant  (1) 

mbere  p and  T are  functions  of  the  tine  t*  following  the  path  of  the 
propellant  portion  considered,  of  course^  if  the  temperature  is 
eonrelated  with  pressure  so  that  f » T(p)»  the  function  f(p,T)  depends 
only  on  the  Instantaneous  local  pres sura.  The  total  time  lag  will  then 
be  given  Iqr  and  the  space  lag  can  be  found  If  the  v^oci- 

ties  along  the  path  are  assumed  to  be  known. 

Let  us  now  look  at  some  of  the  results  obtained  for  the  ease  of 
low  frequency  instability.  This  is  the  simplest  of  all  cases  in  that 
the  pressure  developed  by  gas  generation  at  a certain  location  can  be 
considered  to  be  transmitted  instantaneously  to  all  points  of  the 
chamber,  since  the  wve-travel  time  is  much  amaller  than  the  period  of 
osciUatiims,  Thus  the  pressure  throughout  the  chamber  may  be  assumed 
to  oscillate  as  a whole,  tMch  leads  to  the  substantial  sixplificatlon 
that  the  space  delay  has  no  effect  on  the  results,  and  only'  a knoidedge 
of  the  time  lags  is  necessary.^  Furthemore,  knofd.edge  of  an  average 
value  of  the  time  lags  among  all  pc^tions  of  propellants  is  approxi/- 
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mately  sufficient  for  study  of  the  problem,  and  although  substantial 
deflations  In  time  lag  aay  exist,  their  effects  cm  be  sho%m  t®  be 
smll« 

In  the  first  treatment  of  this  kind  of  instability  by  von  Kasmjan 
and  his  group,  the  time  lag  was  supposed  to  be  a fixed  quantity,  in« 
sensitive  to  pressure  or  temperature  variations.  If  the  rate  at 
^eh  the  pressure  bulld-t^  follows  the  burning  rate  (depending  m 
the  balanee  between  chamber  volume  and  exhaust  area)  and  the  rate  at 
which  the  injection  rate  folloisB  the  pressure  variations  (determined 
hy  the  pi^ical  oimditions  of  the  feeding  system)  are  such  that  fcs^ 
the  given  time  delay  the  burning  rate  is  increased  the  pressure 
is  highsr  than  normal  and  vlce-veraa,  the  ccnditicais  for  \sstable 
eombusti<m  are  created,  theoretical  analyses  based  on  this  assuiip* 
tion,  were  published  by  Qunder  and  Ariant^,  Yaehter^,  and  Sunaerw 
field^,  the  analyses  were  made  poseible  by  a linearisation  of  the 
equations,  following  the  assuxption  that  the  fluctuations  of  all 
quantities  are  relatively  small. 

It  might  be  well  to  include  here  a faw  reenarks  on  this  assusp» 
tion,  which  is  the  one  that  makes  possible  general  theoretical 
dictions  «md  therefore  has  been  adopted  so  far  in  all  treatments  of 
the  problem  of  instability.  The  fact  that  in  actual  unstable  com- 
bustion the  aaplitiides  are  large  has  gwrierated  some  doubts  about  the 

1 •’Stability  of  Plow  la  a Rocket  Motor,**  by  D,F,  Qundsar  and  D.R, 

Priant,  Journal  of  j^pUed  Mechanics,  vol,  17,  no,  3,  1950* 

PP.  327-333 

2 "UlBQUssion  of  tha  Paper  of  Reference  !,••  by  M,  Taohter.  Jour- 

nal of  Applied  Hechaaics,  vol,  18,  no,  1,  19^,  pp,  Uh-llSj 

3 ”A  Theory  of  Unstable  Ccsabustioa  in  Liquid  Propellant  Rocket 

Systams,*’  by  M,  Susacer field  (presented  at  the  Heat  Transfer 
and  Fluid  Mechanics  Institute  at  Stanford  University  on  June 
20-22,  1951),  Journal  of  ths  American  Rocket  S5»eiety,  vol,  S?.® 
no*  5*  195^, 
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suitability  and  the  practical  seaniog  of  such  linearised  analysis*  These 
doubts  would  be  justified  if  the  theory  did  not  predict  unstable  condi«> 
tions,  since  in  this  case  non»linear  effects  coitld  conpletely  change  the 
results.  However^  since  the  linearised  treatment  actually  predicts  in» 
stability  in  certain  ranges,  it  is  certain  that  in  these  ruiges  the  non* 
linear  effects  cannot  substantially  change  the  results,  and  that  the  only 
radical  contribution  of  non*llnearity  would  be  the  appearance  of  addition* 
al  unstable  ranges  for  large  disturbances.  The  inclusion  of  no»*linear 
affects  would  in  this  case  produce  only  a refinexeat  of  the  linearised 
theory,  tdilch  in  itself  provides  a trustworthy  tool.  Of  course,  one  nust 
be  careful  in  applying  the  predictions  of  the  linearised  theory,  which 
apply  only  to  conditions  of  incipient  Instability,  to  check  experimental 
results  obtained  with  folly  developed  unstable  conditions* 

In  the  theoretical  analyses  xentloned  above,  it  was  shown  that  in* 
stability  could  appear  only  if  the  relative  pressure  drop  at  the  iajeo* 
tor  (ratio  between  Injector  pressure  drop  Ap  and  mean  chamber  pressure 
p)  is  below  0.5  « When  this  condition  Is  verified.  Instability  can  be 
generated  only  if  the  time  deli^  is  above  a certain  critical  value  idiieh 
depends  on  the  physical  conditions  of  the  feeding  system  and  of  the  cham* 
ber*  The  frequency  of  the  oscillations  under  conditions  of  iaoipiont  in* 
stability  was  also  determined,  and  could  vary  in  a wide  range  of  low  fre- 
quencies depending  on  the  forementioned  physical  conditions*  It  should 
be  noted  that  in  this  case  the  frequency  of  the  oscillatlous  is  not  re- 
lated to  any  natural  frequency  of  oscillation  of  the  systw,  but  is  de* 
termined  ty  a balance  between  time  lag  and  relaxation  tiEss  of  the  pro- 
cesses follonilag  the  offset  of  the  mass  balance  in  the  combust  ion  cham- 


ber and  of  the  dynamic  balance  of  the  feeding  3yat«s* 
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XT  noir  we  assune  that,  tho  time  lag  le  not  oonatant»  the  analysia 
becomes  somei^t  raox*©  involTod*  Considering  a pressvsre-seBsitlvo  por» 
tion  of  the  time  lag,  as  disonssed  previously,  and  exlstents®  of 
some  correlation  between  pressure  and  teinperature  oscillations,  the 
assumption  of  Equation  (1)  may  be  introduced  and  used  to  formulate  a 
reanalysis  of  ^he  problem#  The  form  of  the  functional  relation 
f(p,T)  is  required  in  order  to  pursue  the  analysis,  and  was  taken  to 
be  proportional  to  the  n^^  power  of  the  pressure,  idiere  n is  called 
the  pressure  interaction  index*^ 

The  fundamental  results  of  the  analysis  of  low> frequency  instability 
considering  a pressure-sensitive  portion  of  the  time  lag  are  as  foUowst 
If  n is  larger  than  0^5,  IngtabiUty  is  possible  as  a result  of  the 
interaction  between  pressure  oscillations  and  burning  rate  even  for  con- 
stant rate  of  injection*  In  order  to  obtain  this  so-called  intrinsic 
instability,  the  time  lag  must  again  be  larger  than  a certain  critical 
value  idiich  depends  on  n,  and  on  the  so-called  residence  time  Qg  of 
the  gases  in  the  chamber  (equal  to  the  ratio  between  the  mass  of  gases 
pr*esent  in  the  chamber  and  the  average  or  steady-state  flow  rate)*  As 
in  earlier  analyses,  it  is  observed  that  the  characteristic  isn^tability 
frequency  may  fall  in  a wide  range* 

If  now  both  time  lag  and  injection  rate  are  variable,  the  resulting 
conditions  are  more  difficult  to  discuss  in  general,  since  the  type  of 
feeding  system  has  an  important  effect  on  the  results.  However,  by  con- 
sidering two  typee  of  feed  systwas,  one  with  a constant  pressure  su^ly 
and  the  other  with  a constant  rate  supply,  it  is  possible  to  detersdae 
the  effects  on  stability  of  variations  in  such  quantities  as  line  length 

and  elasticity*  This  is  done  in  Reference  It,  and  diagrams  are  prescnS'ii.i 

it  ** Aspect 8 of  Combustion  Stability  in  Llqiiid  P»:*opel3j®t  Socket 

by  Lj  Crocco,  Journal  of  tho  American  Rocket  ^ciety,  Nov,  1951  and 
Jan»Feb  1952*  ErS'ffiOTl® 
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iMoh  giv«  the  values  of  the  critical  tiiae  lag  above  ishich  cs»bustion  is 
unstable  and  the  corresponding  period  of  the  oscillation  whicb.  will  exist 
unde**  given  feed  line  lengths,  injector  pressure  drops,  and  presmire  inter- 
action indices*  These  uiagrans  Indicate  the  isportant  fact  that  if  the 
rocket  is  intrinsically  unstable,  it  is  possible  only  in  exceptional  cases 
to  obtain  stability  through  changes  in  the  feed  systm,  and,  in  general, 
cmibustion  will  always  remain  unstable.  The  diagrams  can  be  made  to  apply 
to  nu>re  explicated  feeding  aystaas,  involving  combinations  of  non-aer© 
line  length  and  non-aero  elasticttys  They  may  be  lised  bo^  far  aonopro^ 
pellant  cxbustion  and  for  bipropellant  systems  in  which  certain  relations 
are  satiefied  by  the  gexetry  of  the  two  feeding  systems*  When  these  cer- 
tain relations  are  not  verified  in  the  case  of  bipropellant  systeus  the 
diacussioa  becooes  more  involved,  and  it  is  difficult  to  present  general 
diagrams  or  predictions. 

Sumaarizing  the  discussion  of  low-frequency  instability,  if  the  con- 
bustion  is  not  intrinsically  unstable,  simple  variations  of  the  physical 
parameters  of  the  feed  systxs  and  of  the  chamber  nay  cure  instability* 

If,  however,  the  cxbustion  is  intrinsically  unstable,  changes  in  the 
feeding  system  are  generally  ineffective,  and  the  cure  should  cxe  through 
an  increase  in  residence  time  «r  a decrease  in  time  lag* 

For  more  explicated  systems  with  conplicated  automatic  controls  there 
My  be  interference  between  the  controls  and  the  combustion*  It  seems, 
however,  that  generally  tte  frequencies  of  the  servosystems  and  of  the  un- 
stable cxbustlx  are  in  differraat  ranches,  which  case  Interfarence  Is 
n^ligible*  A pure  instability  of  servreontrols  is  possible  but  this  must 

not  be  considered  as  a typo  of  combustion  inatobilityp  Hi®  possibl®  inter- 

< 

ferenee  of  servocontrols  with  ecmbuoti©B  «ias  suggestcfl  tf^  Toica  mi 
5 ’’Servo  Stabilization  of  Roekot  Ilcfeoy  Oscillations,”  b/  H.S® 

Journal  of  the  American  Boshot  Soeiety,  «ov-Dac  1952 
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Ifarbl®  ft  possibility  of  stsbiUsation  of  unstable  combustion  vken  other 
meant  are  not  effect  ire  or  possible* 

A.  different  method  of  analysis  is  required  for  the  study  of  high 
frequency  instabilityf  In  which  acoustic  resonance  is  the  coordinating 
factor*  Of  the  several  different  possible  types  of  interaction  idiich 
may  produce  instability,  the  one  which  is  believed  to  be  the  most  im» 
portent  for  good«>perforBanoe  rockets  is  the  combust ion»pres8ure»tetsper> 
aturs  interaction,  and  only  this  interactlcm  has  been  analysed*  The 
l^esent  discussion  is  also  limited  to  the  case  of  longitudinal  instabil- 
ity, so  that  a knoidedge  of  the  transversal  distribution  of  combustion 
is  unessential*  Therefore,  a longltudiricl  distribution  of  time  lags  and 
space  delays  is  all  that  is  needed  for  the  study  of  the  problem,  to- 
gether with  a proper  definition  for  the  interaction  mechanism* 

The  distribution  of  space  lags  is  determined  when  the  eonibustion 
distribution  along  the  chaniieF  is  prescribed*  Two  estremc  types  have 
so  far  been  analysed:  a concentrated  combustion  f^ont  corresponding  to 
a single  value  of  the  space  delay  fcr  all  of  the  propellant,  and  a 
unlftmly  distributed  combustion  in  which  the  velocity  of  the  gasos  in- 
orsases  linearly  with  the  distance  from  the  injector  end*  Boring  the 
preseure  insensitive  part  of  the  time  lag,  :ditleh  is  unessential  for  the 
following  considerations,  the  velocity  of  the  propellants  is  sloped  down 
from  the  injection  velocity  until  it  reaches  the  velocity  of  the  main 
flow  of  gases,  and  during  the  pressure  sensitive  part  of  the  time  lag 
the  average  motion  of  the  gas  is  followed  by  the  propellants*  In  the 
eoiieentrated  combustion  case,  no  velocity  of  the  gases  exists  tqtstresm 
of  th®  combustion  front*  Therefore,  the  displacement  of  the  propellants 

during  the  pressure  sensitive  lasrt  of  the  time  Isg  ia  aero  gad  the  values 
"S^o  Stabilisation  of  Low  Frequency  Oscillatiowa  in  a Bipropellant 
Rocket  Mot  or,  ^ ly  F.E®  Marble,  Jovurnal  of  the  jimeriGan  Rocket  Society, 
March-April  1953  RSSi>jSIGTED 
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ftf  th«  pressuraa  under  the  iat%T»l  aign  have  to  be  eoaqjuted  at  the  loca« 
tlon  of  the  front » In  the  c&ae  of  distributed  combuntion  the  velocitjr 
is  t)on-aero>  and  the  ohange  in  location  during  the  integration  has  to  be 
taken  into  account.  Kowevtr,  the  sinpliiying  assunptlon  has  been  »ade 
that  this  displacement  is  small  as  coieoared  sith  the  total  space  lag. 

The  probl«s  is  treated  by  writing  doun  the  equations  of  the  gas 
notion  for  small  non- steady  perturbations  with  the  prescribed  sources  of 
bursoed  gases^  by  supposing  thereafter  an  esqponential  time  dependence  with 
genoral  complex  exponent  and  searching  for  the  neutral  ccmditions.  ^ere 
the  exponent  is  a pure  imaginary  and  the  system  satisfies  the  prescribed 
boundary  cmditions.  The  systma  is  further  sinplified  by  supposing  that 
the  Kach  numbers  of  the  gas  flow  are  sufficiently  small  throughout  the 
chaad>er  so  that,  their  squares  can  be  neglected  with  respect  to  unity. 

The  boundary  conditions  to  be  applied  are  the  followings  at  the  in- 
jector end  the  perturbation  of  the  velocity  must  always  be  sero^  and  at 
the  nozsle  end  a certain  ooi^lex  ratio  between  the  perturbations  of  vel- 
ocity and  pressure  or  density  (analogous  to  the  ordinary  acoustic 
dance)  Is  prescribed.  This  complex  ratio  has  been  determined'*^  for  a 
particular  type  of  nosale  with  linear  velocity  distribution  in  the  con- 
verging partj  as  a function  of  the  Kach  number  at  the  nossle  entrance. 

TIte  phase  difference  between  the  perturbations  starts  fron  sero  at  sero 
frequency^  Increases  to  a snaximum  as  frequency  is  increasad,  and  then 
decreases  again  to  aero  for  very  high  frequencies.  The  ratio  of  aapli- 
tudes  of  velocity  and  density  increases  steadily  fr«m  the  value  -1 

"T^ 
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7 ’’The  Transfer  Function  of  Rocket  Moazles*’,  by  H.S.  Tsien,  Journal  of 

the  Merioan  Rocket  Societyf  May-Jtine  19^2 

8 •’Siapercritical  Gaseous  Discharge  with  High  Frequency  OsciUaticast" 

by  L«  Grocooji  VIXI  Congress  of  Applied  Mechanics}  (to  be  published 
in  *L'Aerotecnica”) 
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THEORKTICAL  SrUDlES  (Gont5d«) 

idiich  i«  characteristic  of  very  low  frequencies  (and  corresponds  to  the 
condition  of  constant  Mach  nuxber)  to  & much  larger  value  characteristic 
of  infinite  frequencies.  Two  casea  have  been  considereds  one  for  an 
ideal  one>dinenaional  nozxle  of  sero  lengthy  and  the  other  for  an  ideal 
nosale  in  which  the  length  of  the  converging  part  of  the  nossle  is  about 
one^third  the  length  of  the  chamber. 

HtuB  results  of  the  analysis  for  the  concentrated  combustion  ease 
are  |n*esented  in  detail  in  Appendix  A.  These  results  can  be  sumnarised 
as  followss 

When  the  combustion  front  is  located  close  to  a pressure  node  in 
the  ohanher  the  corresponding  mode  of  oscillation  ic  alMiys  stable^ 
but  instability  appears  when  the  combustion  front  moves  suffieimtly  far 
away  from  the  nodes.  The  most  dangerous  situation  for  all  laodes  occurs 
idien  the  combustion  is  concentrated  at  the  injector  end  or  at  the  nozsle 
end. 

For  one  particular  value  of  the  interaction  index  use  of  the 
sero-length  nozzle  boundary  condition  indicates  that  instability  will 
exist  over  nearly  all  possible  values  of  time  lag  and  coabustion  loea- 
tioQj  if  modes  of  oscillation  higher  than  the  fundamental,  are  taken  in- 
to account.  However,  this  rathm*  discouraging  situation  ia  alleviated 
when  the  more  realistic  nozsle  boundary  condition  is  applied,  since  all 
nodes  except  the  fundanental  become  stable,  and  the  regions  of  instabil- 
ity of  the  fundanentaX  mode  itself  are  greatly  reduced. 

If  the  interaction  index  a is  increased,  the  stability  conditicns 
become  less  favorable,  and  it  is  cleai>  th&t  n has  an  effect  somefdiat 
similar  to  that  described  in  the  low  Xre<queacy  case.  Indeed,  S.t  can 
actually  be  shown  that  high'-frequaney  instability  can  o<ily  b«  present  If 
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III.  THEORECICAI.  SfUDIES  (Cont*d,) 
n is  linger  than  a value  voey  close  to  0*5  * 

. The  analysis  for  the  case  of  distributed  combustion  is  presented  in 
detail  in  Appendix  B«  The  fundamental  x^sult  of  this  treatment  is  that 
no  matter  hov  long  the  noaale  is  taken  to  be«  the  values  of  the  inter» 
action  index  n for  tdiich  instability  may  appear  are  well  abov@  the  rashes 
to  be  esipeoted^  and  therefore  the  combustion  is  always  stable^ 

It  is  interesting  to  ccmipare  some  of  the  mentioned  results  with  the  ex- 
cellent  experimental  work  by  SerBsan  and  Cheney  of  General  Slectrie^*  In  this 
work  the  ccxnbustlon  has  been  viewed  through  a slot  window  in  a eylin<h:‘ieal  . 
eonbustion  chaiifeer  and  the  wave  motion  actually  observed*  Coeparisons  with 
pressure  oscillations  gave  a full  check*  With  a short  nossls  the  instabil- 
ity could  easily  be  obtained  and  gave  rise  to  high  amplitude  oscillations 
that  degsnerated  into  a shock  wave  oscillatory  processi  however^  this  shock 
type  of  instability  was  always  preceded  by  & wave  type  of  increasing  ant* 
pUtude*  When  the  nossle  was  replaced  by  a longer  one  with  the  same  throat 
area^  the  stability  conditions  seemed  to  be  considerably  improved  and  the 
shock  type  was  never  observed^  and  at  the  same  time»  the  frequency  de- 
creased slightly.  This,  and  many  other  experimental  results,  are  in  sub- 
stantial agreement  with,  the  predict5.ons  of  the  theory  and  seen  to  indicate 

f 

that  the  assumptions  are  substantially  correct* 

Much  theoretical  work  remains  to  be  done,  in  order  to  extend  the 
evaluationa  in  the  high  frequency  case  to  more  general  combustion  dis- 
tributions and  less  i 'rictlve  aesui^tions  about  the  time  lags,  and 
especially  to  develop  an  ar»alogou»  treatment  for  the  case  of  transversal 
modes,  iMch  se«»a8  to  become  of  coiiparable  and  even  greater  importer "le 
than  the  longitudinal  type  idles  the  thr»ust  of  the  rockets  increases  and 

% 

the  transversal  dimensions  become  equal  to  or  larger  than  the  longitudinal 

^ iwii*.wwmrmirfwiwwraii— LJL  W'K'iiu.ai  ' vem*-.c-<waiiwWwsirTaaij!ai  iisj.aj.  TmTfriT-imww-niitnwtT—T  r— ,mrr rirartn  — r-n-;  -■* -vir --l■l£^-^3au-  - 

9 ^Combustion  Studies  in  Rocket  Kotora,**  by  Kr.  fiarriea  and  S,h,  Oheneyg 
A»®tlcan  Rocket  Society  Meeting,  Dee«  195^,  How  York  SSSTHlCTEIi 
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tHEOaETIGAL  3TUDISS  (ContM*) 

A syatsnfttic  oxpttrlMcntal  )*«se6X>ch  in  pKr'nIlal  vith  th«  th«K3(retieaI 
Is  now  being  oonduoted  under  the  tenus  of  the  subject  contr«cb»  mth  the 
definite  purpose  of  detemlnlng  the  on^irlcel  l&we  of  the  tine  lag  and 
eheoIdUstg  the  predict  Iona  of  the  theo(ry»  It  is  only  vith  a dwd  progran 
of  this  kind  that  man  can  hope  to  gain  full  understanding  of  this  coG^li- 
eated  phenomenon  and  of  the  ways  to  prerent  it* 
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17.  AFFARUCS 
A*  Faoillties 

The  rocket  test  cell  and  control  roois  are  now  cos^>let«|  and  InstaUso 
tlon  of  the  nonopropellaht  feed  system,  thrust  starid^  and  control  system 
is  under  way.  Inst nimentat  ion  wiring  has  been  run  between  the  rocket  con* 
trol  room  and  the  central  recording  rom  in  •’B"  Building.  An  lastruiaenta- 
tion  control  console  and  six  vertical  instrument  display  racks  have  beeii 
built  and  installed  in  the  central  recording  rocm,  and  Installation  of 
instrisMntSy  control  coiq}on€iats  and  wiring  is  progressing.  It  is  eiqpected 
that  this  equipment  will  be  sufficiently  cos^lets  to  be  used  on  fe®d  system 
flow  calibration  tests  without  introducing  additional  delays* 

The  hydraulic  and  pneumatic  conponmt  test  panels^  although  lacking 
SOM  components  due  to  extremely  late  deliveriaa,  are  in  current  use  for 
oalibration  and  operational  testing  of  various  feed  system  components. 

It  is  expected  that  application  of  this  facility  will  be  sxtsnded  enor- 
mously idien  the  delayed  c<«poneats  arrive  and  can  be  put  to  use. 

B.  Constant  Rate  Honopropellant  Feed  ^tma 

Installation  of  piping  and  standard  feed  qrstem  c<xaponents  on  the 
veat  wall  of  the  rocket  test  cell  has  been  conpleted>  and  is  shown  in 
Figure  1.  All  instrumentation  and  pneumatic  control  lines  have  been  run 
through  a convenient  bulkhead  for  ease  in  pi'es sure-check  and  calibration 
proceduresy  as  indicated  in  the  figure.  The  concrete  thrust  stand  founda- 
tion for  Bonopropellant  tests  has  been  poured^  and  the  thrust  stand  itself 
has  been  constructed  and  mounted  in  the  test  cell.  Delivery  of  components 
such  as  the  pulsing  unity  propellant  valve  and  emergency  propellant  valve 
actuators,  and  rocket  motor,  all  of  which  are  to  be  mounted  on  the  thrust 
stand,  has  been  delayed,  but  these  components  are  all  expected  to  arrive 
shortly.  The  major  obstacle  encountered,  however,  has  been  the  extensive 
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apparatus  (ConVd.) 

B.  ConstMit  R*t«  Mon<^ropell*at  P«®d  ^st«i  (Oont'de) 
delay  in  delivery  of  the  nonopropellant  tank,  and  this  has  sharply  c«« 
tailed  uork  on  both  tlie  gae  pres eur icing  line  and  tho  llqidd  propellant 
line«  It  is  possible  that«  if  any  further  delays  are  encouiistered  in 
delivery  of  this  itea«  a temporary  taidc  .say  bo  rigged  as  an  interim  meas- 
ure* Possible  delays  in  delivery  of  the  pulsing  unit  have  been  foreseen^ 
and  consequently  an  interim  unit«  consisting  of  a butterfly  mounted  in  a 
shunt  line»  has  been  designed  and  is  nearing  c<M»pletlon  in  the  shop* 

Construction  of  a heat  exchanger  to  act  as  a source  of  heat  in 
maintaining  uniform  propellant  temperature  in  the  tail..:  Is  in  progress^* 

This  heat  exchanger  uiU  be  moxmted  directly  under  the  i'eed-sysima  com- 
ponent panel*  The  control  cosisole  and  coim>anion  panels  aiilch  mill  ssount 
operational  gages  and  meters  are  to  be  delivered  in  December*  Comp^enta 
for  these  panels  are  arriving  dally,  and  it  is  expected  that  tl»y  uili 
all  be  available  ahen  installatlm  of  the  monopropellant  feed  t^staai  is 
complete*  A schematic  of  the  monopropellant  control  system  is  shoan  in 
Figure  2* 

Calibration  tests  of  the  feed  system  cosponents  were  initiated  and 
operati«sial  checks  of  various  standard  coaponents  are  in  progress*  Static 
testing  of  the  first  cavitating  venturi,  which  is  expected  to  provide  a 
constant  rate  of  propellant  floif  for  a ciiamber  pressure  of  300  psi,  has 
been  cwspletad*  The  venturi  performed  according  to  design  specifications, 
passing  a constant  flow  of  1*5  Ib/sec  downstream  preasurea  ranging 
ircM  approxiamtsly  100  to  h50  psi©  Dynajidc  calibration  of  ths  venturi 
awaits  ai*rival  of  the  pulsing  unit  (or  interim  ptasing  unit)  and  coaqjlc-^ 
tion  of  calibration  of  the  pressure  pickups* 

A servicing  trailer  to  facilitate  cthylonfs  orid®  loading  operatiwia 
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IV,  APPAliJgUS  (C0Kt*d.) 

Bs  Goaetant  »at©  Honopropslianfc  Feed  syutOM  (Contod,) 
has  been  designed  and  constr-wcted*  The  trailer  is  fitted  with  « swinging 
bocn  and  winch  to  aid  in  handling  drums  of  ojdde^  and  a service  piping 
system  for  transferring  osdde  from  the  drums  to  the  monopropellant  tank 
without  air  contact  is  ready  for  installation* 

Co  Monopropellant  Rocket  Motor 

The  rocket  motor  shown  in  the  first  quarterly  report  is  now  being 
manufactured^  and  should  be  delivered  in  December*  The  eight  Inject^s 
are  nearing  conqileblony  and  although  severe  welding  difficulties  nay 
alter  the  anticipated  delivery  date«  it  is  expected  that  they  will  arrive 
well  before  the  motor  itself  is  delivered* 

A high>power«  low  voltage  ignition  unit  has  been  delivered  free*  the 
Scintilla  Magneto  Olvisioa  of  Bendix*  This  unit,  the  TIH»10^  has  been 
used  successfully  under  extremely  unfavorable  oondxtions  in  turbojet  en» 
gineSf  and  should  siqiply  satisfactory  ignition  in  tho  monopropellant 
motor*  Scintilla  has  also  supplied  two  adjustable  e:^rimental  igniter 
plugs  for  use  with  the  TlS-10  tmit* 

The  oxygen  system  for  starting  and  the  inerWgas  purge  system  have 
been  desigiied  and  all  parts  preeured*  The  inert~gas  purge  will  start 
automatically  ^en  the  tank  vent  is  closed  and  the  emergency  propellant 
val.v»  opens*  It  stops  when  the  propellant  valve  openss  and  starts  again 
automatically  at  the  end  of  a run  wiien  the  propellant  valve  is  closed* 
OJQrgen  How  is  controlled  by  a magnetic  valve  which  is  coupled  directly 
to  the  electrical  ignition  system  so  that  spark  and  oxygen  are  tusmod  on 
manually  together  and  are  both  shut  «»ff  automatically  by  tm  adjustablft 
time^delay  relay  aotuated  a chenbor  pressuf  o s?}itch  (eae  Fig*  2).^ 
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Vf*  APPMAm  ( COttliM.) 

B*  XnatnMentatdon 

!Ihtt  ma^cn*  affort  on  iastrumantation  dm'ing  th«  period  eo'vorad  l>7  thia 
report  ha»  been  tn  the  derelopment  and  testing  of  the  Prinoet^»MIf  pre&s* 
are  pioSoi^  and  in  d^manle  floieseter  resoareh.  Conetruction  nod  wiring  of 
an  Inetrunent  control  console  and  display  racks  and  operational  cheeks  of 
conaercial  eqtiipaent  received  during  this  period  have  progressed  sati»» 
factoriij;  and  most  of  this  typo  of  work  is  expected  to  be  completed  short» 
ly«  The  Majority  of  ctmnnerclsl  equipsent  ordered  has  be«a  deliv»red«  not» 
able  exceptions  being  the  sonio  analyser^  the  Id  nass  flowsmBtar^  the 
I^ttelman  electrssaagnetic  flotmeter*  and  the  specially-designed  four- 
chuuMl  filter}  all  of  which  are  well  overdue^  but  are  eiqoected  liy  January^ 
1953. 

(1)  Sensing  ElsHcnts 

Three  differential  pressure  pieki;q>s«  one  absolute  pressure  pickupj» 
and  tmip  duRmy  pickup  (identical  to  the  absolute  pressure  unit  except  that 
the  strain  tube  element  is  lacking}  were  delivered  by  Control  Engineering 
Corporation.  Figure  3 shows  a photograph  of  the  differential  pressure 
unity  which  is  the  priate  feature  of  Instrumentation  the  subject  con- 
tract y and  the  sketdh  of  Figure  U illustrates  its  operation. 

The  transducer  is  designed  for  flush-mount  log  in  the  rocket  motor 
combustion  chamber.  The  two  catenaiiqr  diaphragms  shown  ia  Figtar©  are 
rigidly  fastened  together  at  the  Inflection  points  by  a metal  ringy 
and  provision  is  made  to  siqjply  a flow  of  cooling  iat©r  between  the  dia- 
phragas.  A strain  tube  is  mounted  on  the  upper  diapbragMy  and  strain 
w3res  bonded  to  the  tube  are  deformed  when  the  diaphragm  flexes  under 
pressure.  These  wires  form  one  side  of  a liliea’fcstone  bridge,  so  that  if 
the  input  voltage  across  the  bridge  ie  fixed  and  the  bridge  balanced  at 
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IV.  AP?M/gUS  (Cont«d.) 

D.  Instrtunentation  (Ccmt*d.) 

(1)  Sensing  Elenants  (Oont*d.) 

some  pressure*  changes  in  resistance  of  the  strain  wires  (vhen  they  are 
defomed  as  a result  of  changing  the  pressure  on  the  diaphragm)  cause  the 
bridge  to  become  unbalanced,  and  the  unit  is  designed  so  that  the  un* 
balanced  voltage  will  be  directly  propcotional  to  the  magnitu^  of  the 
pressure  change.  In  order  to  obtain  very  high  sensitivity  to  smaU 
pressure  changes  at  high  pressure  levels,  the  chamber  abpve  the  top 
diaphragm  is  sealed,  and  may  be  pressurised.  This  maris  that  if  it  is 
desired  to  measure  small  pressure  fluctuations  in  a chamber  operating 
at,  600  psi,  the  "back  pressure"  on  the  diaphragm  is  set  at  some 
value  in  the  neighborhood  of  6OO  psi«  fhen,  the  strain  gage  mill  rsad 
differences  in  pressure  from  this  reference  value,  and  c<aaaeqaent3y  a 
very  sensitive  strain  gage  Kiay  be  uiedi  one  that  has  a ran^e  of,  say, 

A 

- 100  psi.  If  no  "back  pressure"  could  be  applied,  it  would  ba  neoess* 
ary  to  use  a strain  gage  with  a range  of  700  psi,  and  hence  much  lower 
sensitivity  to  ennall  fluctuati<ms  would  be  obtainable.  One  difficulty 
introduced  by  the  use  of  highly  sensitive  strain  gages  is  the  necessity 
for  overpressure  protection,  and  this  is  included  in  the  transducer  in 
the  form  of  a rigid  metal  ring  fastened  to  the  main  pickup  body  between 
the  diaphragms,  thus  limiting  deformations  to  acceptable  amounts  in  both 
directions. 

In  addition  to  the  differwxtial  transducer  described  above,  an  ab> 
solute  pressure  pickup  of  the  same  basic  douole-diaphragm  design  wag 
delivered  for  use  in  the  measurement  of  combustion  chamber  pressure 
level,  starting  transients,  etc*  Unfortunately,  this  pick\^  was  des- 
troyed during  operational  testing,  as  will  be  described  later,  and  a 
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APPARjItrUS  (Cont'd.) 

D,  Xnstrunentatltm  (Cont'd.) 

(1)  Sensing  Elemsnts  (Cont*d«) 

replacement  will  be  ordered  for  use  in  the  monopropellaiit  aiotor* 

Qperatioiu;!  testing  and  oalibriitlonal  tests  of  the  PrincetonoHIf 
design  hare  been  bc^un»  and  results  of  the  initial  testing  are  reo 
ported  in  a later  section* 

Kuch  research  has  been  done  during  the  period  covered  ^ this 
report  on  the  problem  of  muasuretnent  of  transient  floH  phenomena* 

As  reported  earlier,  three  basic  flowte1;er  types  have  been  invest!* 
gatedt  the  angular-momentum  i^ss  flounet^  designed  by  Dr*  I.T*  Id 
of  the  ICDf  Instrumentation  Laboratory,  the  electrcHaagnetle  floisaeter, 
and  the  hot-vire  anemometer  modified  for  use  with  liquids* 

the  basic  princlpls  of  operation  of  Dr.  Li*s  mass  flowmeter  iw 
indicated  in  the  diagram  of  Blgure  5*  The  fliiid  flofws  through  a D* 
tube  idiioh  is  rotated  at  constant  speed*  A strain  tube,  to  idideh  a 
strain  gage  is  bonded,  ooag):rlse8  a section  of  the  rotating  member, 
as  shown  in  Figure  5>  As  the  mass  flom  varies,  the  angular  momentum 
of  the  fluid  in  the  U-tube  changes,  and  hence  the  torque  required  to 
maintain  its  rotation  at  constant  speed  will  change*  Since  the  strain 
tube  transmits  torque  frt»  the  motor  to  the  U-tube,  its  deformation 
be  proportional  to  the  torque  and  hence  to  the  angular  momentum  of  the 
fluid  In  the  U-tube*  Thu%th*  strain  gage  reading  is  directly  prqpor- 
tlwMl  to  the  angular  nomcntum,  or,  since  the  geometry  ^he  rystem  is 
fixed,  to  the  mass  flow*  The  frequency  response  of  the  flowmeter  de- 
pends primarily  on  the  natural  frequency  of  the  mechanicai  end  hydraulic 
portions  of  the  fluid  passages,  since  the  natural  fTequ^aoy  of  the  elec- 
trical elements  of  the  syst«n  will  be  quite  high*  Preliminary  checks  by 
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SCHEMATIC:  LI  MASS  FLOWMETER 
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apparatus  (Oont’d.) 

U*  Xnstruffientation  (ContM.) 

(1)  Sensing  Klewenta  (C<ait»d») 

the  dosigner  indicate  a natural  firequency  tor  the  unit  of  approximately 
800  cycles  per  second^  which  means  that  the  rei^onae  should  be  linear 
to  at  least  200  cpa  and  may  possibly  be  good  up  to  300  cps«  This  flow^. 
meter  is  to  be  delivered  in  January,  1953*  at  which  time  it  will  be  cal- 
ibrated for  frequency  response,  phase  lag,  and  average  flow  rate,  and 
will  compered  with  other  dynamic  flowmeters* 

'Tho  designs  of  electromagnetic  floweieter  will  be  used  im  flow  cal- 
ibration tests,  the  Hittelman  meter,  designed  and  built  by  Eugene  Mittel- 
man,  Inc*,  and  the  Arnold  meter,  designed  by  Br*  James  Arnold  of  the 
Physical  Science  Laboratory,  New  Mexico  College  for  Agriculture  and  Mech- 
anical Arts,  las  Cruces,  New  Msxioo*  Both  meters  are  similar  in  operat- 
ing principle,  but  have  several  differences  in  construction  and  CM^on- 
ent  design*  The  Arnold  meter  is  to  be  delivered  shortly,  on  loan  l^om 
the  Applied  Physics  Laboratory  at  Johns  Hopkins  University,  and  the 
Klttelman  meter  is  expected  In  January,  on  loan  front  the  Naval  Air 
Rocket  Test  station  at  Lake  Doimark,  N.J. 

The  operating  principle  of  the  electromagnetic  flowater  may  be 
outlined  as  follows^  A magnetic  field  is  set  tq>  around  a non-conduct- 
ing tube  through  which  the  fluid  flews*  If  the  fluid  is  a conductor 
of  electricity,  a current  will  be  induced  in  it  as  a result  of  the  mag- 
netic field,  and  the  magnitude  of  this  current  will  be  proportional  tc 
the  velocity  of  the  fluid  flow,  according  to  Faraday's  principle  of 
electromagnetic  induction.  This  current  is  then  transmitted  through 
electrodes  set  Into  the  nonconducting  tube  to  an  electronic  amplifier 
and  recording  instruments* 
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IV,  APPAHATUS  (Coat’d,) 

X),  Instmmentatlon  (Cont’d,) 

(1)  Sensing  Elements  (Cont’d^ 

m order  to  obtsin  the  highest  possible  frequency  response,  a constant 
ajagnetic  field  (snsdntained  by  direct  current)  would  be  most  advftntageou8| 
however,  this  type  of  field  introduces  a asajor  problem  known  as  polariaa« 
tioa,  A conducting  liquid  will  electroiysq  idien  an  electric  current  is 
passed  through  it,  and  if  this  current  is  always  in  the  samo  direction, 
as  would  be  induced  by  a constant  magnetic  field,  electrolysis  of  the  liq» 
uid  will  cause  gas  to  form  on  the  electrodes.  When  the  nonconducting  gas 
coats  the  electrodes  completely,  the  instrument  is  useless,  since  the  slg* 
nal  can  no  longer  be  maintained.  For  this  reascm,  it  has  been  necessary 
to  use  an  alternating  magnetic  field,  ihen,  however,  the  frequency  res<» 
ponse  of  the  instrument  is  limited,  since  Instead  of  measuring  the  flow 
continuously,  ws  are  obtaining  a number  of  discrete  readings  correspond- 
ing to  the  pulses  of  the  magnetic  field.  For  example,  suppose  the  flow 
were  oscillating  sinusoidally  at  200  cycles  per  second.  Let  us  say  that 
ten  points  on  each  individual  sine  wave  cycle  are  necessary  in  order  to 
delineate  the  waveform  of  the  flow.  Then  it  is  necessary  to  oscillate  the 
magnetic  field  at  2,000  cycles  per  second,  and,  since  electromagnets  poss- 
ess large  values  of  inductance,  the  problem  of  maintaining  linearity  can 
become  rather  severe. 

The  Mlttelman  meter  uses  a magnet  frequency  of  1,000  eps,  and  hence 
can  be  used  to  measure  fluid  flows  oscillating  at  100  eps.  The  Arnold 
flowmeter  is  claimed  to  have  a magnet  frequency  of  000  cps,  and  if 
forthcoming  te;.w«  Indicate  satisfactory  operation  of  that  frequency,  this 
meter  will  be  useful  at  flow  oscillations  of  500  cps,  . 

Sense  research  has  been  conducted  at  Princeton  on  the  possibility  of 
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I¥*  APPAaAT'03  (Coat’d®) 

B©  lnetsmm%&t±m  (Coat’d.) 

(1)  Sensing  £l€8se'it8  (Cont*d«) 

using  direct  eusrent  excitation  of  tho  rsagnet,  and  cancallisig  out  th© 

I 

DsG.  coaaponent  of  the  induced  Gui*r«nt  in  the  oscillatijig  liquid  flow 
by  applying  a reverse  D.C.  potential  across  the  electrode®®  !?hi8  taould 
prevent  polarisation  fross  occurring*  and  at  the  saras  tirae  w&ald  not  dis«* 
turb  the  AeC®  portion  of  the  induced  current  ^dsich  sssasur®®  th®  oscill- 
ating coB^onent  of  the  liquid  flow.  This  would  peraiit  meaaureKmt  @f 
very  high  frequcnciss*  sines  tbs  only  liwitation  i®  the  natural  fr®- 
quenoy  of  the  electronic  system  ooaponents®  Hiafozi;,unt8tely*  the  press 
of  other  work  ims  prevented  any  extensive  development  of  this  schme* 
and  it  has  been  put  off  until  more  time  or  manpower  becomes  available* 
The  hot-wire  type  flowmeter  is  nearly  ready  for  use  in  calihra» 
tlon  of  the  monopropellant  feed  system.  This  instrument  consists  of 
a fine  wire  wrapped  on  a non-conducting  support*  which  is  placed  dir» 
ectly  into  the  stream  of  fluid©  A constant  electric  curr^t  Is  passed 
through  the  wire  to  heat  it*  and  a isheatstone  bridge  is  used  to  measure 
the  reslstarice  of  the  wire.  Now*  as  the  flow  past  the  hot  wire  changes* 
the  heat  transfer  £vcm  the  wire  varies  accordingly  * and  hence  it® 
temperature  wilx  change.  Thin  alters  the  resistance  of  the  mH 
thus*  by  recording  the  output  of  the  Wheatstone  bridge,  a asasuro  of 
the  flow  variations  can  be  obtained*  IJje  major  failing  of  this  type 
©f  flowjaater  is  that  the  heat  transfer  characteristics  of  tho  wir©  arc 
nonlinear  and  subject  to  Euny  other  ©ffects  bssidos  flow  ratQ|  hcaco 
the  aaaaored  aiaplitudec  of  flow  oscillations  are  not  at  all  accus*;'.to. 
?b,Q  hot  wlra*  hot?3v©r*  ia  entirely  Batiai'cjcto;>ry  ao  a flow  phaceicoio?* 

aiase  its  frequoncy  respoaoa  vjith  a sufficiently  €541'j  tdro  ip.iopocd  ia 
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IV,  APPARATUS  (Cont9d.) 

D,  Instrottoabation  (Oont'do) 

(1)  Sensing  Slenents  (Cont^d,) 

I 

liquid  is  highj)  and  tho  nonlinsarity  does  not  affect  phase  measurmentsa 
Its  principal  application  in  connection  with  the  subject  contract  will 
be  for  calibrating  the  phase  difference  bet^leen  actual  flow  rate  and 
dTnamic  pressure  drop  across  the  injector  orifices  as  measured  bv  two 
pressure  transducers*  It  will  also  be  used  as  a check  on  the  other  types 
of  flowmeters  described  above. 

The  nuijor  problem  encountered  in  the  measurement  of  transient  flows 
is  the  fact  that  no  accurate  standard  of  reference  is  available  to  use 
as  a calibrating  datum,  A flowaeter  designed  to  measure  steady  flows  can 
easily  be  calibrated  collecting  and  weighing  the  liquid  which  passes 
through  the  meter  in  a given  period  of  time,  but  no  known  msthod  exists 
for  determining  accurately  the  transient  wave-form  of  an  oscillating  flow, 
' In  view  of  this  lack  of  a calibration  reference,  the  method  to  be  used  in 
connection  with  the  present  contract  will  consist  of  cceg>arlng  tiM  read- 
ings of  all  four  dynsmiic  flowmeters  (the  Li  type,  tho  Arnold  meter,  the 
Kittelman  device,  and  the  hot  wire)  at  flow  oscillation  frequencies  start- 
ing from  aero  (when  all  four  should  read  the  same)  and  increasing  to  th® 

t 

point  at  %tolch  the  readings  all  diverge.  It  will  be  assumed  that  the 
frequency  at  which  this  divergence  occurs  is  the  maximum  response  lAilch 
can  be  accurately  recorded,  although  it  is  certainly  possible  that  om 
of  the  meters  may  be  usable  at  a high«*  frequency.  Of  course,  if  it 
could  be  aesumed  that  the  flow  pulsing  device  were  delivering  a {mown 
i«ve»shapa  of  known  as^litude,  an  absolute  reference  could  be  had,  but 
there  is  no  way  of  knowing  whether  or  not  this  is  actually  tte  ^ ase. 

Hence,  in  order  to  be  conservative,  it  is  necessary  to  operate  only 
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XV.  APPARATUS  (Cont*d.) 

Do  Imrtruffiisntation  (Cont<d«) 

(1)  Sensing  EXassnts  (Ocmt*d«) 

to  frequencies  for  tMch  at  least  t%;o  of  the  flometer  readings  coincide; 
and  either  one  of  these  t«o  flometers  maqr  then  be  used  in  sucooeding 
tests* 

It  is  hoped  that  it  nill  not  be  necessary  to  Include  a d^nanle 
flom»ter  in  the  actual  rocket  aotor  tests;  but  that  it  uill  be  poss» 
ible  to  measure  oscillating  flous  recording  instantaneous  pressure 
drop  across  the  injector  orlfiees*  In  order  to  Terify  the  theoretically 
derived  relationship  between  dS^nan^c  j^essure  drop  and  dynamic  flow;  a 
series  of  calibrations  will  be  run  on  the  monopropellant  feed  S7stan  us- 
ing one  of  the  actual  rocket  motor  injectors  with  two  pressure  transducers 
and  a flouneter*  These  eallbratlons  will  consist  of  measuring  a pulsat* 

Ing  flow  with  one  of  the  dynamic  flowmeters;  at  the  sane  time  recording 
the  instantaneous  pressure  drop  across  an  injector*  According  to  theory; 
both  a phase  and  an  aaplitude  difference  between  flow  smd  pressure  drq;> 
are  introduced  as  a result  of  inertia  of  the  liquid*  Ihe  purpose  of  the 
calibration  tests  wrill  be  to  measure  these  differences;  and;  more  in^or- 
tant;  to  determine  tdiether  or  not  they  are  consistent  at  various  flow 
conditions*  These  tests  will  also  be  used  to  evaluate  the  effect  of  wave 
t raved,  time  in  the  feed  line,  by  ^wrying  the  distance  between  flowmeter 
and  injector  and  measuring  the  t:lme  lag  between  flowmeter  pulses  and  press- 
ure indicator  pulses*  If  the  time  leg  (idiich  is  actually  the  phase  diff- 
erence between  flow  and  pressure)  is  then  plotted  against  distance,  the 
intercept  at  zero  distance  will  be  the  true  phase  difference  between  flow 
and  pressure  at  the  injector  itself,  and  will  be  independent  of  the  dis- 
tance between  the  flowmeter  and  the  pressure  transducers* 
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IV»  APPMUS  (cmt*«a,) 

III  wwttwini  HI— wintiWK— ri 

0*  Insfcrujaentaticm  (Gont'd,) 

(i)  Senatog  Elements  (Gont*d.) 

In  addition  to  the  djnmtoic  floMnetei*a«  the  average  flow  rate  id.ll  be 
recorded  a Pottez*  floi«aeter$  vhS.oh  is  a pressure  balanced  turbine  «hoae 
speed  Is  proportional  to  the  flo»  rate*  Ihe  revolutions  of  the  turbine 
are  counted  and  Integrated  by  an  <^lectronic  circuit  (constructed  at  Brince- 
ton)  idilch  indicates  the  turbine  BPH>  and  hence  the  flow  rate>  on  a re» 
cording  potentloneter*  This  meter  will  be  used  in  conjunction  with  the 
dynamic  floMaeter  testa  desoribsd  above  in  order  to  evaluate  the  average 
flow  indications  of  the  various  dynamic  asters*  Calibration  of  the  Potter 
meter  can  be  easily  accos^liched  with  a stopwatch  and  weighing  tank$  since 
it  is  used  only  to  measure  steady  flow* 


(S)  Interaediate  Elements 

(a)  Hlgb-stabillty  0,0,  diffei’sntial  ai^pHfier 
A 0«C,  an^llfier  was  designed  by  the  Advance  Electronics  Corpor- 
d^lon,  FassaiCf  N.J,  in  cooperation  with  l^inceton  for  particular  appli> 
cation  to  the  type  of  measurements  required  for  the  subject  cmtraet* 

The  amplifier  is  distinguished  by  several  outstanding  features^  includ- 
ing excellent  8tabillty«  which  results  from  the  600  volt  "B"  voltage^ 
an  output  drift  of  only  about  *030  millivolts  per  hour«  and^  most  i»> 
portant,  a circuit  idiich  will  an^lify  the  difference  between  two  input 
signals*  Careful  attention  was  given  to  the  maintainance  of  sero  phase 
shift  in  an  oscillating  inpxxt  -signali;  Six  chann^s*  each  idth  variable 
gain  of  frMt  one  to  twenty,  were  built  in  three  rack-mounted  units  of 
two  channels  eaehr  Initial  testing  of  the  amplifiers  indicates  ex- 

I 

cellent  performance,  conformijig  to  all  design  specifications,  and  as 
soon  as  two  pressure  picktps  have  been  calibrated,  the  anpllflera  will 
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D«.  IjistjnMmtfttiaa  (Cont'd«) 

(2)  Intermediate  Elements  (Cont*d«) 

(a)  High»stability  E,C»  differential  amplifier 
be  obeoked  tmder  operational  conditions* 

(b)  Harrow  band-pass  filter 

one  of  the  major  problems  anconntered  in  setting  xi^  instrament«» 
tioa  for  tims-lsg  measureraents  in  rocket  motors  was  the  design  of  a 
high-fidelity  ultra-low  frequency  band-pass  filter  with  a band  width 
of  the  order  of  20  ops  and  sero  phase  shift  of  the  input  signal*  Sm» 
haustitre  searches  of  available  ecHnmerclal  equipment  revealed  no  single 
filter  or  combination  of  filters  which  could  meet  these  specifications* 
Sveral  bids  involving  special  designs  were  submitted  and  reviewedf  and 
it  was  finally  d<i!«ided  to  accept  that  of  3eva  Laboratorletf  Trenton« 
H.J*»  for  a four-ohamel  frequency-modulated  filter. 

Basically^  the  filter  consists  of  a frequency-selectiva  heterodyne 
anpllfier  continuously  adjustable  from  50  cycles  to  10,000  cycles  per 
sec(wid*  llie  entering  signal  is  modulated  by  the  output  from  a higb- 
frequanoy  oscillator  with  a drift  stability  better  tlian  3 cycles  per 
week*  The  resulting  «b©atw  signal  can  then  be  filtered  with  high 
selectivily  by  a multiplier  ciretiit,  obtaining  a froquenqy  band 
of  well  below  20  cps  outside  of  idiich  the  cutoff  is  extremely  sharp* 

The  frequency  signal  Is  then  demodulated,  and  the  filtered  output  passes 
on  to  the  recording  instrum^t* 

PreHmlnary  tests  by  the  manufacturer  indicate  very  small  phase- 
shift,  ^.ch,  if  it  exists  at  all,  is  nonstant,  and  hence  can  be  com- 
pensated for*  The  noise  output  of  the  instrument  was  claimed  to  b® 
negligible^  Delivery  of  the  four-channel  unit  is  ejspected  In  Jaaiaary, 
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I?.  APPARATUS  (CoKt'd.) 

T|[  ■>  III  ll«»lHWIIIIf  HW  — IITO 

D.  Instrmeirtatioa  (Cont'd,) 

(2)  Intermsdiate  Elesttsnts  (Oc»at*d*) 

(fe)  Ham>w  ba»d-paes  filtea? 

and  perrormnce  tests  itdll  be  conducted  at  that  tisne* 

(©)  Othecr  interraediato  equipment  consists  of  a tw-ciiannel  A,C. 
audio-^u^lifieTj  designed  and  constructed  at  Princeton^  and  t90  Integra** 
tor  circuits,  one  designed  by  Reaction  Motors,  Inc,,  the  other  by  the 
Naval  Air  Rocket  Test  Station  at  Lake  Denmark,  (both  constructed  at 
I¥lnc9toa)  for  use  with  the  Potter  flovneter,  A bridge  circuit  and  an- 
pUfler  for  the  hot-idre  anemometer  were  on  hand  at  Princeton,  and  have 
been  modified  for  the  particular  application  of  this  instrument  as  des- 
cribed previously*  Muoh  aiuciliary  equipnent,  consisting  of  power  sup- 
plies, calilrating  and  testing  instnsaents,  multiple  switches,  rr'»ays, 
potentiometers,  etc*,  have  been  purchased  or  constmcted  ly  Princeton 
to  service  this  and  other  contracts* 

(3)  Indicating  and  Recording  Elements 

AH  recorders  and  indicating  instruments  xised  in  the  measurements 
of  rocket  stability  parameters  are  commercially  available  types,  and 
were  described  briefly  in  the  first  quarterly  r^ortj  however,  a 
slightly  more  detailed  account  of  each  instrument  might  be  of  interest 
here* 

(a)  Dual-bean  oscillo8c<9e 

The  dual-beaiB  oscilloscope,  nanufa'tursd  by  Dumont y is  th®  same 
as  an  ordinary  laboratory  cathode-ray  oscilloscope' except  that  all 
faeilltias  necessary  to  observe  two  input  signals  are  included*  Thus 
this  instrument  is  an  excellent  means  of  copparlisg  two  t-rsnsient  or 
osolllatliig  ffisasuroraents*  The  oscilloscope,  unfortunately,  is  not 
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D.  Xnstxttiientation  (Cont^d.) 

O)  Indicating  and  Recording  El<^»e&t0  (Oont^d.) 
(a)  Dual-bean  oacllloscope 


adequate  for  quantitative  neasrurements  of  the  required  accuracy^,  partly 


due  to  the  high  ratio  of  trace  ;ddth  to  the  scale  of  coordinates^  and 


partly  due  to  its  high  aero  drift*  However*  its  extremely  high  fre- 


quency response  makes  it  an  excellent  means  of  studying  and  con^iaring 


pheacm«m  such  as  transient  inressures  in  the  injector  and  chaaber  dur« 


ing  ignition*  usitig  a Dumont  oscilloscope  camera  to  record  its  data* 


The  osoilloscope's  extreme  flexibility  also  makes  it  a vital  tool  in 


ths  qualitative  comparison  of  tso  mutually  dependent  quantities  (e«g** 
Injector  and  chamber  pressure*  injector  pressure  drop  and  flow*  dy- 


namic flow  aa  aeasurad  by  two  different  meters*  etc*)*  which  is  a 


necessity  in  order  to  determine  such  things  as  the  usefulness  of  the 


data*  the  selection  of  optimum  rocket  motor  test  conditions*  selection 
of  the  range  and  type  of  high*accura(y  recording  device*  etc* 

In  addition  to  the  above  direct  applications  to  rocket  tests*  this 
instrument  has  been  put  to  good  use  for  calibration  and  operational 
testing  of  various  other  pieces  of  equipment*  and  this  will  probably 


continue  to  be  one  of  its  major  fields  of  operation. 


(b)  Recording  potentiometers 


Six  of  these  general-purpose  Leeds  and  Northrup  instruments  (two 
of  Which  were  purchased  by  the  subject  contract)  sre  available  for 


recording  steady-state  fleasureaieats  such  as  average  chaittoer  pressure* 
average  propellant  flow  rata*  average  injector  pressure  drop*  etc. 

In  addition  to  these  direct  applications,  tho  recording  potentiometers 


have  been  invaluable  fcr  general  calibraticsa  instruraents  due  to  their 
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17,  APPARATUS  (Cont’d.) 

D*  In3trunentation  (Qont<d,) 

! 

(3)  Indicafciag  and  Recorddjag  Eleaents  (Cont'd,)  1 

(b)  Recording  potentiometers  i 

excellent  accuracy,  the  inaxiraum  overall  error  in  their  readings  being 

only  - 0,3^  of  full  scale, 

(c)  T«o>channel  magnetic  tape  recorder 

1 

The  tape  recorder  Hill  be  perhaps  the  most  important  single  unit 
in  the  electronic  instromentati^  system.  It  is  a highly  specialised 
model  designed  particularly  for  recording  rocket  combustion  data,  and 
has  several  features  Which  make  it  unique  in  this  applicaticsi. 

The  basic  recorder  unit,  amnufaetured  by  Aspex,  is  a high-fidelity 
instrument  with  less  than  \%  tape  flatter,  designed  for  industrial  re- 
search applications  utilising  electrical  iiput  signals.  The  basic 
unit  has  been  modified  so  that,  data  may  be  recorded  at  a tcpe  speed 
of  inches  per  second  and  played  back  at  either  60  or  3 inches  per 
second.  This  feature  provides  for  an  amplification  of  the  tine  scale 
of  an  oscillating  signal  by  a factor  of  twenty  to  one;  for  exaople, 
if  it  is  desired  to  record  a signal  oscillating  at  ^,000  cycles  per 
second,  use  of  the  low»speed  playback  actually  reduces  the  flreqi^noy 
to  250  cycles  per  second  without  distortion  of  the  signal  in  any  way, 
and  hence  the  high-frequency  input  signal  may  be  recorded  or  observed 
on  a low^firequency  instrument.  Furthermore,  if  it  is  desired  to  ob- 
serve the  wave  form  of  the  signal  on^^  say,  an  electromagnetic  osolHo- 
griph,  the  amplification  of  thae  scale  provided  by  the  tape  recordes^s 
lowospeed  playback  feature  actually  stret cites  out  tho  oscillograph 
trace  and  hence  peinnits  a clear  pictw©  of  tho  signal  to  b©  prasentod® 

This  feature  thus  multiplies  the  effective  tape  writiag  r,peeds  of  ■ 
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IV.  APPAEarUS  (Coat'd.) 

D®  lastrumeafcation  (Coat'd.) 

(3)  Xndieating  aad  Reoordiag  Elesneats  (Conti^d®) 

(o)  Two-chaaaol  magnetic  tape  reoordea? 

the  oociHogJPaph  a factcir  of  twenty  wltlMsut  lose  of  fidelity* 

One  of  th  j most  In^jortant  adymntagee  of  usiag  the  Anpex  machine  to 

record  data  is  that  .it  prorldos  a permanent  "signal  source"  ^Moh  can 

be  played  back  as  many  times  as  desired  into  different  types  of  display 

% 

instruments j|  with  different  filter  bandsj»  different  degrees  of  ai^}lifiea«> 
tion,  etc*  thus*  one  may  study  several  different  phenomena  of  a given 
test  carefully  «id  repeatedly*  selecting  optimum  instrument  ranges* 
rath^  tim  being  obliged  to  use  data  recorded  at  a single  set  of  in- 
strument settings  as  would  be  the  case  if  a non-reproducible  recorded 
trace  (such  as  the  oscillograph  tape)  were  used* 

The  Jatptex  has  arrived  recently*  but  unfortunately  one  of  the 
electrln  circuits  was  in  a damaged  condition  and  had  to  be  returned 
to  the  vendor  for  repairs*  Consequently*  only  brief  operaticnal  checks 
have  been  performed  to  date  on  this  instrument* 

(d)  Phasemeter 

!!he  phasemeter  is  an  instinment  idiich  measures  the  phase  difference 
between  two  input  signals*  This  phase  difference  may  be  observed  cn  a 
meter*  and  at  the  same  time  may  be  delivered  as  an  electrical  signal  to 
a recording  device  such  as  a recording  potontioneter  (for  a constant 
phase  difference)  or  an  eleetrccagnetie  oscillograph  (if  the  phase  diff- 
o?ence  varies  rapidly  with  time)  5 This  instrusent  thus  perforiis  olee° 
tronioallja  and  hence  aocuratGlj^**  m operation  Tshich  would  b©  both 
tteo-conoumiKg  and  probably  not  too  acoiu'ato  if  tho  phase  diffCTcaoc 
had  to  ba  msasured  sanually  freoj,  s©y*  ooeiliegraph  traeoe  ©,?  the 
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SEOTOWnratolON 

I?*  APPilRATgS  (Cont«d.) 

D*  Inst ruiaeniat ion  (Cont’d*)  * 

(3)  Indicating  and  Recording  Elem^ts  (Contsd®) 

(e)  Sonic  Analyser 

Delivery  of  the  sonic  analyser  from  Panoramic  has  been  delayed  tm-® 
til  January,  1953»  This  instrment  performs  electronically  a Fotiries* 
analysis  of  an  oscillating  electrical  sig^iol,  displaying  an  oacillo* 
scops  trace  (uhlch  can  be  recorded  by  a standard  oscilloscopt  camera) 
of  amplitude  versus  freqmncy*  The  Panoramic  analyser  nlH  be  used  as 
a olnuiiMn*  or  injector  pressure  monitor  for  most  of  the  e:^erim3ntal  runs, 
to  detect  any  oscillations  occurring  at  frequencies  other  than  the  im» 
posed  frequency  and  to  indicate  the  order  of  magnitude  of  the  ratio  be> 
tiiesn  the  imposed  pressure  oscillations  and  extraneous  combustion  tran» 
sients*  Thus,  this  Instrument  perforas  a valuable  function  in  that  it 
indicates  the  probability  of  obtaining  usefnl  data  from  a given  test 
prior  to  the  time-consuming  acourate  analysis*  Further,  it  provides  a 
record  of  all  frequencies  occurring  daring  any  test,  and  thus  permits 
selection  of  those  frequency  ranges  which  should  be  subjected  to  dd- 
tailed  investigation* 

(f)  Six-channel  recording  oscillograph 

The  terminal  display  instrument  for  all  dynamic  recorded  data  is  a 
six-channel  Hathaway  oscillograph  id^ilising  electromagnetic  light-beam 
galvanometers  idiich  trace  a signal  on  photosensitlvo  paper*  The  galvano- 
meters have  a flat  fre  ncy  respoase  up  to  1^00  cycles  per  oecond  (nat- 
ural frequency  about  3500  cps),  and  thus  provide  fear  accurate  aeasursmient 
of  oscillations  up  to  30,000  cps  idion  used  in  conjunction  with  the  Ampex 
tiqse  recorder.  The  usual  typo  of  synchronous  time  aarto  is  oispl<^ad^ 
and  the  chart  speed  may  bo  varied  in  steps  ;f>om  O9I  to  200  Inches  per 

second*  mm  the  high  chart  .'^ecd  is  iw3d  to  record  (iornbuatieii  dctta 
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X?c  APPARATgS  (Cont»d.) 

D,  Instrunentatien  (C<mt»d*) 

(3)  Indicating  and  Recordii^  E3.emants  (Cont’d®) 

(f)  Slx«channel  recording  oscillograph  (cont*d*) 
taken  on  the  tip»  recorder#  a stnglis  cjrclo  of  a signal  oscillating  at 
ijjOOO  cycles  per  second  will  occtspy  a fall  inch  of  oscillograph  tape# 
and  thus  the  pressuro  wavs  forms  occurring  during  the  high«»fr(i5quenc3r 
instability  may  be  carefully  analysed* 

Basic  instrumentation  research  has  progressed  to  the  point  at 
idiich  it  is  now  possible  to  set  up  a fomaliaed  test  program*  Of 
course#  developments  during  the  course  of  testisig  may  require  esr»> 
tain  changes  to  be  made#  but  the  Imtsic  program  will  be  undertaken  as 
followst 

!•  Calibration  of  individual  sensing  and  recording  elements#  de- 
fining linearity,  range,  accuracy#  repeatability,  frequency  response# 
and  applicability* 

2*  Comparative  calibration  of  the  dynamic  flowmeters  as  discussed 
earlier  in  this  section  (see  Fig*  6)* 

3«  Bstermination  of  dynamic  flow  as  a function  of  dynamic  pressure 
drop  across  an  injector  orifice  (see  Fig*  7)* 

(&)  Detenaination  of  wave  travel  time  in  the  feed  system  and  its 
effects  on  dynamic  flow  and  pressure  measurements* 

It*  Calibration  of  all  instrumentation  as  arranged  for  rocket  motor 
tests* 

5*  Measurement  of  combustion  time  lags  in  the  Konopropellaat  racket 
motor  system  (see  Pig®  8)* 

As  has  been  mentioned  in  previous  reports#  b.  very  largo  percentage 
of  the  effort  in  this  invest igatS.on  has  been  ovpended  on,  instrmmt-at.ion# 
since  the  s^sitiTO  Keaaurcmont  of  conburt,ion  p;3ra>'*-.yterc  i‘oi*5:?.3  the  tesio 
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IV.  APPARjgUS  (Cont«d.) 

D.  Instrumentation  (Cont’d.) 

(3)  Indicating  and  Heoording  Elements  (Coat’d.) 

(f)  Six»channel  recording  oscillograph  (Cont’d.) 
for  any  experimental  work  concerning  rocket  motor  stability.  It  is 
necessary  that  this  tj^pe  of  research  continue  without  intermpti«m  for 
at  least  the  next  six  months  in  order  to  fulfill  the  requirements  for 
accurate  rocket  combustion  data«  and  consequently  instrumentation 
studies  will  be  pursued  both  in  connection  with  and  in  addition  to 
monopropellant  rocket  motor  operation  during  the  forthcoming  report 
period* 
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?<,  Um'Mm'i-idli  mD  mi‘k 

Xmtvment  GiUlbratiMB 

(3.)  ri€ae?“al 

Tha  mjes"  eisphasis  during  t^ia  pariod  mti  plass’.d  m salibr&tiCK© 
aad  opaffatrionel  ieafts  o£  the  jH'ossis’*®  t.raQBdKss«»j>  whisfe 

«dll  ba  doscribed  bal0j«  at  seme  (|)«ratlc»tel  @h@3ks  «£  all 

intermodlate  and  re‘2c*s‘ding  mr&  made,  tka  exception 

of  the  four-ohannel  filter  and  the  Panorando  sonic  analyser,  i^hieh 
hare  not  yet  been  dellv^ed«  Preliadnery  testing  of  the  Anipex  tape 
re$<^'der  revealed  that  one  of  the  oirouite  had  been  delivered  in  a 
daiieaged  oandltion,  and  this  section  of  the  instruiaent  aas  returned 
to  the  ■aimfaoturer  for  r6palrs« 

Hone  of  the  (^aiale  flo^aneteo^s  have  been  received  as  of  the  date 
of  this  reports  Strength  tests  on  the  hot  aire  sensing  elsaent  have 
been  eondacted,  and  a method  for  laount  ig  the  hot  %dre  so  that  it  will 
not  break  at  high  flow  rates  under  pressure  has  been  devised*  Modlfi* 
0£tlon  of  available  electronic  aTnriliary  equipment  has  been  completed, 
md  the  hot  wire  is  ready  for  operational  testing,  awaiting  delivery 


of  ""he  flow  Bulsix^  unit* 

B?eliHinary  testa  of  the  Petto®*  flow  Kater  were  mde  using  an 
alsctrmle  integrsitor  d®gign©d  by  Eaactlca  Motors,  Ince,  and  built 
Pris^setontj  Thesa  teste  had  to  b©  run  at  low  flow  rates  under 
^'mtar’  pressure,  sine©  late  de3.iverlo8  of  high-pressure  equipisent  de- 
.layed  oC'^lotisn  of  th©  eosponeRst  .tost  panelj  however,  the  tests  ia«* 
dicatf4  tliat  performance  could  be  liiq>ro'^8d  by  uso  of  a new  type  of 


lKte,^/alGrs)  dosignt^d  racently  at  t,h-3  Rocket  Test  Ststios*, 

'):his  ftciij  mil-  K;vs  just  bofn  coiuplcitcd,  and  th©  Fottor  no  ter  will  be 
o;  OAh-.-i/i.vd  I'-rmi  con^iotion  of  efefe^r  touto  now  utilising  tho 
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V»  IKFORMATION  AND  DATA  (CQat»d.) 

w 

A*  Xn3trttmer»t  Calibrations  (Cont«d«) 
(1)  General  (Cont’d.) 
test  stand. 


(2)  pressure^senaing  elements 

feats  on  the  Princetai-MIf  pressure  transducer  were  dlidded  i.nU>  two 
categoriest  (a)  Operational  tests  under  actual  rocket  ssotor  operating  con« 
ditioaas  (performed  by  NACA  staff  Members  at  the  NACA  Rocket  Laborati^  in 
Cleveland) y and  (b)  Static  cala.brations  of  the  differential  unit  at  prince^ 
ton* 

(a)  Operational  Tests 

Before  Installation  of  the  dumsy  picktq:>  in  a rocket  chamber,  a short 
series  of  tests  was  made  to  determine  the  water  flow  rate  through  the 
pickup*  It  was  found  that  the  flow  rate  obeyed  the  usual  orifice  squatlont 

M 

w = s-VAP 

where  K = *00558 

W = witer  flow  rate  (Ib/eec) 
and  AP=  water  pressure  drop  across  pickup  (psi) 


For  a water  flow  rate  of  0*035  pouads  per  second^  corx'eoponding  to  a 
pressure  drop  of  36  psi  across  the  pickup,  the  allowable  heat  transfer 
across  the  diaphragm  for  a coolant  temperature  rise  of  100®F  would  b© 

23*3  Btu/s0c/sq»  in*  This  is  based  on  the  diaphragm  area  of  0*15  in® 
If  one  iacludos  the  enbijS’e  pickup  face  (0*37  sq.  in*)  the  acceptable  hsat 
transfer  would  still  bo  Btu/sq,  ins/sec*,  which  ma  believed  to  bo 
adequate  cooling  for  rocket  cSaaifl>a¥'  usso 

The  pickup  units  vcro  fhoa  tectad  in  a rocket  engine*  The  rocket 
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INfOBMATION  AMD  DATA  (Coat'd.) 

A.  lastrumeat  CallbrAticms  (Coat'd.) 

(Si)  Pressuro-nensing  elenoati  (Coat'd.) 

(a)  Operational  Teats  (Oent'd.) 

osed  'MS  a wat<H»«cooled  ohamber  of  l^«ineh  diaaeter  SU  inolMS  long»  equipped 
with  a water*eooled  noszle  whloh  was  rated  at  ^ pounds  thrust  at  300  psia 
ohanbar  pressure.  Two  t;^es  of  injeetors  were  used.,  one  a showsrhead  type 
ai3d  one  an  iiaplnglng  Jet  type,  ^e  showerhead  injector  was  equipped  with 
a pressure  tap  to  afford  a measure  the  chamber  pressure*  whereas  the 
impinging  Jet  injector  was  net.  A boss  drilled  and  tapped  for  m 18  milli- 
meter spark  plug  thread  was  Installed  in  the  eyliadrioal  part  of  the  cham- 
ber 5«3A  inches  from  the  nozsle  throat  (2-lA  inches  from  start  nossle). 
The  dumsy  piolnqp  was  installed  in  the  boas  so  that  the  diaphragm  wms  rs« 
ssssed  about  1/32*  from  the  chamber  wall  surface. 

The  propellant  combinatiooa  used  vrere  liquid  oxygen  - ammonia  and 
liquid  caygeo  - JP3.  Previous  ejqperioace  with  the  rocket  had  indicated 
heat-transfer  rates  in  the  chamber  of  about  O.U  Btu/sq.  in./see  under 
normal  (non-screaming)  operations  and  2.1  Btu/sq«  in./sec  under  screaming 
operations  with  ammonia  as  fuel.  The  performance  based  on  ratio  of  ob- 
served to  theoretical  specific  impvilse  was  about  1$  percent  under  non- 
screaming  conditions  and  about  97  percent  under  screamirtg  conditions  with 
no  correction  for  heat  Injeotion.  The  accuracy  of  the  performance  figures 
is  about  ^ percent. 

liflth  JP3  as  fuel  the  heat-transfar  rates  wejro  about  50  percent  higher 
than  with  anmonia,  and  the  performance  figuros  were  substantially  unaffected. 

A series  of  ten  runs  was  made.  The  dummy  pickup  was  used  in  the  first 

seven  and  in  the  ninth  and  tenth  r^ms«  The  first  six  of  these  used  aamonia 
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a'FOPJ»M’ICM  jySD  MA  (Cont«d,) 

4»  Iu.»fer«&wfc  Cslteatioas  (Cosat*d*) 

(S)  Br©38\ir®«senslag  elements  (Confc»de) 

(a)  Operational  Tests  (C«%t*d») 

as  the  fnelg)  the  others  JP3*  The  dmssy  siu^ned  no  evldenoe  deterioration 

with  either  the  asmtonla  or  JF3  fuels  prior  to  the  last  rtm.  ¥lth  J!P3  fu^ 

( 

a eonslderable  fihe  of  carbon  was  formed  on  the  diaphragm  (see  Fig«  9)» 

This  could  b?  removed  easily  (Fig*  10)  and  apparently  did  not  harm  the 
piektp*  The  measured  heat  transfer  to  the  picki^,  as  determined  frcei 
the  teaye'fe.txsFe  tho  ester  teto  jsad  oust  of  the  pickwg  wa,«  Btn/«ec/ 
aq«  in*  during  a noroial  run  with  the  amaonia  fuel*  A non»soreaBing  run 
with  J?3  fuel)  indicated  a heat-transfer  rate  of  0*75  Btu/sq*  ine/see  * 

The  "live"  picki^  was  installed  on  the  eighth  run*  This  run  broke 
Into  screaming  after  3/h  seconds  of  nomal  operation  and  the  pickup  me 
destroyed  after  another  1/2  second  (see  Fig*  U)*  Visual  observation  of 
the  oscilloscope  traces  indicated  very  satisfactory  pressure  response  up 
to  the  nKsaent  of  failure*  A further  test  under  similar  cperatixig  con- 
dltionS;»  again  using  a dunay  pickup  (run  nine)#  resulted  in  no  scre^uning 
and  a heat  transfer  rate  of  0,75  in. /sec*  A seccad  ran  (run  10) 

again  broke  into  screaming.  This  destroyed  the  dumy  pickup#  ^ch  ap- 
parently failed  immediately  after  screaMng  started*  The  run  my  have 
been  an  unusually  severe  screaming  run#  since  the  chamber  wall  burned 
threugh  at  about  the  midpoint  of  the  chaitter#  and  this  chamber  had  with- 
stood many  leagtl^  screaming  runs  without  previous  dasaige®  Th®  pickup 
run  conditions  and  obaoi^vations  are  tabulated  iii  Table  1* 

It  ia  quit®  difficult  to  estimate  accurately  th©  he?,t  transfer  t© 
the  plcb'po  dui'iag  screoFlng  nHao®  Jueb  bofora  failuT©  th©  t©c!> 

ps’-aturn  out  of  tho  live  piekuu  (*nm  8)  acmad  to  lovol  off  at  a teR® 
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FIGURV.  9 


TRa;r;iUC  'R  AFTI®  NURMfJ'.  aoCK*^!  MOTOR 

'ii'm  r..iQGiT)  c:aGT',N-jF3 


can  be  mad 


section  twf  this  vm  ms  3«2 


' (b)  io  CsUb^stion 

t^sliaiinsrjr  cheeks  of  the  thr««  different  isl  trsnsducws  delivered 
to  Priaoston  revealed  the  exist  moe  of  leakage  from  the  reference  pres<» 
sure  ohsKber  above  the  \3pper  diaphragm  (see  Plg»  U)  in  all  three  units. 
Teo  of  the  picki^s  msis  retumod  to  MIT  for  correction  of  this  oondition, 
end  the  third  one,  idiich  had  <»il^  a slight  leak,  vras  temporarily  r®» 
tained  at  Princeton  for  calibration» 

1!he  single  pickup  ms  pressurised  with  air,  separate  c^tr<^$  par- 


fh»  pressures  wsr®  ©bsorved  on  12*^  Heise  Bowdon  tube  gages  with 
a combined  reading  and  instrument  error  of  - 1 psi  for  the  range  used® 

A Eubicon  hand-balanced  potentiometer  with  an  estimated  maxtem  error 

‘ RCTHICFED 


TABLE  I,  PRlNCETOfJ -M:T  PRESSURE  PICKUP'"'/-  OPERATIONAL  DATA^  NOV.  4-10,1952 

OXIDANT  = ^IQl'iD  OXYGEN,  FUEL  ^ AS  NOTED-  AREA  OF  PICKUP  DIAPHRAGM  ^ TAKEN  AS  0,15  SQ.  IN. 


RUN 

FUEL 

USED 

CHAMBER 

PRESSURE 

{PSlAy 

THRUS'r 

(LB) 

TOTAL 
PROPELLANT 
WT  FLOW 

ilb/sec.) 

0 

F 

R.ATiO 

ACTUAL 

SPECIFIC 

IMPULSE 

(SEC.: 

EFFICIENCY 

(%) 

CHAMBER 

TEMPERATURE 

(Op) 

PICKUP 

COOLANT 

TEMPERATURE 

T ...  (*  F ^ 

IN  ^ 

PICKUP 

COOLANT 

TEMPERATURE 

WATER 
COOLANT 
FLOW 
THRU 
PICKUP 
(LB. /SEC. 

HEAT 

TRANSFER 

ACROSS 

DIAPHRAGM 

(BTU/1N2 

SEC) 

HEAT 

TRANSFER 

ACROSS 

CHAMBER 

IBTU/ 

IN^SEC.) 

PICKUP  I 
USED 

U — ' 

1 

NH  2 

270^ 

488 

3.28 

2.01 

: 49 

b?.:^ 

2030* 

0,034 

— 

• 

DUMMY 

i 

2 

NH  ^ 

275^ 

49c 

3.35 

1.07 

148 

bC.4 

1650 

0.031 

.3 

NH3 

270^ 

510 

3.33 

1.48 

153 

62.5 

1910 

77.3 

79. 4 

0.03  4 

0.048 

• 

1) 

4 

NH3 

270^ 

264 

2.2  2 

2.47 

107 

4 4.6 

890 

0.03C 

)| 

5 

270^ 

530 

2.79 

1.97 

1 90 

84.4 

3420 

y 

: 

'■  3' ' '.'L'r  ■ 

6 

7 

NH3 

27C^ 

3C5 

L.08  . 

1.21 

147 
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3.20 

2 20 
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56.0 

0.037 
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1.  Chamber  temperature  was  obtained  by  multiplying  the  theoretical  comou^tion  charnber  temp^^rature  for  a particular  fuel  - ox.danf  comojsr.on  by  the  square 

of  the  effinien'y.  .tx  u j 

2 Snower-heQQ  type  injector  witr.ouT  pressure  tap  wes  .sed  witr.  rnese  runs,  affording  comcastion  cr.amber  pressure  records. 

3.  Impinging  - jet  injector  without  pressure  tap  was  used  witn  these  runs. 

4.  Screaming  runs  resulting  in  pickup  failure, 

5 Area  of  the  pickup  diaphragm  was  taken  as  0.15  sq.  in;  This  s the  area  of  the  diaphragm  Tself,  not  the  proje:ted  area  of  the  diaphragm  end  of 
the  pickup.  Diaphragm  is  7/16  in,  in  diameter,  and  basic  diameter  of  the  pickup  body  at  diaphragm  end  is  9/16  inches. 

6.  Length  of  chamber  is  24  inches chamber  and  nozzle  were  water  cooled. 

7 Pickup  wus  locoted  5-3/4  inches  from  the  nozzle  throat,  or  2-  1/4  inches  from  junction  between  the  nozzle  and  chamber, 
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ssctiRSFOTSffitfioi 

?e  IMFOEMAKOH  MP  DATA  (Goaisd.) 

A®  lasfcruinsnt  Calibrationa  (Ccmt«d«) 

(2)  Ps’essur©“®easiag  ©lemeat®  (Coat'd®) 

(b)  Static  Calibration  (Ccait»d») 

©f  «»  Oal  sdllivolts  ms  used  t©  measure  the  transducer  output®  Precis* 
ioa  resistors  were  used  for  the  oxtemal  bridge  drovtitj,  md  the  pickup 
bri«%®  is^ut  current  j moult  wed  a sdlliassmetor,  ms  matotained  with* 

la  * 0*1  Billias^er©  (0s3^)  al3,  runs®  A diagram  of  th®  in®trms«tS“ 
tim  is  included  la  Fig«  12® 

Th@  reference  pressure  was  fixed  at  certain  values*  and  picks^  out* 
put  wa®  reccsrded  as  a function  of  simulated  chamber  presswe®  A t^ical 
plot  of  the  data  at  one  value  of  the  reference  pressure  is  givea  te  Fig® 

13®  It  ms  found  that  the  slope  of  the  linear  portion  of  these  plots 
is  sensitive  to  input  current,  rather  than  to  ir^ut  voltage,  and  that  a 
minteum  warmiig)  time  of  20  minutes  is  required  for  the  transducw  electri- 
cal sjstera  in  order  to  prevent  zero  drift  frran  occurring® 

Coiasisteaogr  of  operation  between  two  pickups  could  not  bo  checked, 
since  only  one  unit  was  available  fcr  test,  but  repeatability  of  readings 
was  within  provided  input  current  ms  maintained  constant®  A discus- 
sicsa  of  the  data  spears  in  the  following  section®  . 
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FIGURE  12 

TRUMENTATION  FOR  STATIC  CALIBRATION 
OF  PRESSURE  TRANSDUCER 


DISOOSSION 

A«  GaXl'bx>atiai» 

t 

!•  f^essure  seiuiisg 

(ft)  (^Wfttional  tftfttft 

l%ift  fftll'ore  'of  the  S^inc^oe^HXf  tmtsducsr  mdw  sos°@a»isig  ' 

eoaoSitlOKift  in  & rocket  netor  d^nstrates  that  this  in@tnm@nt  is  not  m 
gr«t  sfttislftertory  for  the  Invest Igatlon  of  instabillt;^  phonmena  in  gmeralf 
honevftr^  the  data  of  fable  1 Indicate  its  ability  to  withstand  ai^  teiHpem« 
tTires  and  pressures  likely  to  be  moom^ered  in  the  aonoprc^tllaiit  rocket 
motor  tests  nsiog  ethylene  oxide»  llhnsj^  the  nonopropellimt  program  sill  , 
utilise  the  transducer  in  its  present  t&emg  but  additional  developn^t 
«LU  be  oarrled  on  in  the  meantiae  in  order  to  increase  the  tsa^rature 
range  of  the  plekiq»  for  application  in  the  proposed  bipropeHant  eoa- 
bustion  inTestigftticn* 

« 

Tha  probleoi  of  leakage  fron  the  baek-pressure  chamber  of  the  pickup 
is  beiag  attacked  by  both  PrlnoetoK  and  MX7*  fhe  difficulty  lies  in 
obtcdnlng  an  adequate  around  the  electrical  terminals  at  the  top  of 
the  transducer  (see  Fig.  I;),  In  the  present  instrument,  sealing  was  at» 
tespted  by  the  use  of  red  glyptal  in  the  very  small  clearances  betufesn 
the  piokiq^  body*  insulators,  and  aetal  terminal  rods,  but  this  method  ap» 
pears  to  be  satisfaetoxy’  only  at  very  low  back  pressures.  Several  possible 
remedies-  have  been  considered,  and  the  cooperative  effort  of  Priacetog?  and ' 
Mlf  is  egqpeeted  to  resolve  the  problem  ulthout  too  much  difficulty^ 

Several  important  characteristics  of  the  transducer  are  revealed  by. 
¥5^i)  13  and  similar  plots^  Hynteresic  in  the  linear  range  is  aero  (to 
about  accuracy).  Although  bysteresis  line  bem  observed  in  the  non<» 
linear  portions  of  some  eali.brations,  the  ntaximw  spread  of  the  linear 
regi^as  of  the  ascmding  and  descending  data  cisrves  always  remains  lidth- 
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71.  PIEOli^SIOH  (G<mt»d.) 

; Ap  Inatyupent  Caliteations  (Cont*d*)  ' 

1*  Fr assure  sensing  elements  (Cont'd*) 

(i)  Operational  testa  (Ccnt*d*) 

in  e;q3erimental  error*  The  linear  regim  remains  nearlj  eonstant  re« 
gardXesa  of  the  nagnitude  of  the  back  pressure^  covering  the  range  froH 
about  $0  psi  to  about  300  psi  pressure  difference  across  the  diaphragms* 
This  portion  of  the  picku$»  characteristic  is  limited  on  the  lowmr  side 
uhen  the  dlaphragn  moves  cwpletelj  awaj-  from  the  strain  tuhe«  and  on 
the  tqjpar  side  bf  idiat  appears  to  be  the  action  of  the  orerload  saf«t;f 
atop  (see  Fig*  1|.)*  The  ?50  pel  linear  range  of  the  pai’tioular  pickup 
used  for  theee  tests  is  probably  insufficient  f<»r  Its  projected  appU«> 
eatloo»  but  indications  are  that  the  range  can  be  increased  kqr  adjust- 
Bent  of  the  safet7  stop  iMchanisii*  This  belief  is  based  on  the  fact 
that  one  of  the  rejected  transducers  was  tested  at  eero  back  pressure 
and  exhibited  a linear  range  of  about  bOO  psi^  and  also  that  the  V9xy 
snail  novement  of  the  strain  tube  under  load  requires  aliaost  Ineredibly 
fine  tolerances  in  the  safety  atop  ecnpcncentsy  making  likely  the  pos- 
sibility of  too  early  applicaticn  of  the  stop* 

The  construction  of  the  pickup  also  requires  that  In  order  to 
measure  pressure  variations  about  a given  mean  value«  it  is  necessary 
to  set  the  reference  pressure  not  at  the  asaan  value^  but  at  sooe  pres- 
sure well  below  it*  This  is  a resTilt  of  the  fact  that  the  strain  tube 
la  not  preloadedi  i,e*,  a slightly  negative  pressure  difference  across 
the  diaphragm  a them  out  of  contact  with  the  strain  tube*  Hence^ 

in  order  to  obtain  readings  ^ddch  oscillate  about  the  center  of  the 
pickup's  linear  range,  the  reference  pressure  wotild  be  set,  for  en^le, 

at  about  175  psi  below  the  mean  chamb(^  pressure  for  the  particulsr 
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’ QI5C0SSI0N  (Cont*d.) 

[ A*  Instnunent  CaUb2*atioiiL»  (Cont*d») 

X«  Br®s8ure  sensing  elements  (Co®fc«d.) 
j (a)  Operational  tests  (Cont*d») 

V 

I picktq;)  of  Fig.  13*  This  behaTior  introduces  no  difficulties  whatsoewr# 

but  ear«  must  be  taken  that  the  linear  rar  ^e  of  each  piokig}  used  in  the 
investigation  is  well  defined  at  all  required  iralues  of  back  pressure 
by  adequate  oallbratlon  tests, 

> As  was  mentioned  ear Her j reproducibility  of  data  was  accurate  to 

I within  1^  provided  the  input  current  was  maintained  constant*  Sero  drift 

I was  negligible  with  proper  warmup  time,  but  the  effect  of  tNg>erature 

t 

J 

j changes  has  not  yet  been  evaluated*  The  frequency  response  of  the  in* 

, < strument  must  be  measwed  under  operating  conditions*,  and  the  ef fisct  of 

I 

I varying  cooling  water  pressure  on  the  transducer  output,  if  it  exists, 

I must  be  detendned* 

i y 

Summarizing the  differential  transducer  in  its  present  form  appears 
' to  be  satisfactory  for  the  early  phases  of  aonopropellant  testing,  with 

' the  foHowiiig  provlsionst 

(a)  The  reference  pressure  chsAfeer  leakage  must  be  eliminated* 

I 

; (b)  The  linear  range  must  be  extended  by  about  $0%  (probably  by 

I adjustment  of  the  safety-stop  device), 

J . ' . 

I (c)  The  effect  of  temperatiare  variation  on  zero  drift  must  be 

i 

; evaluated,  ' 

I Consultation  with  MIT  on  the  first  two  items  is  expected  to  resolve 

i 

’ ^ them  shortly,  and  the  third  will  be  checked  during  the  next  report  period 

J * The  manufactui'er  quotes  a natural  frequency  of  about  26,000  cycles 

pea?  second  with  no  coolant  flow  between  the  diaphragms* 
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LIST  OF  SYMBOLS 


1*  dimensionless  instantaneous  value  of  the  time 

lag  reduced  by  the  characteristic  time 

^ f X 


X 

2.  'tJi  = dimensionless  steady  state  value  and  instantaneous 
' ’ value  of  the  part  of  the  total  time  lag  which  is 

insensitive  to  the  pressure  oscillaticn 

: 3.  ‘T.'T  3s  dimensionless  steady  state  value  and  instantaneous 

value  of  the  other  part  of  the  total  time  lag 
? which  is  sensitive  to  the  pressure  oscillationo 

a 4.  'n.  s exponent  in  the  pressure  dependence  of  the  time 

I lag  = pressure  index  of  interaction 

I . 5.  - rate  of  injection  of  the  propellant  per  unit 

i cross  sectional  area  of  the  combustion  chamber 

! ft 

i 6.  s rata  of  the  generation  of  the  hot  gas  from  the 

I combustion  of  the  propellant  per  unit  cross 

I sectional  area  of  the  combustion  chamber 

•I 

7.  U = combusion  chamber  length  from  injector  end  to  the 
entrance  of  the  de'Laval  nozzle  - reference  length 
scale  for  non-dimensionalization 

j %*  s distance  from  injector  end  along  the  combustion 

chamber 

j 9.  dimensional  time 

i 10,  U*  dimensional  mean  flow  speed  of  the  gas  along  the 

I . combustion  chamber  axis 

ii 

11,  f*>  siod  - dimensional  instantaneous  values  of 

* pressure,  density,  temperature  and 

speed  of  sound  in  the  burned  hot  gas 

12,  ~jf  and  c*  = values  of  f*,  t***  c5' 

I **  * ® at  injector  end 

:!  ‘ s reference  quantities  for  non- 

dimensionalization 

13,  *^1  ^ 1*  = steady  state  values  of  T* 


« ■ 


t= 

X s X*/  Lt 
u,  U = 

M.  M 

|t» , ^ , X and  C 
V , f > X and  c 


- characteristic  time  s time  required  for 
a sound  wave  to  travel  the  entire  length 
of  the  combustion  chamber  filled  with 
stagnant  burned  gas 

= dimensionless  time 

r dimensionless  length 

= dimensionless  velocity  of  the  gas  in  un- 
steady arid  steady  state  operation 

Mach  niimber  of  the  gas  flow  in  unsteady 
and  steady  state  operation 

=:  dimensionless  instantaneous  values  of 
pressure,  density,  temperature  and  speed 
of  sound 


and  c = dimensionless  steady  state  values  of 
^ ^ ^ , and  c _ 

f',  t’  and  c’  = dimensionless  instantaneous  perturbations 
' ' ' over  their  respective  steady  state  values 

c<  = V 1 tO  « root  of  the  characteristic  equation  with 

the  dimensionless  time  as  the  independent 
variable 

X = dimensionless  amplification  coefficient 

lO  s dimensionless  angular  frequency 

ri  = absolute  value  of  the  angular  frequency  = «>  / ® 

If  = adiabatic  index  of  the  combustion  gas 

s M,'  <.-dt) 

^ ” dimensionless  distance  of  the  concentrated  com- 

bustion front  from  the  injector  end  expressed 
as  a fraction  of  the  characteristic  length  L 

^ = reduced  angular  frequency  of  the  oscillation 

= angular  frequency  divided  by  the  velocity  gradient 

in  deLaval  nozzle 


= reduced  velocity  parameter 


% 


>1 


■» 


s the  ratio  of  fractional  variation  of 
velocity  to  fractional  variation  of 
density  at  combustion  chamber  exit 
or  entrance  to  deLaval  noszle 

34.  R.S  r real  and  imaginary  parts  of  Ifp.u) 

35 • K,  W = integers  characterizing  the  modes  of  the 

oscillation 

36.  subscript  X or  t means  partial  differentiation  with  respect 
to  X or  t 

37.  ^ = substantial  derivatives  along  the  path  of  a 

^ propellant  element 

subscript  1 or  2 denotes  the  quantities  evaluated  in  the 
flow  field  1 or  2 

39.  Of  ) = the  order  of  magnitude  of  the  quantity  in  the 
bracket 


I.  INTRODUCTION 


«-  X » 


Rough  combustion  as  a result  of  large  pressure  oscilla- 
tions i-n  the  combustion  chamber  of  a liquid  propellant  rocket 
motor  has  been  observed  under  different  circumstances  in  two 
distinct  ranges  of  frequencies,  the  low  frequency  range  of 
less  than  one  hundred  cycles  per  second  and  the  high  frequency 
range  of  several  hundreds  or  several  thousands  cycles  per 
second.  Such  rough  combustion  not  only  gives  fluctuating  per- 
formance but  also  shortens  the  life  of  the  rocket  motor-  An 
understanding  of  the  basic  mechanism  of  producing  unstable 
pressure  oscillations  that  lead  to  rough  combustion  is  there- 
fore of  great  practical  importance. 

It  has  long  been  recognized  that  for  unstable  operations 
the  oscillation  of  the  chamber  pressure  and  the  oscillation  of 
the  rate  of  hot  gas  generation  produced  by  the  pressure 
oscillation  must  be  properly  out  of  time  phase  so  that  an  in- 
crease of  the  rate  of  hot  gas  generation  occurs  at  an  over- 
pressure period  and  further  increases  the  over-pressure  in  the 
combustion  chamber.  This  time  phase  difference  is  originated 
from  the  fact  that  the  propellant  element  does  not  burn  im- 

4 

mediately  after  being  injected  into  the  combustion  chamber, 
but  burns  after  a certain  time  interval  called  the  "time 
lag"  during  which  the  fuel  and  the  oxidizer  particles  mix 
properly  and  absorb  the  necessary  amount  of  activation  energy. 
In  references  1,  2 and  3,  the  low  frequency  oscillation  has 
been  analyzed  based  on  the  assumption  that  the  time  lag  is 
constant  and  independent  of  the  oscillations  of  the  gas  system 
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f 


( 
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in  the  combustion  chamber.  In  these  analyses,  a pressure  sen- 
sitive feeding  system  which  provides  a varying  injecticn  rate 
under  the  pressure  oscillation  in  the  combustion  chamber  is 
assumed  to  be  the  self-exciting  mechanism  which  creates  the 
variation  of  the  rate  of  hot  gas  generation.  The  senior 
author  of  the  present  paper  has  pointed  out  in  reference  4 that 
the  pressure  sensitivity  of  the  feeding  system  is  not  a neces- 
sary self-exciting  mechanism  for  producing  unstable  pressure 
oscillations.  By  assuming  a varying  time  lag  X which  depends 
on  the  chamber  pressure  ^ in  the  manner 


i dt^  = constant 

h-r 


it  is  shown  that  unstable  oscillations  of  both  the  low  fre- 
quency range  and  the  high  frequency  range  can  be  produced  even 
if  the  injection  rate  is  constant.  This  kind  of  combustion 
instability  cannot  be  eliminated  by  properly  designing  the 
feeding  system  but  is  intrinsic  in  the  nature  of  the  combustion 
processesy 

The  case  of  low  frequency  intrinsic  combustion  instability 
is  extensively  studied  by  the  senior  author  in  reference  4 while 

the  high  frequency  case  is  only  briefly  discussed  using  a 

« 

simplified  model  of  a single  concentrated  combustion  front  near 
the  injector  end  and  for  the  particular  value  of  n = 
where  is  the  adiabatic  index  of  the  burned  gas.  As  n 
represents  the  extent  of  interaction  between  the  pressure 
oscillation  and  the  combustion  processes,  n must  be  an 


important  parameter  as  a stability  criterion.  In  practical 
cases,  a large  part  of  the  corabnstion  often  takes  place  in  a 
narrow  region  somewhere  in  the  combustion  chamber  but  usually 
not  too  close  to  the  injector  end.  Therefore  the  analysis 
in  reference  4 is  extended  in  the  present  paper  to  the  case 
with  arbitrary  location  of  the  concentrated  combustion  front 
and  for  arbitrary  values  of  n , The  effect  of  distributing 
the  combustion  along  the  combustion  chamber  axis  presents  a 
more  difficult  problem  and  is  analysed  in  reference  5®  The 
problem  of  shifting  part  of  the  concentrated  combustion  at 
the  injector  end  to  arbitrary  axial  location  is  briefly 
studied  in  the  present  paper  with  a view  to  some  indication  of 
the  effect  of  distributing  the  combustion  axially. 

The  unsteady  supercritical  flow  in  a deLaval  nozzle  with 
linear  steady  state  velocity  profile  in  the  subsonic  portion 
has  been  studied  in  references  6 and  7 from  which  the  boundai’y 
condition  for  high  frequency  oscillations  can  be  deduced. 

This  boundaiy  condition  will  be  used  in  several  specific 
examples.  In  the  limiting  case  of  a very  short  nozzle,  this 
boundairy  condition  is  equivalent  to  the  boundary  condition  of 
constant  flow  Mach  number  at  the  entrance  of  the  nozzle  which 
is  used  in  reference  4«  This  short  nozzle  boundary  condition 
is  simple  and  admits  analytical  solution  of  the  characteristic 
value  problem.  The  results  with  these  two  boundary  conditions 
will  be  compared. 

II.  FORMULATION  OF  THE  PROBLEM 
2ol  Assumptions  and  Simplified  Models  of  Gas  Flow  System 

The  combustion  chamber  of  a liquid  propellant  rocket 


motor  is  often  a straight  duct  of  constant  ci'^oss  sectional 
area  and  is  filled  with  the  hot  burned  gas.  The  propellant 
elements  injected  into  the  combustion  chamber  ar^e  mostly  in 
the  form  of  atomized  liquid  droplets  suspended  in  and  carried 
along  by  the  burned  gas  stream  without  occupying  app"*eciable 
volume  of  the  combustion  chamber.  The  hot  gas  generated  from 
combustion  is  recirculated  actively  to  the  region  near  the 
injector  end  where  the  hot  gas  supplies  the  activation  energy 
to  the  unburned  propellant  elements.  The  recirculation  and 
the  mean  flow  patterns  of  the  hot  gas  are  extremely  com- 
plicated, depending  largely  on  the  design  and  arrangement  of 
the  injectors.  For  the  present  analysis,  we  shall  consider 
the  gas  flow  inside  the  combustion  chamber  to  be  an  one- 
dimensional  flow  of  the  hot  gas  only  with  the  unburned  pro- 
pellant elements  suspended  in  and  carried  along  by  the  hot 
gas.  Oscillations  in  the  transversal  plane  normal  to  the 
chamber  axis  is  not  being  considered.  The  hot  gases  all 
over  the  combustion  chamber  are  generated  from  the  combustion 
of  the  same  propellant  under  slightly  different  pressure  and 
thus  are  essentially  at  the  same  stagnation  temperatiire.  The 
combustion  process  in  a liquid  propellant  rocket  motor  does 
not  primarily  increase  the  specific  energy  of  the  gas  system. 
The  heat  released  by  the  combustion  of  a propellant  element 
is  used  to  raise  the  temperature  of  the  combustion  products 
of  this  element  to  the  temperature  of  the  burned  gas  system. 
Therefore  5 from  the  point  of  view  of  the  flow  of  the  burned  • 
gas,  the  combustion  process  is  essentially  a process  of 
generating  or  introducing  new  mass  of  the, burned  gas  into 


the  flow  system  while  the  specific  energy  of  the  gas  system 

remains  substantially  constant.  As  the  combustion  chamber 

is  a duct  of  constant  area,  the  mean  flow  velocity  of  the 

burned  gas  must  increase  whenever  additional  mass  of  burned 

gas  is  introduced  into  the  flow  system.  The  "new”  biJimed  gas 

and  the  "old"  burned  gas  are  assumed  to  mix  intimately  and 

accelerate  together  carrying  with  them  the  suspended  unburned 

propellant  elements.  This  process  of  mixing  and  the  process 

of  accelerating  the  suspended  particles  give  rise  to  an 

•v«n. 

entropy  variation  in  the  gas  flow  field ^if  the  specific  energy 
release  due  to  combustion  is  assumed  constant.  This  entropy 
variation  is  shown  in  reference  5 to  be  a higher  order  small 
quantity  if  the  square  of  the  Mach  ntunbers  of  the  ga:s  flow 
system  is  negligibly  small.  Thus  the  gas  flow  system  can  be 
considered  as  isentropic  to  the  proper  order  of  approximation. 

A concentrated  combustion  front  in  such  a gas  flow  system  does 
not  separate  two  regions  of  gases  having  significantly  different 
thermodynamic  states,  but  is  only  a sharp  discontinuity  of  the 
mean  flow  velocity  of  the  gas.  The  mean  velocity  distribution 
along  the  combustion  chamber  axis  indicates  the  distribution  of 
combustion  while, the  pressure,  the  density  and  the  temperature 
of  the  gas  in  steady  state  operation  are  essentially  uniform 
throughout  the  combustion  chamber  within  the  proper  order  of 
approximation. 

The  shape  of  the  velocity  profile  or  the  distribution  of 
combustdon  in  a liquid  propellant  rocket  motor  varies  con- 
siderably, It  is  often  found  that  most  of  the  combustion,  is 
concentrated  in  a narrow  region.  Therefore,  as  a rough 
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approximation,  we  consider  the  combustion  as  a sharp  dis« 
continuous  front..  In  reference  4,  this  concentrated  com- 
bustion front  is  assumed  to  be  near  the  injector  end  which 
means  the  length  of  the  combustion  chamber  is  much  longer 
than  what  is  necessary*  In  the  present  paper  we  shall  con- 
sider the  concentrated  combustion  front  to  be  located  at 
arbitrary  axial  position  and  the  problem  will  be  formulated 
for  a model  with  two  concentrated  combustion  fronts,  one 
situated  near  the  injector  end,  and  the  other  at  arbitrary 
axial  position.  Thus  we  have  a two-step  velocity  profile 
which  can  be  reduced  to  different  simpler  limiting  cases  of 
special  interest. 

Consider  an  idealized  liquid  propellant  rocket  motor 
whose  injectors  provide  at  a constant  rate  two  uniform  streams 
of  propellant  elements.  The  propellant  elements  in  the  first 
stream  have  a common  small  va,lue  of  total  time  lag  so  that 
these  elements  that  are  injected  into  the  combustion  chamber 
at  the  same  instant  will  burn  simultaneously  at  a place  very 
close  to  the  injector  end.  The  propellant  elements  in  the 
second  stream  are  assumed  to  have  a uniform  total  time  lag  which 
is  much  larger  than  that  of  the  elements  in  the  first  stream 
and  will  burn  at  a distance  '%  from  the  injector  end.  In 
steady  state  operation  the  hot  gas  generated  from  the  first 
concentrated  combustion  front  moves  downstream  with  a velocity 
uf  in  the  region  1 bounded  by  the  two  concentrated  combustion 
fronts.  At  the  second  concentrated  combustion  front  a new 
mass  of  burned  gas  is  introduced  into  the  system  and  thoroughly 
mixed  with  the  burned  gas  from  region  1,  They  move  downstream 
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as  a single  unit  with  velocity  U;j  in  the  region  2 bounded  by 

the  second  concentrated  combustion  front  and  the  prlt  of  the 

combustion  chamber.  Under  the  present  simplified  model  there 

is  no  combustion  taking  place  anyv/here  else  except  at  the  two 

— * — *• 

concentrated  combustion  fronts.  Thus  u,  and  Uj.  are  constants 
in  the  two  regions  indicated  by  the  subscripts  1 and  2 as  shown 
in  Figure  1 (a).  Both  when  U(  — *■  and  when  ^ o 

as  shown  in  Figure  i (b)  and  1 (c)  we  have  the  case  of  a single 
concentrated  combustion  front  at  the  injector  end.  When  u © 
as  shown  in  Figure  1 (d),  we  have  a single  concentrated  com- 
bustion front  at  arbitrary  position  \ . The  dimensionless 
velocities  u,  and  obtained  by  dividing  U|  and  Wi" 
through  the  stagnation  sound  speed  c*  as  a reference  quantity 
are  assumed  to  be  so  small  that  u'  ind  Cf'  , coinciding 
practically  with  the  square  of  the  Mach  numbers  of  the  gas 
flow,  are  negligible  compared  to  unity. 

2.2  The  Time  Lag  and  The  Buming  Rate 

It  is  explained  in  reference  4 that  during  the  total  time 
lag,  the  propellant  elements  undergo  a series  of  complicated 
processes  which  ultimately  lead  to  complete  combustion  of  these 
propellant  elements.  Some  of  these  processes  like  the  atomiza- 
tion of  the  fuel  and  the  oxidizer  and  the  proper  mixing  of  such 
atomized  particles  are  rather  insensitive  to  the  pressure  and, 
the  temperature  oscillations  of  the  burned  gas  in  the  combustion 
chamber.  Many  other  processes  like  the  vaporization  of  the 
propellant  elements  and  the  activation  through  ordinary  heat 
transfer  or  other  means  are  rather  sensitive  to  both  the 
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pressure  and  the  temperature  oscillations  of  the  burned  gas. 

The  total  time  lag  is  therefore  composed  of  a constant  or 
insensitive  part  and  a varying  part  T which  is  sensitive 
to  the  oscillations  of  the  burned  gas  in  the  combustion  chamber. 
Without  any  precise  knowledge  of  these  processes  taking  place 
during  the  period  of  the  time  lag,  we  have  to  assume  some  form 
of  the  dependence  of  T on  the  pressure  and  the  temperature 
oscillations  of  the  burned  gas.  For  the  case  of  small 
oscillations  of  the  burned  gas  about  the  steady  state  con- 
ditions, it  will  be  assumed  that  the  temperature  and  the  pressure 
oscillations  are/ correlated  and  that  the  effect  of  the 
temperature  oscillations  can  be  expressed  in  terms  of  the 
pressure  oscillations.  The  average  rate  of  the  rates  of  the 
different  local  processes  is  assumed  to  be  proportional  to  a 
constant  power  n of  the  local  gas  pressure  acting  on  the 
propellant  element.  Thus,  the  relation  defining  the  pressure 
sensitive  time  lag  X of  an  element  burning  at  the  instant 
t is  given  as  in  reference  4: 
t 

t'  ] dt’  « Con3-^x»nf  C 
t-x 

(2.2.1) 

where  t']  is  the  gas  pressure  acting  on  the 

propellant  element  at  the  position  and  at  the 

instant  t*  . Both  the  constant  value  C of  the  integral  and 
the  constant  pressure  index  n of  interaction  are  characteris- 
tic constants  of  the  propellant  under  the  operating  steady 
state  chamber  pressure.  The  index  n is  assumed  to  be  con- 
stant throughout  the  time  lag  period  and  therefore  represents 
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the  average  extent  of  interaction  between  the  combustion  pro 
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cesses  and  the  pressure  oscillations e The  effect  of  the 
temperature  oscillations  is  also  included  in  the  index  irv  • 

Our  knowledge  of  the  kinetics  of  the  individual  process  is 
not  sufficient  for  a theoretical  prediction  of  such  an  overall 
parameter.  This  constant  Tt  for  given  propellant  can  be  deter- 
mined through  experiments  only. 

Now  we  shall  see  how  this  pressure  sensitive  time  lag 
can  lead  to  varying  burning  rate  in  the  simplified  model  as 
explained  in  the  previous  section.  Consider  for  simplicity 
the  case  of  a single  concentrated  combustion  front.  Let 

be  the  total  amount  of  the  propellant  injected  from 
the  beginning  of  the  operation  up  to  the  instant  t and 
be  the  total  amount  of  the  propellant  burned  up  to  the  same 
instant  t . At  this  instant  t i the  propellant  injected 
into  the  combustion  chamber  during  the  interval  to  t 

has  not  burned.  Therefore 

-mbCt)  = TO  (2.2.2) 

The  instantaneous  burning  rate  is  hence 

dt 

(2.2.3) 

As  we  are  considering  the  intrinsic  type  of  combustion  in- 
stability in  the  idealized  liquid  propellant  as  described  in 
section  2.1,  the  injection  rate  is  a 

constant  equal  to 
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Differentiate  equation  (2»2d)  with  respect  to  t to  obtain 
\ “ ^ and  substitute  \ - ^ into  equation 

(2,2.3) • The  burning  rate  is  thus  given  as 


r ^ ctK  ti 

t-'T  )'] 


(2.2«4) 


Since  the  gas  pressure  varies  with  time  and  position^  the  burn- 
ing rate  is  sometimes  larger  and  sometimes  smaller  than 
which  is  equal  to  the  burning  rate  in  steady  state  operation. 
This  variation  of  burning  rate  produced  by  the  pressure 
oscillations  in  the  combustion  chamber  is  the  self-exciting 
mechanism  of  producing  the  intrinsic  combustion  instability. 

In  the  case  of  low  frequency  oscillation  when  the  length 
of  the  combustion  chamber  is  much  smaller  than  the  wave  length 
of  the  pressure  oscillation,  the  pressure  in  the  combustion 
chamber  is  almost  uniform  at  any  instant.  Thus  j>  [x  W-t) , -t'T] 
is  approximately  equal  to  . In  the  case  of 

high  frequency  oscillations  when  the  length  of  the  combustion 
chamber  is  of  the  same  order  of  magnitude  as  the  wave  length 
or  several  times  larger  than  the  wave  length  of  the  pressure 
oscillations,  it  is  not  obvious  that  js  [v  Ct-t ) , i ” t ] can 
be  replaced  by  XCtJ.t-'tl  as  is  done  in  reference  4« 

However,  it  is  shown  in  reference  5 that  the  spacewise  varia- 
tions of  the  pressure  in  the  combustion  chamber  has  a contri<» 
bution  which  is  a higher  order  small  quantity  in  the  stability 
calculation  compared,  with  the  contribution  of  the  time-wise 
variation  of  the  chamber  pressure  at  a given  location  if  the 
pressure  sensitive  time  lag  T is  much  smaller  than  the 
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pressure  insensitive  time  lag  'T^  . In  the  simplified  model 
which  we  are  considering  now,  T is  much  less  than  Ta  t 
the  spacewise  distribution  of  combustion  is  approximated.  We 
shall  hence  put  1 1 ^ ^ 

and  the  burning  rate  tt)  is  given  as 


(2.2.5) 


where  both  pressures  are  evaluated  at  the  same  location  “X 
where  the  propellant  element  burns.  It  should  be  noticed 
that  when  the  total  time  lag  varies  with  the  gas  pressure, 
the  position  where  the  propellant  elements  bum  also  varies. 
We  shall  neglect  the  oscillation  of  the  concentrated  combus- 
tion front  about  its  steady  state  position  and  evaluate 

|»(t'  and  jaCt-T)  at  the  steady  state  position  of 
the  concentrated  combustion  front. 


! 2.3  Small  Perturbation  Equations  for  the  Gas  System 

In  analyzing  the  low  frequency  oscillations,  the  com- 
bustion chamber  pressure  is  assumed  to  be  uniform  at  any 
instant  but  fluctuates  as  a whole.  This  is  justifiable 
because  the  period  of  the  lovf  frequency  oscillation  is  very 
large  compared  to  the  time  required  for  a pressure  wave  to 

I travel  the  length  of  the .combustion  chamber  so  that  before 

! 

; 4.  the  periodic  oscillation  of  the  gas  has  produced  any  appre- 

ciable changes  in  the  gas  properties,  the  pressure  wave  has 
travelled  many  times  back  and  forth  and  has  made  the  flow 

i 

properties  nearly  -unifoimi.  Thus  the  consideration  , of  mass 
balance  is  sufficient  to  formulate  the  low  frequency 
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oscillation  problem.  In  the  case  of  high  frequency  oscilla- 
tions, the  characteristic  time  for  the  pressure  wav©  to  travel 
the  entire  length  of  the  combustion  chamber  is  comparable  to 
the  period  of  the  oscillations.  Therefore  the  wave  propagation 
phenomena  will  have  to  be  considered  along  with  the  combustion 
process.  The  gas  dynamic  equations  of  continuity,  momentum 
and  energy  will  then  be  used.  From  previous  discussion  in 
section  2,1  the  energy  equation  will  be  replaced  by  the 
equation  of  isentropic  change  of  state  which  is  justified  with- 
in the  proper  order  of  approximation.  Thus  we  have 


I ft*  * 
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(2.3.1) 


where  all  quantities  are  dimensional  and  subscript  \ or 
* 

m means  the  partial  derivative  with  respect  to  the 
corresponding  variable. 

Use  the  following  scheme  to  make  all  these  quantities 
dimensionless, 


%r 


(2.3.2) 


with  superscript  * indicating  that  the  quantity  is  dimen- 
sional and  subscript  ^ indicating  that  the  quantity  is 
evaluated  at  the  stagnation  condition,  is  the  sound 

speed  in  the  stagnant  gas  and  U is  the  length  of  the  com- 
bustion chamber.  The  characteristic  time  is  defined  as 
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the  time  required  for  the  sound  wave  to  travel  the  combustion 
chamber  length  in  a stagnant  gas,  i.e,,  * 

This  characteristic  time  is  one  half  of  the  characteristic 
time  defined  in  reference  4 under  the  approximation  u <<  s « 
Equations  (2,3.1)  when  expressed  in  dimensionless  quantities 
are: 

• " ■!  V«  (2.3.3) 

I (.  f V , , 

These  equations  govern  the  tinsteady  flow  of  the  gas  in  region 
1 and  in  region  2 respectively,  separated  from  each  other  by 
the  velocity  discontinuity  at  the  concentrated  combustion  front. 
For  the  study  of  the  small  oscillations  in  the  gas  system, 
we  shall  consider  the  flow  as  a small  perturbation  over  the 
steady  state  flow.  Thus  define 


where  the  mean  quantities  ^ and  are  practically  unity  in 
dimensionless  form  under  the  approximation  w I , 

The  mean  velocity  U is  constant  either  in  region  1 or  in 
region  2,  Both  u,  and  are  assumed  to  be  small  such  that 
u is  much  less  than  unity.  Introduce  these  perturbations 
into  equation  (2.3o3)  and  linearize  the  equations  with  respect 
to  small  pertui’bations,  we  have 


u, 


o 


h 


z o 


{2.3.4} 
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These  equations  are  recognized  as  equivalent  to  the  simple 
wave  equations  and  admit  solutions  of  the  type  a=  u’(j> 
and  f ' " ?’  ( ^3  with  * t ~ ay  ^ Substituting 

these  into  equation  (2.3»4)»  one  obtains 

/ ( I - u a ) - a.  f =0 

I - ^ 0-  ua.)  = o 

In  order  to  have  non-zero  solutions  of  and  we  must 
have  = o Therefore  only  two  values  of  o- 

are  possible,  i.e.,  » 'Tira  o^s  *”7"^ 

where  -i-  is  the  speed  of  the  dovmstream  moving  wave  and 
is  that  o'"  the  upstream  moving  wave  as  observed  from 
the  combustion  chamber  wall.  Thus  the  general  solution  of 
equations  (2.3.4)  is 

1 u,'  5 -«•  (2.3.5) 


Put  equations  (2.3.5)  into  equations  (2.3.4),  separate  the 
upstream  and  the  dovmstream  moving  waves  and  integrate  with 
the  boundary  condition  that  in  steady  state  both  f’  and  a' 
must  vanish.  We  find  the  relations 


Hence  equations  (2.3.5)  become 


. j’*  u;  - u,'  (t-OsX) 

' u'-  a5(t-ajX3 


Let  UvS  investigate  the  stability  of  periodic  solutions 
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of  exponential  type 

« ft- 1 

1 u s . tt-OsxO  (2.3.7) 

where  and  are  integration  constants  and  ot  ^ >vlu) 
with  X = amplification  coefficient  and  angular  frequency 

of  the  wave.  For  simplicity,  let  us  also ‘write  the  perturba- 
tions as 

u'«  (X)  o|>(at3  , S (V)  (o(t;  o««j  )p  s <|>  (x)  e.x|,  (dt ) 

Then  the  solutions  for  the  functions  and  ? are 

(2.3.S) 

^ tJ»3  « “•  Cj 

' These  solutions  apply  to  regions  1 and  2 respectively.  In 
region  1 the  ratio  of  the  two  integration  constants  c^,  /c,j 
can  be  determined  by  the  boundary  condition  at  the  injector 
end.  In  region  2,  the  corresponding  ratio  can  be  determined 
by  the  boumdary  condition  at  the  combustion  chamber  exit. 

The  two  sets  of  solutions  in  the  two  regions  will  have  to  be 
matched  at  the  second  concentrated  combustion  front  as  re- 
quired by  the  boundary  conditions  at  such  front.  This  match- 
ing of  the  two  sets  of  solutions  defines  completely  the  com- 
plex quantity  , If\is  positive,  the 

disturbance  will  grow  exponentially  v/ith  time  and  therefore 
is  unstable.  If  X is  negative,  the  disturbance  will  die 
out  exponentially  with  time  and  is  stable.  The  boundary  between 
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the  stable  and  the  unstable  regions  as  a relation  between  the 
characteristic  constants  of  the  gas  flow  system,  t ^ 

etc,,  is  obtained  if  we  put  % . The  determination  of 

such  stability  boundary  is  one  of  the  major  objects^  of  the 
present  investigation. 


2.4  The  Boundary  Condition 

The  boundary  condition  at  a concentrated  combustion 
front  will  be  investigated  first.  As  has  been  explained  in 
section  2,1,  the  concentrated  combustion  front  is  not  a 
discontinuity  of  pressure,  density  and  temperature  but  is 
only  a discontinuity  of  velocity  of  the  flow.  Hence  the 
boundary  condition  at  a concentrated  front  consists  of  two 
parts: 


1,  The  steady  state  values  as  well  as  the  small  per- 
turbation values  of  the  gas  pressure  and  the  gas  density  are 
continuous  at  every  instant  across  the  concentrated  combus- 
tion front.  That  is  * f/  . 

These  are  equivalent  to  the  conditions  . 


f ^ - ^1  , * Sy  Gvid 

' S 0* a 3 (d  s X "V 


(2.4.1) 


The  last  equality  in  equations  (2.4.1)  indicates  that  the 
oscillation  frequency  and  the  amplification  rate  are  the  same 
on  both  sides  of  the  concentrated  combustion  front. 

2,  The  fractional  increase  of  the  difference  of  the 
mass  flow  rates  across  the  concentrated  combustion  front  is 
equal  to  the  fractional  increase  of  the  burning  rate  at  the 
concentrated  combustion  front.  The  fractional  increase  of 
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the  burning  rate  )/m,i  can  be  obtained  from  equation 

(2.2.3).  But  in  the  analysis  of  the  simplified  model,  we. shall 
use  equation  (2.2,5).  Thus  neglecting  higher  order  small 
quantities,  we  have: 


- ■wv; 


A'C 
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The  second  boundary  condition  at  the  con  cent rat ed  combustion 
front  separating  region  1 and  region  2 is  obtained  by  equating 

to  the  fractional  increase  of  the  difference 
of  the  mass  flow  rates.  Thus  replacing  the  small  perturbations 
o(«erntdrtl/etB  by  their  periodic  form  and  cancelling  the  common 
factor  ex^{4t)  , we  obtain 


^1  (2.4.2) 

In  equation  (2.4.2)  we  have  replaced  't-  by  x in  the  coefficient 
of  the  small  perturbation  S“,  , neglecting  the  difference 

X ~T!  as  a higher  order  small  quantity. 

If  the  concentrated  combustion  front  is  located  at  the 
injector  end,  the  upstream  side  of  the  combustion  front  has 
no  oscillation.  Hence  the  boundary  condition  at  x=  0 
is  obtained  by  putting  the  disturbance  and  the  mean  velocity 
of  the  upstream  flow  in  equation  (2.4.2)  to  zero. 

-1?,  u,  (»-  JTrt>  S',  ■»-  uj  ifn.  S',  t-dT  )'^  o (2.4,3) 


Now  we  come  to  the  boundary  condition  at  the  combustion 
chamber  exit  where  the  gas  enters  the  converging  section  of 
the  deLaval  nozzle.  The  reflection  of  a one  dimensional 
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pulse  at  the  entrance  of  the  nozzle  is  essentially  three 
dimensional  and  will  lead  to  a complicated  problem  somewhat 
like  the  problem  of  Mach  reflection  from  a wedge.  The  result 
of  such  analysis  even  if  it  could  be  obtained  in  reasonably 
simple  form  is  not  suitable  as  a boundary  condition  for  the 
one  dimensional  flow  in  the  combustion  chamber.  To  be  con- 
sistent the  boundary  condition  must  be  obtained  from  one 
dimensional  consideration,  neglepting  completely  the  two 
dimensional  effects.  The  physical  boundary  condition  at  the 
nozzle  entrance  is  that  there  are  no  discontinuities  of  the 
gas  flow  properties  at  this  vStation.  In  other  words,  the 
iinsteady  one  dimensional  motion  of  the  gas  in  the  combustion 
chamber  and  that  in  the  nozzle  must  be  joined  continuously 
at  this  entrance.  In  reference  6 the  unsteady  flow  in  a 
detaval  nozzle  with  linear  steady  state  velocity  in  the 
subsonic  portion  is  determined  for  arbitrary  frequency  and 
the  boundary  condition  at  the  entrance  to  the  nozzle  is 
presented  in  graphical  form  as  the  ratio  of  the  fractional 
variations  of  velocity  and  density  at  the  entrance.  This 
boundary  condition  when  ^ = \.zo  is  reproduced  in 
figure  Z with  I s "TT^  ® R+<S  plotted  against  the 
reduced  frequency  ^ » cO  • _ u)  '^here  is 

the  length  of  the  subsonic  portion  of  the  nozzle  as  a 
fraction  of  the  combustion  chamber  length. 

If  the  length  of  the  subsonic  portion  is  very  short, 
i.e,,  “S'  O the  values  of  ^ corresponding  to  the 

fundamental  or  the  first  few  modes  of  oscillation  are  very 
close  to  zero.  Hence  for  very  short  nozzle  one  has 
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(2.4.4) 


if  the  one  dimensional  result  holds  good  in  such  limiting  case. 
This  boundary  condition  has  been  shovm  in  references  6 and  7 
to  correspond  to  constant  Mach  number  of  tho  gas  flow  at 
the  entrance  of  the  noaale. 


2,5  Final  Formulation 

Having  established  the  boundary  conditions  for  the 
solutions  in  different  flow  regions , we  can  proceed  to  for- 
mulate the  equation  for  the  determination  of  the  complex 
quantity  o(  > \ -vloU  for  a given  system. 

For  region  l,the  boundary  condition  at  the  injector  end 
as  given  in  equation  (2.4.3)  can  be  used  to  determine  the 
ratio  of  Cyi,  and  C„  , Call  this  ratio  -A  . 


-A.  - -1-- (2,5.1) 


The  solution  in  region  1 can  hence  be  written  as 


. I?,  W s c,,  exy>  A \ - sCCoa.-  a»j  x>3 

I - d CaA~a»)x>] 


In  a similar  manner  the  boundary  condition  at  the  combustion 
exit  Xsb|  can  be  used  to  determine  the  ratio 

Call  this  ratio  ~ B ex']|>  )] 


% 
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When  the  nosale  is  very  short  the  boundary  condition  given  as 
equation  (2,4.4)  is  used,  we  have 


B (o,rx  \ » 


i (BT— < ) /g, 

I - )a/i 


(2.5.4) 


When  the  noazle  is  long  B is  a fiinction  of  ^ and  w.  with 
lCp,M.b  given  in  figure  2, 


0.i=  (2.5^5) 

\ - 1 3 

In  general,  we  can  write  the  solutions  in  region  2 as: 

c-  B exb{(A<«^-as)0~x)i  ] 

• . (2.5.6) 

Q^Cx)  a -Csj,  [ I -V  8 tOA-a*)  ('-“M)  ^ 

Now  the  two  sets  of  solutions  (2.5.2)  and  (2,5.6)  are 
to  be  matched  at  the  second  concentrated  combustion  front 
X * ^ , The  oscillation  of  the  second  concentrated 

combustion  front  about  the  mean  steady  state  position  ^ 
is  neglected.  Equate  S',  (§)  and  as  given  in 

equations  (2,5.2)  and  (2,5.6)  we  have 

—*•  = A <txj»\~oUQA-a>)gl  ^ H- (>a.cA|! 


T 
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where  a,  s ^ *'  been 

substituted. 

By  introducing  equation  (2,5»2),  (2«5«6),  and  (2.5.7) 
into  equation  (2.4.2)  and  dividing  the  resulting  expression 
with  , we  have 

I-  B ^ - A 

BeK^taoid-^)]  A 

' (2.5.S) 

s ( “tj  j](» i-  ^n.  €»j>  ^ *-olf 


This  is  the  final  form  of  the  equation  for  determining  the 
complex  quantity  of  a combustion  system 

with  two  steps  of  concentrated  combustion,  one  at  the  injector 
end  and  the  other  at  ^ . For  the  determination  of  the 

■4 

stability  boundary,  we  put  \«o  ot-  o(  = then 
separate  the  real  and  the  imaginary  parts  in  equation  (2.5.^) 
to  get  two  real  equations  from  which  we  can  eliminate  u)  • 
The  eliminant  is  the  equation  defining  the  stability  boundary. 
We  see  that  in  equation  (2,5,^)  both  A and  B are  in  general 
complex  quantities.  The  separation  of  the  real  and  the 
imaginary  parts  of  equation  (2.5,^)  is  quite  laborious  and 
the  resulting  real  equations  cannot  be  handled  easily. 

For  the  case  of  a single  concentrated  combustion  front 
at  arbitrary  position  ^ (fig.  l»d)  we  put  a,  » «s  in 
equation  (2.5.^).  Thus 


~ Uj,  lit ^ e-K|s(~  ctt  >1 


(2.5,9) 


where  A becomes  unity.  When  the  noazle  is  short, 
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a real  constant  given  by  equation  (2.5*4) • Thus  o(  is  the 
only  complex  qviantity  in  equation  (2.5*9).  The  solution  of 
the  problem  is  therefore  greatly  simplified  algebraically. 
When  the  nozzle  is  long  the  value  of  B is  given  by  equation 
(2.5*5)  where  I is  a complex  function  given  in  figure  2. 

Then  the  determination  of  the  stability  boundary  will  have 
to  be  done  graphically.  If  ^ o , equation  (2.5.9)  be- 
comes 

0- Vn  e'-b  (-at ) = (o  k 

This  case  corresponds  to  a single  concentrated  combustion 
front  at  the  injector  end.  This  equation  (2.5*10)  can  also 
be  obtained  from  equation  (2.5*^)  by  taking  the  limit  of 
either  ^ ® or  » m»  (fig,  lb  and  Ic).  When 

Ktt  a I , equation  (2.5*10)  reduces  to  a form  identical 
with  equation  (12,15)  in  reference  4* 


III.  SOLUTION  WITH  SHORT  NOZZLE 

3*1  Single  Concentrated  Combustion  Front  at  Injector  End 
Revrrite  equation  (2.5*10)  and  drop  the  subscript  2. 


Jfw,  u eTt^(~elT ) ) u 5! 


(3*1*1) 


where 

\ -V  ex’-! ) /a, 

1 •"  (?~i  ^2«5.4) 


■i 


4 


<“  23 


Putting  o(h  yii*>  for  the  neutral  oscillations  and 
separating  the  real  and  the  imaginary  parts  of  equation 
(3,1.1) , we  have 


Cwu)^  = “ t'*'^n)U4  — 
yn.  u ^YTi  tOT  s 


I-  s 


4-  afe  Ce*  aM3 

a<  avO 


(3.1.2) 


I -«•  B*'4  xB  C»s*m3 


With  the  approximation  i*'  < < ) 

- V ~ t<y~< ) a /x I*'  ^ ^ 

V \ -4  [ly-Ow/a]*  “ 

Thus  equations  (3.1.2)  become 

t 

Y»uC«*«JT  r - (l-yn.)  - 

(3.I./4) 

I"  a u c*JT  s 4»w  mJ 


we  can  write 
(3.1.3) 


When  ”(fn.  is  of  the  order  of  unity,  4^nu)  is  of  the 
order  of  XI  ; therefore  we  can  neglect  -4^r.*'  as  com- 

pared to  unity.  Thus  the  first  relation  in  equation  (3.1.4) 
immediately  gives 


Ccs  lOT  - 


't-*- 1 - 


(3.1.5) 


and 


S*n.  s 2:  r »-  ( I-  ^(11  I 
*“  a]fn  ■'  J 


>1 


(3.1.6) 


* - t'-ci-  !C±lf] 


(3.1.7) 


Vx 


From  equation  (3.1*6)  or  (3.1.7)  we  see  that  only  when 


(3.1.^) 


can  we  have  real  solutions  of  a?T  and  u)  for  neutral 
oscillations.  Therefore  the  index  -rv  must  be  bigger  than 
the  minimum  value  combustion  system  is  to 

have  unstable  pressure  oscillations.  As  the  values  of  ^ 
for  most  of  the  combustion  products  of  the  common  rocket 
propellants  are  usually  close  to  1,2  or  1.3 » the  minimum 
value  of  n is  about  0.45  and  is  the  same  for  all  modes  of 
the  high  frequency  oscillations  when  the  constant  Mach  ntamber 
boundary  condition  is  used. 

From  equation  (3.1«7)  we  obtain  the  frequency  of  the 
neutral  oscillations  as 


tO  s k*TT  i t 


(3.1.9) 


where  W • I-  3,  — 


The  fundamental  mode  of  the  high  frequency  oscillation  is  ob- 
tained when  k = 8 ; solution  with  Is  identified  to  • 

be  the  low  frequency  solution  in  reference  4 and  is  therefore 
discarded  in  the  present  investigation.  Equation  (3.1.9)  shows 


that  the  frequencies  of  the  neutral  oscillations  are  close  to 
the  integral  multiples  of  “TT  which  in  dimensional  form  are 
the  natural  organ  pipe  frequencies  of  the  gas  system.  The 
fractional  deviations  of  the  neutral  frequencies  from  the 
corresponding  natural  organ  pipe  frequencies  are  of  the  order 
of  magnitude  of  the  Mach  number  of  the  gas  flow. 

Now  we  shall  demonstrate  on  which  side  of  the  stability 
boundary  is  the  oscillation  unstable.  Taking  B '=  \ and 
differentiating  equation  (3. 1,1)  with  respect  to  -t  we  have 
on  the  stability  boundary, 

(s|o(  . Ju) 

^ — Ail/i  I 't?  — —————  ^ I 

IfftU  Ces*’*o'f-  ^rt,uCos*'»Ox 


ix  . , 

au9  ’ ~ <t  J'nu  - Co* 


(3.1.10) 


When 
sign  of 


, then 


Co#  vox  <0 


, thus  the 


is  the  same  as  the  sign  of  — 4 ir> 


the  sign  of  ~ tO 


When 


is  positive  on  the 


stability  boundary,  the  region  which  is  reached  by  increasing 
<0  from  the  boundary  is  the  \instable  region,  (From  the 
usual  argximent  of  continuity,  this  statement  would  be  expected 


to  hold  good  for  -xv  > 

H 


without  making  an  effort  for 


detailed  proof.)  Therefore,  the  unstable  range  of  the 


frequencies  are 

k-rr  "■  u ^ 1 ^ 

< kTT  ->■ 


(3.1^1) 
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It  is  therefore  concluded  that  all  the  frequencies  of  the 
unstable  oscillations  are  very  close  to  the  natural  organ 
pipe  frequencies. 

The  corresponding  ranges  of  the  values  of  T for  such 
unstable  oscillations  are  given  as 

k-rr-.-  3rnu[<-  (t-  i|l  ^ ^ (3.1.12) 

where  \c=-«,a,  3, indicates  the  successive  modes  of 

the  high  frequency  oscillations  and  o,  \ . x. indicates 

the  successive  higher  ranges  of  the  values  of  t for  un- 
stable oscillations  of  a given  mode.  The  unstable  range  is 
shaded  as  given  in  figure  3a  for  u = o.a<3  . We  see  that 

the  unstable  ranges  of  the  values  of  *t  for  a given  set  of 
values  of  k and  W increases  with  increasing  . 

When  *it/  equals  the  minimum  value  -2^-,  equations 
(3.1.11)  and  (3.1.12) 
gives 

tO  * W-TT 

I u)T  = I j-TT  (3.1.13) 

The  oscillation  is  neutral  and  the  range  of  the  unstable 
values  of  T vanishes  at  this  minimum  value  of  'TU, 

3.2  Single  Concentrated  Combustion  Front  at  Arbitrary  Axial 
Position 

Rearrange  equation  {2^5)9)  and  drop  the  subscript  2, 


•*  -I 


{ 
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I I. 


I j 

i'  ' t 


' i 


T ® [u  en^  (Aip]{; 


(3.2.1) 


Separate  the  real  and  the  imaginary  parts  of  equation  (3.2.1) 
for  the  neutral  oscillation  vdth  XtO 


ytiu  Cl»»*ri'T  = -It-yn^u  ♦ -j- — — 


1 >»ta^;v,uJT. 

Ce*^u)  ( 


(3.2.2) 


Using  equations  (3.1.3)  we  have 


. Vtu*.  Co<  <ait  « - ( '•y*'')  - 

' " ACe,*0-;^uO 


V.  Ynw  vox  a 


■^5«.  tO 


(3.2.3) 


Co*  (\-|^NU> 


By  squaring  and  adding  the  two  equations  in  (3.2.3)  one  ob- 
tains an  equation  from  which  we  can  solve  for  the  frequencies 
of  the  neutral  oscillations.  For  qualitative  discussion  of 
the  results,  we  write  this  equation  as  a quadratic  of 

Cos*(w^jvO  and  solve  for  .iCo**  . Thus 

■4  m «o 


sf  L__.±  . Cv-^ytHy-o  ^ 

^ Coi'^wO  ij-n--!  S’**  J 


(3.2.4) 


♦ ^ u~vw)(».^i)  1^, 

ijn-l  jujn-l  J 

+ Orill! 

ay>x~» 

The  minus  sign  before  the  square  root  is  dropped  because 

is  non-negative.  Introducing  the  following 

two  inequalities; 


t/) 


rt 


C 


f ..  •.  1 


« 
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into  equation  (3.2,4)  we  obtain  the  following  necessary  con- 
dition for  any  possible  real  solutions  of  . 

O-yMcr-o  . o-if  ?'4 

».ra-l  t*“  '■Xyi-I  ■*  s.Yw“l  I > ^ 

(3.2.5) 

which  checks  with  the  minimum  value  of  n as  given  in 
equation  (3.1.S). 

From  equation  (3,2,4)  we  have  another  necessary  con- 
dition for  real  solutions  of  a)  if  only  the  inequality 
X :>>  is  introduced  into  equation  (3.2.4) 


h"”  (3.2.6) 


As  n is  of  the  order  of  unity,  we  must  have 
$;>,  w)  = i ^ • 0(0 

ox  to  » Wtt  Jt  "wi  • o(  •; 


(3.2.7) 


Since  the  unstable  regions  of  for  the  case  of  s o 
is  defined  in  equation  (3.1.11)  we  conclude  that  the  fre- 
quencies of  the  neutral  and  the  small  unstable  oscillations 
for,  all  possible  values  of  ^ are  always  close  to  the 
natural  organ  pipe  frequencies  l<*ir  , 

The  right  hand  side  of  equation  (3.2.6)  vanishes  when 
or  when  C**&tO=:o  . When  either  of  the 

two  conditions  is  satisfied,  must  be  aero.  The 

result  with  V-  is  the  same  as  the  result  obtained 


h 


♦ ^ 
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from  equation  (3.2.5).  The  case  Cos  leads  to  an 

important  restriction  on  the  regions  of  ^ that  admit 

real  solutions  of  u)  , When  both  Sio  and 

are  zero,  Cos  must  vanish,  then  equations  (3*2.3) 

cannot  admit  any  real  solutions  for  , 

From  equation  .(3.2.7)  we  know  tO  is  approximately  w-tt 
therefore  the  zeroes  of  are  approximately 

^ * aU  . ^ $ £ I • These  positions 

correspond  to  the  nodes  of  the  pressure  oscillations  of  the 
k"’*'  mode.  Around  each  of  these  nodes  there  is  a range  of 
positions  ^ of  the  combustion  front  which  is  always  stable. 

The  critical  positions  which  define  such  stable  range  around 
each  node  are  those  where  the  oscillations  are  neutral  with 
» kTT  and  zero  unstable  ranges  of  x . Equation  (3.2.4) 
then  gives  as  . 


I ^ 


' 9 Sc  = 


\ - 


(3.2.^) 


For  the  fundamental  mode,  , there  are  two  values  of 

lying  between  0 and  1 symmetric  with  respect  to  i.  For 
the  mode,  we  have  xk  values  of  between  0 and  1,  de- 

fining W stable  regions  about  the  k nodes  of  the  mode 
of  oscillations.  The  results  are  plotted  in  figure  4 with  a 
v.s.  for  K « 1,  2,  and  3 and  V«  j.A,o 

It  is  observed  that  if  -vt=  the  critical  values 

of  are  0,  ^ ^ etc.,  and  the  entire  region 

of  ^ from  0 to  *1  are  stable  except  at  those  critical  positions 
where  the  oscillations  are  neutral,  and  that  if 

exist.  This  result  agrees  with  equations 


no  real  values  of 


f 
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(3.1.S)  and  (3*2. 5),  and  in  addition  reveals  the  fact  that  if 
the  combustion  is  concentrated  at  the  injector  end  ^ » o 
the  combustion  system  is  most  liable  to  have  unstable  high  fre» 
quency  oscillations.  The  positions  ^ ^ etc.,  are 

approximately  the  anti-nodes  of  the  k"'*'  mode  of  pressure 
oscillations  where  the  amplitude  of  the  pressure  oscillations 
is  the  largest. 

Equation  (3.2.4)  can  be  used  to  solve  the  frequency  of 
the  neutral  oscillation  for  arbitrary  values  of  ^ . But 

with  equation  (3.2.7)  a rapidly  converging  iteration  pro- 
cedure can  be  used  to  solve  the  critical  values  of  'T'  and 
from  equations  (3«2o3)  and  (3.2.4).  By  replacing 

and  with  (.t- ^ ) Ww  and  |,W"ir 

we  can  calculate  -5m  from  equation  (3.2.4)  and  get 

Wtt  dt  J?  u , Using  u)'*  in  the  place  of  , 

fXt.'S  ^ 

we  can  determine  to  « Vctt  ± x ■ and  repeat  the  process 

till  the  necessary  accuracy  is  acquired.  Then  equation  (3.2.3) 
gives  the  values  of  , For  the  cases  calculated,  two 

iterations  are  sufficient.  Calculation  is  carried  out  for 
W = \ . X , and  3 and  U - o . i and  K with 

■yii.  and  U » o. o . The  results  are  plotted 

as'  given  in  figure  5;  These  curves  are  not  symmetric.  The 
shaded  regions  for  a given  mode  of  oscillation  represent  the 
values  of  t and  ^ that  make  this  mode  unstable.  We  see  that 
there  are  only  a few  spots  where  all  the  three  modes  are 
stable. 

From  this  calculation,  it  is  clear  that  the  position  of 
the  concentrated  combustion  front  along  the  combustion  chamber 
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axis  has  considerable  importance.  If  the  combustion  is  mostly 
concentrated  in  a narrow  region  whovse  width  is  only  a small 
fraction  of  the  distance  between  two  adjacent  positions 
of  a given  mode  of  oscillation,  the  stability  behavior  of 
this  mode  can  be  analyzed  with  the  simplified  model  of  con- 
centrated combustion.  If  the  width  of  the  combustion  zone 
is  larger  than  the  distance  between  two  adjacent  of  a 

given  mode,  it  is  hardly  possible  that  the  simplified  model 
could  be  used  satisfactorily.  As  the  distance  between  two 
adjacent  positions  decreases  rather  fast  when  U 
increases,  the  approximate  simplified  model  becomes  less 
satisfactory  for  the  higher  modes  of  oscillation.  Consequently 
if  most  of  the  combustion  in  a liquid  propellant  rocket  motor 
is  concentrated  in  a narrow  region,  for  example  one  tenth  of 
the  length  of  the  combustion  chamber,  the  stability  behavior 
of  the  fxmdamental  and  the  second  mode  of  the  high  frequency 
oscillations  can  be  determined  by  using  the  simplified  model 
of  concentrated  combustion.  But  the  stability  behavior  of 
the  higher  modes  as  analyzed  by  using  the  simplified  model 
should  not  be  considered  too  seriously, 

3.3  Two  Step  Concentrated  Combustion  with  . “TL  = y , 

For  the  study  of  the  nature  of  the  solution  of  the  two 
step  concentrated  combustion  with  one  at  the  injector  end, 
we  shall  consider  the  representative  case  of  "yVi.  » \ for' 
the  purpose  of  simplicity.  Equation  (2.5.^)  becomes 

_ D exj,  ))-  ti- ^ 

Si  ( -- e«p  t-git ) 

(3.3.1) 
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This  equation  can  be  rearranged  in  the  form 

4 c - ^ - B «x|»  < XoO"] 

~ C”ae^^)Cl  + ft€>«pCX<i)](Uf*  u,)^Xlp  ) 

5 ]^  » + B exptxa(l-2^)] 

The  coefficient  of  (-;i*'t ) which  is  of  the  order  of 

(.Ui-wO  u.,  will  be  considered  as  negligibly  small 

compared  with  the  coefficient  of  fyt^C-<^T3  . Thus, 

s, 

CX|9(.alx) 

+ ^ (Ui-Sot  4x|>f  iloco-  f>)l]  D 

(3.3.2) 

Let.  o(s  AuJ  for  the  neutral  oscillations  and  separate  the 
real  and  the  imaginary,  parts  of  equation  (3.3.2)  with  the 
approximation  of  equations  (3.1.3)  , 

C*,  u)t  * ~{|f-.i)a*ti(2x^  ui) 

-Sin  u>t  = 

(3.3.3) 

The  equation  for  determining  the  frequencies  of  neutral 
oscillations  of  the  system  is  obtained  as; 

A"  tO  » A ^ *^*''*’  ^ S 1^" ' ) A M>  Ows® 

"*  if.  i vtj)  Wj  C«^ * A '^ve^ t«j  Ces*  <*33 

Co»  A^?0  •*■  Coj  m)  Co^ (i- 

(3.3.4) 
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It  is  not  easy  to  solve  for  w from  equation  (3*3*4) 
directly-.  But  a method  of  successive  approximation j, 
analogous  to  the  method  used  in  section  3»2  can  be  used. 

Prom  equation  (3»3*4)  we  s®e  that  if  u, - u,  is  small 

compared  to  ui  , we  have  as  a zeroth  approxima- 


.to> 


tion  CO""  a U-rv  Uj 
reference  4.  Introducing 


which  is  the  solution  given  in 


(X> 


.to) 


into  the  right  hand  side 


of  equation  (3.3*4),  we  can  calculate 

M 4-  -0)  ...  I , rr«  — ui  »*' 


u) 


to 


4^^  s'  4 cr*  - ( 


from 


(3.3.5) 


The  solution  of  the  frequency  of  the  neutral  oscillations 
can  thus  be  written  as 


u> 


Wtr  t f.  V4..J, 


(3.3.6) 


where  x is  a function  of  ^ and  ^ . From  equation  (3.3.5) 
it  is  obvious  that  the  coefficient  Z will  be  less  than  unity. 
At  the  nodes  of  the  l<  mode  of  pressure  oscillations,  i.e., 
au  etc.,  equation  (3.3.4)  gives  the  minimum 

value  of  as  provided-  u,  is  not 

too  small. 

Consequently,  the  frequencies  of  the  neutral  oscilla- 

\ 

tions  of  the  system  are  not  significantly  different  from  the 
natural  organ  pipe  frequencies  when  a small  fraction  of  the 
concentrated  combustion  is  shifted  from  the  injector  end  to 
arbitrary  location  along  the  axis. 

Having  determined  u)  , we  can  obtain  tjDx  from  equation 
(3.3.3).  The  result  of  such  calculation  is  given  in  Fig,  6, 
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It  is  noticed  that  when  the  second  concentrated  combustion 
front  is  located  at  the  nodes  of  the  U • mode  of  oscilla- 
tions,  i«e.,  ^ = > C.osui^ 

or  a iPtUfOnrr  ? C"?  - ) 

This  is  identical  with  the  result  for  the  case  of  a single 

concentrated  combustion  located  at  the  injector  end*  However, 

owing  to  the  slight  decrease  of  the  value  of  , the 

ranges  of  the  values  of  t for  producing  unstable  oscilla- 

tions  when  the  second  concentrated  combustion  front  is  located 

at  the  nodes  are  slightly  decreased  as  compared  with  that  of 

a combustion  system  having  a single  concentrated  combustion 

front  at  the  injector  end. 

« 

From  figure  6 we  see  that  the  shift  of  a small  fraction 
of  the  concentrated  combustion  from  the  injector  end  to  any 
arbitrary  locations  along  the  axis  is  slightly  stabilizing. 

The  stabilizing  effect  is  most  significant  when  the  small 
fraction  is  shifted  to  the  node.  Even  in  this  case  this 
shift  does  not  give  rise  to  unduly  large  effect  on  the 
stability  behavior  of  the  system.  It  may  be  inferred  that 
the  distribution  of  a small  fraction  of  combustion  from  a 
given  concentrated  combustion  front  will  not  change  the 
qualitative  picture  of  the  stability  behavior  of  the  system 
as  analyzed  with  the  single  concentrated  combustion  front p 

IV.  SOLUTIONS  WITH  LONG  NOZZLE 

4.1  The  Boundary  Condition  at  Combustion  Chamber  Exit 

The  proper  boundary  condition  at  the  combustion  chamber 
exit  is,  as  mentioned  in  previous  section,  the  continuity  of 
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the  disturbances  in  the  combustion  chamber  and  the  disturbances 
in  the  deLaval  nozzle*  When  the  steady  state  velocity  dis- 
tribution in  the  nozzle  is  linear  from  the  .entrance  to  the 
sonic  throat,  the  ratio  I of  the  fractional  variation  of  the 
velocity  to  the  density  disturbance  at  the  nozzle  entrance  is 
obtained  in  reference  6.  The  real  and  the  imaginary  parts 
of  I when  uxo  are  plotted  against  the  reduced  frequency 
parameter  for  different  valtxes  of  the  reduced  steady  state 
velocity  parameter  » s (fig. 2),  The  reduced 

frequency  jSi  is  defined  as  the  angular  frequency  of  the 
oscillation  divided  by  the  steady  state  velocity  gradient 
in  the  subsonic  part  of  the  nozzle  i.  e. 

Here  is  the  length  of  the  subsonic  part  of  th© 

nozzle  as  a fraction  of  the  combustion  chamber  length. 

For  the  fundamental  mode  of  high  frequency  oscillation, 
we  know  «0  is  about  'Tf  . The  value  of  p corresponding  to 
the  fundamental  mode  is  approximately  « 1 

Since  u is  relatively  small,  the  value  of  ^ for  the  funda- 
mental mode  depends  primarily  on  the  value  of  . In 

practical  eases  is  often  close  to  */3  . If 

of  the  nozzle  is  very  short,  the  value  of 
corresponding  to  fundamental  mode  is  very  small.  Therefore 
if  the  one  dimensional  picture  could  be  extrapolated  to  such 
limiting  case,  the  boundary  condition  of  constant  Mach  number 
would  apply  satisfactorily  to  bhe  case  of  a nozzle  with  a 
very  short  subsonic  part.  On  the  other  hand,  if  is 

reasonably  long,  •'■•he  stability  behavior  of  the  fundamental 
mode  can  be  considerably  different  from  that  determined  bv 
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the  constant  Mach  number  boundary  condition. 

We  shall  consider  a practical  example  vrith  "“1 

such  that  <0  • , Thus  the  fundamental  mode  of 

acoustical  oscillation  corresponds  to  ^ » V and  the 

mode  corresponds  to  ^ - U .It  should  be  noted 
that  for  different  values  of  u , the  condition  oO.s 
requires  that  nozzles  of  different  lengths  are  used  with 
the  combustion  chamber. 


4*2  Solution  for  the  Case  of  A Single  Concentrated  Combustion 
yroni  ' ' ■ ' 

Equation  (2.5.9)  can  be  rewritten  as 

a (l-  yn)-*-  )f»v  exj»  ( -ro/T  ) (4.1.1) 


with  o(»  Au)  for  the  neutral  oscillations  and  B* 

i "•  J lA 

vhere  1=  is  given  graphically  in  figure  2 for  y - \.zo 

as  a function  of  Z ( or  u ^ and  jS  (of  «>)  , The  stability 
boundaiY  is  determined  numerically  for  two  particular  cases, 

Z » 0.0  5 and  0»10  corresponding  to 
u = o.z|3  and  0.301,  For  a given  value  of  ^ and 
a series  of  values  of  the  left  hand  side  of  equation 
(4»1.1)  is  calculated.  Call  this  quantity 


X.  * *T 

J_  f 1 - ^ ^ ^ H 

" wt  I + B exy  I 


(4,1.2) 


Then  the  values  of  n corresponding  to  the  series  of  values 
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of  lO  are  given  by 


-n.  * — +■ 


4 

and  the  corresponding  critical  values  of  T by 


(4.1.3) 


5, 


(4.1.4) 


X, 

gn. 

t-  X 

Ifri*' 

plotted  as  shown  in  figures  7,  d,  9»  10  and  11.  For  the 


The  value  of 
with  C*»  t4T  » (- 


is  taken  in  the  quadrant  consistent 
, The  calculated  results  are 


purpose  of  comparison,  these  curves  are  also  calculated  and 
plotted  when  the  botindary  condition  is  used. 

' It  is  noticed  in  figures  7a  and  7b  that  curves  of 
h vs  with  the  two  different  boundary  conditions  are  of 
similar  shape  and  vary  with  the  same  trend.  The  curve  with 
the  boimdary  condition  I _ I for  long  nozzle 

is  shifted  toward  larger  values  of  n and  smaller  values  of 
p as  compared  to  the  corresponding  curve  with  the  constant 
Mach  nxiraber  boundary  condition  for  very  short  nozzle.  As  a 
result,  there  are  many  apparent  but  important  results  with 
respect  to  the  effect  of  changing  the  ratio  of  the  length 
of  the  subsonic  portion  of  the  nozzle  to  the  combustion 
chamber  length. 

1,  The  minimum  value  of  n compatible  with  any  un- 
stable pressure  oscillations  determined  by  the  long  nozzle 
boundary  condition  increases  rapidly  with  increasing  ^ or  , 
This  is  shown  in  figure  9.  V*hen  the  constant  Mach  number  or 


short  nozzle  boundary  condition  is  used,  the  minimum  value 
of  n for  all  the  different  modes  of  oscillation  is  the 
same  constant  , slightly  less  than  i which  is  the 

minimum  value  of  n for  the  unstable  low  frequency  oscilla- 
tion corresponding  to  W » . 

2,  When  ^ increases,  the  rate  of  increase  of  the 
minimum  value  of  n determined  by  the  long  nozzle  boundary 
condition  is  larger  than  the  corresponding  rate  determined 
by  the  short  nozzle  boundary  condition.  Both  results  are 
plotted  in  figure  4 as  the  dotted  and  the  solid  curves  re- 
spectively. 

3.  For  a given  value  of  n the  unstable  region  of 
X and  ^ as  determined  by  the  long  nozzle  boundarj^  con- 
dition are  smaller  than  the  corresponding  unstable  region 

as  determined  by  the  short  nozzle  boundary  condition.  Owing 
to  the  rapid  increase  of  the  minimum  value  of  n for  in- 
creasing fcO  there  is  an  upper  limit  of  the  frequency  of 
possible  unstable  oscillation  as  determined  from  figure  9 
for  given  value  of  n , This  is  quite  different  from  the 
result  given  in  section  3.2  and  figure  5 where  there  is  no 
such  upper  limit.  Calculation  for  the  case  of  y = e.|33 

and  u,  =o.a.i3  using  the  long  nozzle  boundary  con- 

dition shows  that  the  fundamental  mode  will  be  unstable  when 
the  values  of  •r  and  ^ are.  in  certain  ranges  while  all 
the  other  modes  are  always  stable  (fig, 10).  For  the  case  of 
-yt-I  I both  the  .fundamental  and  the  second  mode  may  be 
unstable  but  not  the  higher  modes  {fig.ll). 

Consequently,  in  every  respect,  the  combustion  system 
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with  long  nozzle  is  more  stable  than  the  one  with  very  short 
nozzle. 

It  should  be  noticed  that  the  rate  of  increase  of  the 
minimum  value  of  n compatible  with  unstable  mode  of 

oscillation  increases  rather  fast  as  Vc  increases.  The 
actual  magnitude  of  n of  the  common  liquid  propellants  is 
not  known.  However,  being  an  exponential  index,  n is  not 
likely  to  be  much  larger  than  unity.  Therefore  for  practical 
cases,  the  higher  modes  of  pressure  oscillations  are  expected 
to  be  stable.  Furthermore,  even  if  these  higher  modes  become 
unstable,  they  are  not  likely  to  build  up  rough  combustions 
owing  to  the  more  effective  viscous  damping  of  the  high  fre- 
quency components.  It  is  the  unstable  fundamental  and  the 
first  few  high  frequency  modes  that  result  in  rough  combustion 
and  hence  are  of  practical  Interest. 

A comparison  of  the  results  in  figures  7a  and  7b  and 
in  figure  9 for  the  two  cases  with  different  u shows  that 
by  increasing  u from  0.213  to  0,301  while  the  nozzle  is  so 
modified  as  to  maintain  the  relation  u)  « , the  minimum 

value  of  n compatible  with  unstable  oscillations  are  slightly 
decreased,  that  is,  the  system  becomes  a little  more  unstable. 
The  change  is,  however,  rather  small  as  produced  by  4055  in- 
crease of  the  flow  Mach  number.  Whether  the  effect  of  in- 
creasing u is  destabilizing  or  stabilizing  is  as  yet  uncer- 
tain because  this  small  difference  in  n can  be  easily  com- 
pensated by  the  texmas  of  the  order  of  u*"  which  have  been 
neglected  in  the  present  analysis.  It  seems,  however,  that 


\- 
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the  stability  behavior  is  rather  insensitive  to  the  change  of 
flow  velocity  in  the  combustion  chamber  when  ~ ] 

is  kept  unchanged. 

As  has  been  mentioned  in  section  3,  if  the  combustion 
is  reasonably  concentrated,  the  simplified  model  of  concentrated 
combustion  front  can  be  used  to  analyze  the  stability  behavior 
of  the  fundamental  and  the  first  few  high  frequency  modes  of 
oscillation,  but  cannot  be  used  to  analyze  the  stability  be- 
havior of  the  higher  modes  with  the  constant  Mach  number  boundary 
condition  at  X*  I • Now  with  the  boundary  condition  for 
long  nozzle  at  x * I all  that  we  are  interested  in  is  the 
stability  behavior  of  the  fundamental  and  the  first  few  high 
frequency  modes  because  all  the  higher  modes  are  expected  to 
be  stable*  Consequently,  the  simplified  model  of  concentrated 
combustion  is  a convenient  model  for  analyzing  the  stability 
of  the  high  frequency  oscillation  if  the  combustion  zone  lies 
in  a region  which  is  small  compared  to  the  combustion  chamber 
length. 

From  the  calculated  results  with  T ~ we  see  that 

if  the  combustion  zone  lies  either  between  » 0 3'='' 

and  ^ or  between  ^ a and 

^ a *.7^  , all  the  high  frequency  modes  are  stable  for 

any  values  of  T when  n is  not  larger  than  1.  The  quantita- 
tive results  as  obtained  in  this  calculation  will  be  somewhat 
different  when  the  steady  state  velocity  profile  in  the  subsonic 
portion  of  the  nozzle  is  not  linear,  but  the  qualitative  results 
as  shown  in  the  present  calculation  is  expected  to  remain  un- 
changed, • 
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an  anti-node  of  all  the  modes  of  oscillations  and  the  con- 

/ 

figuration  with  concentrated  combustion  front  at  the  injector 
end  is  most  liable  to  become  unstable. 

5.  When  the  concentrated  combustion  front  is  at  the 
injector  end,. the  minimum  value  of  n for  the  W""  mode  to 

be  unstable  increases  rapidly  with  increasing  k except  when 
the  nozzle  is  very  short  and  the  constant  Mach  number  boundary 
condition  is  to  be  used  at  the  axial  exit.  The  higher  modes  of 
oscillations  are  relatively  more  stable  than  the  lower  modes 
in  practical  cases. 

6.  Rockets  with  very  short  nozzle  are  more  unstable 
than  similar  rockets  with  long  nozzle.  A rocket  with  short 
nozzle  has  lower  minimum  value  of  n and  larger  unstable 
ranges  of  t . The  qualitative  stability  behavior  of  the  two 
rockets  are,  however,  similar  and  the  results  as  determined  by 
short  nozzle  boundary  condition  is  a helpful  qualitative  guide 
in  the  numerical  determination  of  the  stability  boundairy  with 
long  nozzle  boundary  condition. 

7.  The  stability  behavior  of  a system  with  concentrated 
combustion  is  rather  insensitive  to  the  change  of  the  flow 
Mach  number  of  the  gas  in  the  combustion  chamber  if  the 
velocity  gradient  in  the  subsonic  portion  of  the  nozzle  is 

not  changed; 

,d.  If  the  combustion  is  mostly  concentrated  in  a 
region  whose  width  is  only  a small  fraction  of  the  combustion 
chamber  length,  the  stability  behavior  of  the  fundamental  and 
the  next  few  higher  modes  of  oscillations  can  be  satisfactorily 
analyzed  by  using  the  s5.mplified  model  of  concentrated  ' 
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c imbustion  front.  The  higher  modes  of  oscillations  are  ex- 
pected to  be  stable  when  the  long  nozzle  boundary  condition 
is  used  at  the  combustion  chamber  exit.  If  the  combustion 
is  distributed  so  that  the  combustion  zone  covers  considerable 
portion  of  both  the  stable  region  and  the  unstable  region 
of  ^ of  the  fundamental  mode,  there  is  no  obvious  position 
of  the  concentrated  combustion  front  in  the  simplified  model 
which  can  be  satisfactorily  used  for  the  analysis  of  the 
stability  behavior  of  the  actual  system. 

Since  the  analysis  is  made  on  an  one  dimensional  basis 
in  the  axial  direction,  all  results  apply  to  longitudinal 
oscillations  only. 
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tine  lag,  that  is,  the  proBeure  ineensltire  epaee  lag 


T. 


i 


m 3(1)  • th«  i^crtal  dl«trlbutl^  function  of  p:^pllc2i^t  with 

protfuro  lasoasltl'TO  spupe  log  loa»  th«uo  or  o^biil 
*0 

'4 

Sd*  for  tho  hot  g«LS  oaeiXlfttlon  in  th«  ocadbaatlcm 

ohtoibor 

97a  StthMrlpt  X or  t awoitf  pxrtlxl  dlfforontixtioa  with  respoet  to 

X or  t 


38«  ^ •■  etthstamtiol  dorlmtlvoii.  along  tho  path  of  a propollflait 

' ol«noa^ 

S9»  0 ( ) ■ of  tho  ordor  of  aagoitudo  of  tho  quantity  . in  tho  hraokot 
40a  9i(  ) ■ roal  part  of  tho  quantity  in  tho  hraokot 
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. X ISffiOIKKSnCH  . i 

Cimtast'  .it  i@.  liquM  »s*i!>p9ll«tat  rook*t  me^Qrt  oitm 
' S'm#  oiMuatataxi^aa  m a.  result  of  too  largo  prosvwro 

oMlllatiMS^  i%  tte  ooa^uetioa  ohaaiboro  Thero  aro  two-  digtlaot  ra^gos  of  < 

of  rm^  o«wbugtim^  th@  Iw  frsquonoy  rmg®  @f  loa®  thtm  IW  I 

p»r  MooiSg,  aud  th«  Itigh  froquoaoy  rasgo  of  hmidrodo  or  ! 

thottsaMs  oyolos  p«r  eoamSe  Iho  physieal  roaom  of  tho  solf* 
omittog  ogoill&titma  aad  tho  stability  oritorim  la  euoh.  a oosdiastisii 
ey@t®a  ar®  of  groat  pmotioal  is^|»ort©ao@  m&  haw  baaa  a sabjoot  of  9x.<» 

during  r«iO@at  ^arsa  a^tg  atsthors  haw  asal^od  tMs  ^robloas 

\.  1 

aasaniag  & ocmatakt  tim  dolay  fr^  tho  iastaat  «h«a  a propsllaat 

lajootod  iato  the  oee^uetloa  shassher  to  the  iastaat  vdtiea  this  propsllaat 
elsaeat  burQs«  The  self»«a«itiag  mohaaicH  Is  atsuaMid  to  be  ths  seaaltitlty  j 

■ i 

of  the  fesdiag  system  to  the  >'*ea8ure  oseillatioa  ia  the  eoabustim  ohaaibere  ; 

(leferoaostt  1,2,3 )•  ?he  frequeaoias  of  the  uastable  oseillatloas  of  ths 

ooistbastioa  ohaaibar  ob*tala@d  la  those  refareaoss  aro  la  the  lotr  fre^ueaos^ 

raaga  and  the  poseibllity  of  ha'viag  metable  hi^  frequeaoy  oaoillat^oaa 

has  a©t  beea  iadioat*^* 

!fh@  geaior  aiathor  of  the  pro  seat  paper  has  poiatod  out  the  faot 
that  a pressure  seasitlw  feedlag  systen  is  aot  a aeoessary  Mohaaisra  of 
produeiag  s«lf<»e£Oitlug  osoillatloas  (Sefeireaoe  4)»  fhs  iatemotioa  be»  ; 

tmm  ths  rate  @f  oeabustioa  aad  the  pressure  oeoillatioa  oaii  produee  both 


of  eoi^stim  i^tabllity  has  been  eaXled  intrinslo  instability  sinoe  it 
saaaot  be  eiSjaimtod  by  proper  design  of  th®  feeding  aystem®  1e  Refereaoe  d® 


the  oase  of  Iw  fre^ueaay  ©seillatioaa  la  th®  oojBbustion  systesa  ha®  beeai 
extenslwly  studied  uhll©  the  oas©  of  high  frequeaay  oaeiU&tioas  has  been 
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f »..  < -I 

< .. 


i@4  sodol  witk  & fii&sl<si 

'««Hbttftiion  front  n««ir  th@  iaj@otor  m.A  n ^ • Th«  tr«tt»(Kt  ia 

u 

RtftrfttM  4 hat  bfMn.cxttnM  in.  HMf^rfaoc  5 t®  mor«  gmornl  iM^ts  tsstog 
th«  Mmw  fimplifledl  mod®!  of  oosissatrotoS  ocaibuistiQn  fro&te.  a^ftl 
RoeltlOR.  of  oonMatmtod  'tatestioa  fres&t  hiuB  tooa  fo\m&  to  Iso  iac* 
portont.  ia  dotomlaiae  tho  «t«bilit^  b«h«^or  of  tto  4iff«ir#at  a®d@§  of 
hig^  froquoBoy  otolllotioase  tt  is  thoroforo  euspootoS  that  th« 
bebatior  of  a oyaton  tvith  i«i4«ly  diatributad  oondsustlim  my  b®  ©i^idormbly 
difforoati  tram  that  of  a syataa  with  omamtratod  oeaBbtastiom*  primary 
bb joat  cf  this  paper  ia  to  look  Into  the  offaot  of  dlitrlbKtiag  o«s&*> 
buatloB  along  the  axis  of  the  ocedsastien  ohnBher« 

Xa  the 'present  in's^atlgatlono  the  oeelllatiag  eysten  is  eon» 
aidered  as  a gas  eolona  with  oeatixmous  nass  and  energy  addition  both 
detemiaed  by  the  looal  buraiag  rate*  The  fuel  aad  the  oxldiaer  are  ia^ 
^ted  into  the  ooidsnetion  ohaaiber  ia  the  fom  of  atoodsed  liquid  droplets 
at  a eoBstant  rate  rsgardlsss  of  the  pressure  oaoillatieas  ia  the  ooatbnstioa 
ehsaiber*  The  pattern  of  the  distrihutioa  of  suoh  atonieed  liquid  droplets 
is  aeauMid  to  be  iadepealeat  of  the  pressure  osolllations  so  that  the  fuel 
eleaeat  aad  the  exiditer  eleaeat  that  bt?ra  together  in  steady  stata 
operation  will  still  bum  together  when  the  operation  Is  wopteady*  .Another 
iaportaat  simplifioatioa  ooasists  ia  the  aseuaption  that  eaoh  of  sueh  laixsd 
eleaeats  bums  suddsnly  after  a oertaia  delay  instead  of  gradually*  The 
presence  of  the  gradual  prooesseo  aad  the  effeets  of  the  pressure  sad  the 
temperature'  variations  oa  the  rates  of  the  gradual  proeesses  are  taken 
into  aooount  by  assuming  that  the  delay  ie  affeotod  by  the  variaticms  of 
the  gas  pressure  and  the  gas  temperature  in  the  oonaibustion  ohasaber*  For 
oaoillatlons  of  small  amplitudes » the  temperature  ehange  may  be  related  to 


5:  ajf 


V. 


i. 


th«  preisur*  oh«MS«  of  t;b«  ^s«  The  reXatlsa  beiw&eu  thi@  tim«  d§^y  cmi 
'^ht  pressure  . eheiip  of  the  gas  is  ^ssuKsd  to  he  th@  sem  ae  doeoidheS  in. 
Hef«r«uM  'l«i  M'  a result^  the  total  lag  of  eaeh  eltmmt  defied  m 
the  tine  iattrval  tfm  the  iustaut  of  ihjeetlou  to  the  imtamt  of  ooe^loto 
oisuihustioa  is  iaflu<m9#d  by  the  pressure  ®s@illatious«  The  ehimg®  ®i  tto 
huruiug  rate  resuitiug  frsa  the  ehasige  of  the  total  tiise  lags  of  iifferasit 
eltasmts  uiider  pressure  oioillatioa  is  ths  self-eaceitlmg  msehaaiss  of  the 
unstable  ossi nations* 

In  the  present  inoestigati^  the  problm  is  formulated  M,th 
arbitrary  distribution  of  ooishustion  or  arbitrary  steady  state  wlooity  , 

profile  in  the  oenbustion  ohsnberp  and  au  explioit  approximate  solution 
of  the  final  integro«differential  equation  uith  a retarded  triable  is 
earried  out  for  ^e  ease  of  a linear  steady  state  ^looity  profile  through* 
out  ^e  oonbustion  chsniber  vith  small  final  Maoh  number*  This  linear 
steady  state  selosity  profile  iS  of  partioular  interest  as  it  gives 
approximately  equal  burning  rate  at  all  axial  looations  and  therefore  rep* 
resents  the  extreme  ease  ^i9rm  the  effeot  of  distributing  the  oonbustian 
axially  on  the  stability  behavior  of  the  eoaibustion  system  is  stost  pro* 
nouBsed* 


IX  fORXJUaiQS  OF  les  feobush 

2*1  LLniarised  Small  Perturbatioa  Squatloas 

We  shall  oonsider  the  gas  flow  in  the  oonstsnt  area  omahustien 
ohanhor  to  bo  one  dimensional*  !Rie  injeotors  Ipoated  at  the  end  of  the 
eonibastion  ohamber  provide  at  a ocmstunt  rate  a oontinuous  supply  of  pro« 
pellants  mostly  in  the  form  of  atomised  liquid  droplets  that  bum  at 
different  stations  along  the  oonbustion  ehauiber  axise  It  is  ass^^sasd  that 


>?  *. 


th@  sagliglbl#  rotas  md  &r«  saap»oi#d  im  mi 

sanfta  mlwg  bf  th«  bum«4  f liwiag  lasi4®  tha  aombmtiea 

tha  a;dt«  Fsirt  of  tha  hot  $m  gasaemt^^  irm  etmsibuatim  is 

aotiraiy  to  th«  upsts^ma  ragim  rhara  It  @uppii@@  tha  aotiimtim  m»rg^  t® 

tba  mb:2s%iad  propaOimt  alama.ts«  fha  tamp^ratisra  of  th®  bumaS  gas  throu^- 

I 

out  tha  aefflbuatlosi  ohaabar  is  nasrl^  ooiustmt  &»  all  tha  gas  is  {.smmtad 

t 

froBSi  the  saaa  pfapallmt  uaiiar  essantiany  game  praseui^e  So  S&r  as 
the  £lm  of  i&a  gas  aystiB  is  ooKoemad^  tha  eoHbuiitioa  is  a preoasa  of 
generating  or  ints’oduoing  new  mss  of  tha  bunted  gas  into  the  flow  system 
ettiile  the  apooifio  energy  of  tha  gas  ramalns  nearly  oonstant* 

fhs  thraa  basis  oonsermtion  laws  of  masSe  mmasntm  and  ersrgy 
girs  ths  following  thrss  gas  dyasmio  squatioas  in  dlamsionless  form* 

’ |M^4  yuux  s - «^Ik  f 


(2a.i) 


vt»T9  U|  is  tha  ditaensionlais  relooity  of  the  unbumed  propallat^ 
elementa;  If  the  shift  of  the  dissooiation  reaotione  due  to  the 


pressure  variations  in  the  oombustion  ehsaber  is  neglsated,  the  steady 
state  energy  integral  is  found  to  bt 


(2«leS} 


The  speoifie  entropy  variation  along  a given  burned  gas  elamsnt  Is  found 


^ a u (W -Ujj)  •+  M w <x^] 


(2.l»s) 


r UrUt 


% -t 
^ 1 


$h«  stat#  oratiauiti’  «qui,tio®  »hms  tlmt  it  &£  th»- 

offdi«?  ffif ■ u it  3r«i^®t  t©  k d©«e  a®t 

;9toai|;©  th©  ©r4@r  ©f  of  th®  quantity^  fh@  ©f  tte 

prq||Milltt%ti  is  of  th«  smm  ordtr  ©f  as  th»  n^m  gaa 

u smr  tha  i^ajaetor  an^i.  fhe  ©f  the  sposifi©  ontiropy.  of 

-«  B 

a giwa  pe  alaxmt  at  the  steady  state  is  theseif ore  of  the  order  of  > u 

15  ^ 2i 

©r  ssmilero  Sis&oe  ^ i®  steady  state  ^ tha  speaifie 

entpogy  ^adieat  in  the  ^ - dlreetieas,  * i®  of  ths  order  of 

• thm  steady  state  pressure  aad  density  of  tho  gas  oan  ho 

obtaiaed  ae 

X X 

^ s ) 2 [ ^ ^ ^ f o 1 

?*  * [(-’  ^a’'3^e»^Co(a')] 

(2a*4) 

Thus  it  the  steady  state  gas  YSlooity  u Is  so  saeill  that  Is 

negllglhle  ooapared  to  unity,  the  floor  of  the  gas  in  steady  state  eas.  be 
omsidered  as  isentroplo  and  the  gas  prsssure,  density,  and  temperature  otn 
be  tideen  as  unlfom  throu^out  the  oenbustioa  ohaaber*  Under  the  same  apprexi- 
nation,  the  total  aaouzd:  of  propellants  burned  before  station  x 

par  unit  time  is  approxifflately  ui  « 

Under  the  lialtations  of  small  pertusbations  the  order  of  ma^tuie 
of  these  quantities  mill  not  be  ohanged  by  the  prosenoe  of  the  small 
esoillations#  Thus  the  speoiflo  entropy  variation  in  unsteady  state  opemtlon 
is  negligibly  smXl«  Sinoe  the  speoiflo  entropy  variation  is  not  ooasidered 
as  ^e  basis  meohaniam  of  produolng  self-exoltlng  osolllatlons,  eie  Avail  oon« 
elder  the  one  dimensional  flow  of  the  burned  gas  in  the  rooket  ooabastlon 


f 

fr 


'vditi  omtimaws  mss  additiw  as  ies^itroplo*  7h«  m%r^  «qwti<^ 
in  th@  syetm  gas  dpiastla  aquations  vilX  replaoed  hf  th«  isontropls 
rflatloxi  of  ohan^  of  statse 


t ) 


^>x 

* ~f  (>* 


( 

p f 


Wo  shall  investigato  th@  stability  of  ^11  pos^iodia  dlstufbmsoo 

iatroduood  into  tho  steady  state  gas  flov  gystmc  Gmpl^x  safiablo  is  bss4 
for  semrsAlbaeoii  ¥hus  1st  tho  s»all  portuzbatioas  bo 


jp*  s <iC  t e»|p  (oit.)^ 

f'~  ^ C S (X)  eK|5  ( el-tll 

^ u's  V(x)  exp 


ahoro 


and  V are  ecaploz  funotioos  of  ^ and  ^ s \ 4 lu). 


Yho  dlstttrbsnoos  are  stable,  nsatral  or  unstable  aeoording  as  X ^ ^ 

Thus  the  llnsarissd  smll  perturbation  equatiaas  fc?/  the  eystSK 

of  equations  (2*1«S)  is  eamplsx  Tariables  are 

_ ^ 

s S^-*  s e" 

|lol4axW-*.ju^<4UM*S  s 

Tf  i ^ 
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(2.1.7) 

Where  w(x,  is  th®  perturbations  of  th®  total 

btunaing  rate  before  statioa  • fh®  tem  exp(-ott)  sill  bo 

shorn  later  to  bo  Independont  of  tisss»  By  oliaiinating  and  S in 

oquatiomi  (2*1*7)  and  negleoting  higher  order  smll  quantities  aooordijig  to 
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mm 


iS*'-aM4i7Pi 


•ff 


m ^ a* 
X 

th«t  appreatiatttieiffi  u ■<  < 1 

th»  uidai««a  funotitm  ^ 

V = * |y 

= - If  [( o' » *««  I ^ <■  “ ] 

{2»1.e3) 

In  aqmti^  (3el«3),  Is  still  an  unkn^m 

j^otion  li^ioh  has  to  b«  rolatsS  to  tho  pjrsssurs  vnriatioBi  thirott(|h  tbs 
prcssurs  asitsiti-fity  of  tho  tiao  lag» 

2 .2  fho  gla^  lag  m&  Bgae®  lag 

fli®  prooessos  taking  plaas  during  th®  total  tlat@  lag  Xt  ^^7  ^ 
ditldad  into  two  olas8@s«  the  prouesses  'ehidli  ar®  praotionlly  unaffected  by 

prsssure  and  t«ap®nitur®  -variations  suoh  as  tha  Mining  of  tha  propallantsp 
and  tha  prooassas  ahloh  &ra  Influasioad  by  tliasa  -variations » suoh  aa  haat 
transfer,  -vaporisation,  ohamloal  raaotlonso  As  has  bean  shoan  la  Bafarane®  4, 
aa  shall  approxifflata  -tha  ooMplioatad  sl-tuatlon  by  writing  T 

vd&ar®  represents  a pressure  insansitlve  tins  lag  and  T tha 

prassura  sansltl-va  tisie  lag*  Wa  shall  relate  tha  variation  of  T Only 
to  tha  varlatlona  of  prassura,  aosuolng  that  the  tsnparatura  -varlatltKoo  and 
their  ef fasts  my  ba  oorralatad  with  tha  prassura  variatims*  Ibis  last 
assmption  Is  Justiflabla  espeolally  for  small  osolllationa*  Anothar 
assmption  is  that  tha  rates  of  tha  prsoessas  during  tha  period  of  tha 
prassura  samitl-va  tlm  lag  are  proportional  to  Idia  same  psfwar  n of  tha 
inatantanaous  prassura*  This  exponent  n rapreeents  tha  a-veraga  -value  of  tiie 
different  prooassas  and  Inoludas  the  tmparatura  affeots*  The  relation  wSiloh 
daflnaa  -tho  value  of  X of  an  alomant  that  bums  at  the  instant  t is 


, we  obtsdn  tha  fellwolng  aquation  for 
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Mtai«d  in  i^f«x«no«  4 

rt 

|»  "t*  3 dt  a 

C^aa) 

nh«ff«  th«  pr«®euy®  |j  [ x'ttM  , •<:' ] is  smlKiatsi  al@n|;  th®  gntli  K^Ct*) 

of  the  oleamta  @onsM«r@d«  fh®  omstant  in  oqncitim  (SoSol,)  th® 
pomtat  % nr®  tc  he  th®  same  fer  all  p^psllent 

For  etm-emi^e  -w*  shall  timefom  agnation  C^»Sol)  in  its%a  @f 
the  spatial  mriahle  instead  of  tlme«  this  @f  mriahle  wa 

he  aade  if  tue  1mm  the  aottail  ^leoity  of  the  unhnmed  pr@pel.lant 
during  the  period  • fhe  setlon  of  -^e  eaporlsed  or  the  weil  atondsed 

propellant  eleaents  leasing  the  in^otor  noeele  rolatise  to  the  burned  gae 
near  the  in;)eotor  end  is  mrf  soesi  damped  outc  frm  there  m.  the  motion  of 
the  unbwsned  -prepellant  el@ssnts  &.j  ^ oonsidered  as  approoclnately  the  sene 
M idle  asefage  notion  of  the  jgae  nildioat  being  aff toted  by  idle  hi§^  frequeasy 
eeeillatione*  thus 


u £ (*7 1')  * u 


(2.3o2) 


This  assumption  is  sery  orudei  but  it  has  the  adsaatage  of  allowing  einplsr 
subsequent  developHientsc  TTita  this  assuaptlsn,  we  oan  relate  the  tine  lag 
aad  the  spate  lag  as^lioitly*  C«ssider  a propellant  elesMat  that  ha«  a 
total  tine  lag  ^ *»d  hums  at  the  station  x,  we  hawe 

hy  dofinitlott 


dx 

MU’) 


or>d 


t • 


(2.2.3) 
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sfeere  x*  is  tha  integration  variable  for  x end  | is  the  station  w^ere 


t ^ 

f « r ■■ 


jpvtftttir*  tiw  tog  b«g4zti^s 


Sa  ©s4®r  t©  Imw  di*tslb«ted  ©imbuetiiw,  totol  tla®  lf'4  Xt 
@f  difffswfe  ©l,«Bs®ab®  m«8t  b«  Mffer«ato  MX  i®  by 

wtoMm  (8®^ el)#  th«  mlum  ef  mi  I am  diiiasm^  f@r 


dlff«r®at  gTOg^ltoat  Let  im&M  th«  t®feal  ®mwsrt  ©f  p^©» 

'pfll^’Ss  tisat  haw  b®©@a®  gr®®*!*?®  seasitlw  btfesE’®  th®  statim  | • 

Slaa®  Xi  <si  a gi-mii,  prspeltoffife  is  th®  ®sm®  fw  mr 

stftt®  ©pssmtiea#  | of  th®  will  also  b*  ■a®  ssM'^  3h®»® 

fo7«  5(|S  is  th®  sen®  for  bo^  steady  md  amts&dy  ®tat®  opsimtios#  that 
is  3<^>  * J idiso  J a ^ • H«a9®  th®  fuastlsm 

is  giwn  by 


^C|)=  '&  a in) 


(ZM) 


liMM  X «aA  % aw  r®tot®d  tquatiflos  (2e2el)  sM 


C.otr%'ffv>nf~ 


(S«2»6) 


2e5  Th®  Buraiag  Rat® 

W®  want  to  find  at  Tdtat  station  | in  unet«ady  state  operation 

baglas  th®  pressure  ssnsitiw  tte  tog  of  th®  propellant  elenents  burning  at 

th®  station  x and  at  t^  instant  t«  As  th®  pr®ssar®  sensitl-r®  tin®  tog  T 

is  detomiaod  by  squat  ion  (2«S«l)  with  tto  inte^l  «imluat@d  along  the 

p*ith  s*(t*)  of  th®  propoltont  ®l«Eientg  oonsldored#  th®  pressure  ,t‘| 

is  to  be  ewluatod  at  any  instant  t*  at  th*  station  x*  (t* ) where  pro® 

pellant  elaasRts  were  at  that  Instant  t’#  In  terns  of  the  iptoe  mriabl® 

with  the  instant  t’  is  A£.  earlier  than  t, 

S I , m(c) 
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M >f 
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L vS  ( X ' ) 


cix^  * Cow-S+or»+ 


S (K'J 


5 * (2oSa) 

d| 

Fi«oa  this  •qufttlc90i«  it  oxa  be  ahoim  that  is  aimys  This 

Mans  that  all  slsRents  ahloli  beoona  pressure  sensiti-re  before  the  st&tiem  | 
must  bum  before  the  station  %•  Heno®  if  g(  | ) is  the  mlm  of  the  funsticn 
giwa  b:/  equations  C2*Sa)  and  (2e£«S)  «ith  the  value  of  | defined  by 
equation  (S«S«1),  idie  rate  of  the  total  emount  of  Idie  propellant  eleMsnts 
burned  before  the  station  x at  the  instant  t la  given  by 

W<x,t)  = 9<?)  (2.S.2) 

Proa  equations  (2«2c4t)  and  (2«S*2)  vnt  have 

wfx,ti  ~ wcxis  

ioder  '&#  hypothasis  of  aaall  perturbations  ire  aust  have  If  *■  f l<<  <1 

ISxtts-  the  hi^er  order  Mall  tesMt  in  ^ylor  series  expansion  are  ne 


(2.S.9) 


Pron  equation  (2«8«1)  vith  | W -<  1 ^ found  that 


5-5*1111  C2-3.4) 

^ ^ yn)  I a (XV 

Froa  equations  (2«1«6)  (2«2«i)  (2*2»5)  (2o8*9)  end  (2*5«4)  with  the  appro3d« 

7i  * ^ ^ -1  - Mti  etviS 


— n-l  , , 
■nb  (X  / 


nation 


u <<\ 


, vw  f iad 


-dtt  ^ ^ X 

-e  - TV  u n ^ 


,rw  r #-,1 

^ N’  fcT<5^  1 


^ M I >1 


Xh@  of  th«  total  mto  is  thus  «%plioltl;r  >ss«d  to 

ta»a  of  th«  ps’asstttro  pertus>batl<m  « fbia  i«  th«  quaatltjr  te  1»« 


toio  aqaatlm  (&«1*8)*  It  is  nm  olaar  that  th«  tarn  e>^p(-^t|^X/t  j 
la  todttpmdaat  of  ttoo* 


2*4  yinal  Fogaulatlto 

Bvaluat®  exb(-Jiti  ^ exte(-oa\'5L^ 

' -oX  1 -e  X* 

from  oquation  (2*S«5)  and  anbatitute  thosa  raanlta  into  aquation  (2*7,;c8), 
aa  obtain  tha  following  toto^o-diffarantial  aquation  with  tha  ocnplan 
aaplitttda  of  tha  praaau?®  parturbation  aa  tha  only  dapandant 

aariabla  and  x as  tha  only  indapandant  -tariabla* 


Thia  aquation  ia  too  oonpliaatad  to  ba  aolwad  to  gonasal  tama*  Tha 
quantity  exb  T-  ^ 1 originatad  fron  tha  praaonoe 

! * J M <5  ^ 

of  the  retarded  triable  i - [ , and  aorraapoeds  to  tha 

quantity  €x^  (*  o(T  ) whieh  ia  aseantial  to  dstamtotog  tha  in- 

atability  in  tha  oaaaa  insaatigatad  in  Bafaranoe  4 and  6»  Tha  quantity 

axp  t 1 ii  ^a  only  poeaiblo  eauaa  of  orating 

« 

toatability  of  hi^  fraquanoy  oaoillations  to  tha  oosabustian  aystaso  It 
oan  ba  tbom  that  tha  aolutiona  of  equation  (2«4«l)  with  wmniahtog  ri^t 
hand  aide  hae  no  mistablo  osoillationaa  411  tha  quantities  related  to  oai* 

j 

buetion  are  oentainad  in  tha  right  hand  aide,  and  the  loft  hand  aide  of 
aquation  (2«4«1)  merely  rapraaanta  a one  dtoanalonal  wets  motion  in  a duot 


nmh  m this  of  'ftm  m with  'msx* 

.Ing  erose  mm&<»  fh^@®  ofi^vseti«m  terns  &ni)  m9mmr%lf  of 

or&ox*'  of  hoo&use  the  ox^or  &t  of  a ii  ctiU  Xoft 

?h«  oomhattlon  toms  sr«  quite  eoaplioutoA*  Horn  wo  shall  %ot 
. attoaspt  the  ooluticaji  for  the  ease  with  & g«n®ml  velooi%  dlstrlbutioxi}  hut 
Qoofiao  our  attontlon  to  th«  sinplost  oaso  of  llhoarl^  dlsitfibutod  stswtd^ 
stat*  wwlooity*  This  partieular  ipsleeity  dlstrlhutii^  hriugi  am®  sl»plifl« 
•atioa  la  tb®  otssimstlon  toms  at  the  asaa  tisss  it  stands  as  ea  Saportsat 
sxtswns  oass  just  opposite  to  tho  ono  with  ooaooatsutod  odaibustimi  taoatsd  ia 
Bsfsroasos  4 sad  S«  tho  liataf  walooity  dlstributioa  oaa  ho  writtoa  m 
u (y)  = MX  with  M equal  to  the  Msoh  ausher  of  the  f,as  flow  at  the 

exlt«  y»-l  , uadez*  the  appx'oxiaatioa  0^  <:<.  1 e Thus  the  l@€al 

huraiag  rate  » M is  eonstsat*  ;Hi.ia  speelal  oaiie  thereforf 

eorrespoads  to  the  oae  wllii  ualfenBly  distributed  eosbustlen* 


III  UinFORHLT  SISTRIB13TBD  OO^USTI^ 


5*1  Simplified  Fern  of  the  Peyturbatlsa  Equatl^ 

With  the  steady  state  Tsloolty  prefilo  ixi  » x 
•quatioa  (2*S«S}  beoames 


» jL 

{-  > ^ 
' % 


^ Cx'> 


dx' 


(s.ui) 


The  steady  state  walue  of  the  spaoe  lag  ^ is  fouad  trm  equatisa  (2*E«3) 

08 

I 2 X (i,i®3) 

Siaoe  ths  liaoh  auabsr  of  the  gas  flew  at  th^  oosbustioa  otabef  exit  is  staall|» 


auod  ■X  it.  of  th©  oM©*»  ©f  ^ , it  fts^*  oXoa©  to 

unity*  Ihttt  X~?  * v\]v-  «atp<L~MT^'] 

it  » tnnll  qutatity  o«i|>t}s°©i  to  x«  With  th«  li^mr  'voiooity  dittilbutim 
ti  IX)  * iv)  X «d  th©  ft®9«agrfei0Ei  Ua  »-  ti  , th«  mlu©  of  ^ » T'x 

iMooROf  infinit©*  Stit  it  & omsoqueaoo  of  nogXootins  th©  initinl  diffomio© 
botweon  th©  volooity  of  th©  inj©ot@d  propellant  @lmmtQ  and  a©erags 
©•looity  Of  the  gat*  The  atemptlon  that  x it  nf  th©  order  of  unity 
i©  '^©rafor©  ©quiTaleat  to  assts&pti^  ^t  only  a ©sinll  fraotion  of  th© 
total  tia©  lag-  or  the  epao©  lag  i©  pr©8»ur©  semiitl©©*  With  t of  th© 
.order  of  unity  or  {.v>  1 1 c<  * th©  int©gral  inrolvttd  in  ©quati«m 

(8«1.*1)  ©an  b©  approxinat©ly  •©aluat«d  hy  ©xp«ading  (x'^  in  Taylor 
a©ri©8  about  x* 

If  a©  r©«trl9t  th©  pr©i«nt  ©nalysle  to  t©m©  of  ord©r  IS  and  if 
th©  ord©r  of  n 1©  not  larg©r  than  , all  th©  t©m©  with  HIm  o©©ffi«i«ot 

er  a!^  sea  h#  atgleetfd®  Mfferentiatiag  ©fuatiois  (8«lel)  with  r©>p©et 
t©  X ©•  fiauft 


Tfc<o)s  exp(“<AT)'] /o4 

(8*1*8) 

1 

1 

(8*4*1)  it  reduotd  to  th©  fom 

<=>.  dx  “ O-*- ‘'O  *o 

(8*1*4) 

©1th 

1 M-  ei 

(S«l*5) 

Squation  (S*l*4)  is  virtually  th©  equation  obtained  froa  equation  (2«4*l)  aaad 
(8el*l)  iduKi  ^ (%•)  is  taken  to  b©  <yc*)  , that  1««  th©  ©patial 

variation  of  th©  pr««eur©  p©rt'.isisbatien  hat  negligibl©  oemtribution  in 
©mluating  th©  variation  of  th©  burning  rat©  as  eonpared  vrith  th©  eoatribu» 
tlon  of  th©  tloeuis©  variatiott  of  the  preetur©  at  a given  atatitm* 

Iquatlcn  (8*1«4)  oan  b©  reduo«d  to  ttaMard  fom  of  W©b«r 


'I 


I I ^ 


with  kaossa  ©raanrgeat  ®®rl»8  «ad  &s;?aptotl® 
«oltttio»8'e  H®wwr«  th®  rsit®  of  eeawrgea®©  ®f  if  p©o«* 


ia  the 


©f  iatereot*  ^eref ©re  m transform,  aquatien  (3«1«4:)'  by  tb« 


tfsaiforaation 


($a»6) 


into  a nonaiaesr  Hieoati  Iquation 

V|/  it"  ol  |[  — C%-4-  t ° (?;®1®T) 

This  fonotion  4^  (x)  anst  satisfy  pbysiosi  boimtey  oon&itlm  at 
X a:  o that  the  neon  f lov  Tolooity  u and  the  Aisturbsaoe  'raloeity 
oust  ba  sero  at  any  instant*  frm  equations  (S«l*7)„  so  hat*  ^ 

but  ^(.0  ) , Thus  tho  boundary  onXue  for  ths  funoticn  4 

at  y ~o  ■ is- 


*4/(0  ) . =s  o 


(S.1.8) 


Asymptotic  serlas  solutions  of  aquation  (5*1*7)  in  tarns  of  posars 


of  ^ are  obtalnsd  ast 


3?5^  ' ^ 


(,...  o - 


(S*l«9) 


loithor  of  ths  tno  solutions  satisfy  tho  boundary  oooditim  By 

superposing  the  two  oorrosponding  solutions  of  equation  C3*l«d)« 


the  solution  of  equation  (§el*7)  satisfying  tho  boundary  ooniltlon  (3*le3) 


\i  i 

i ' 

I 

1 t 
] ' 


% i» 

« .; 


til  founiA  b« 





b.  f 4 1^'  ” ^ (sa*io) 

/e  * 


SutesMtntijafi  •quaticms  (3*1*9)  la*®  (S«»1*|,0)  tod  dfbppi^  tb® 


«rd«r  »i»ll,t«ni!S  fto«3®y»iins  as  M*'<< 


■»»  sbiiaia  idi®  soltttiflft  of 


«^bat£«i&  aatisiyiag  th#  bBimdary  mlw»  of  »«•  o foll'osBi 


I •" 


« 

?hl8  soXtttiea  nb®a  ®wluat»d  at  k=1  sa»«t  satisfy  th«  bmmdaffy  9m^ 

dltloc  at  t|u9  eonbustioa  ohaabar  ailt*  Tha  ralatlozi  ©btatead  la  this  ji^aaay 
adll  suffle®  to  datasiaiao  tha  cmplaa:  quaatlty  4,  ^ ' ofa 

givia  syst«Be  For  tha  datarmi nation  ©f  tha  stability  feooiidaifyii  ise  shall 
j^t  V ae  'Shus  Mtpendlag :the  solatl.ac-|3*i;»ll)  sad  lattlag:  x~l 

■m  haao 

^0  ^ * I M [ ^ WtuJ . ^ [ 3 -W*U>3  ^ 

^ ' 


»♦■  -H-  '4t»i7'  uJ  ( < - Co4  >joX  ) • 

xi*i  *A> 

^ ( \ - 4«"’u3  ) »*)X  o-  0»t  USX 

&u3 


I -w«u)+  ^ (3 


" ) -♦■  ■V«k*%*'vO  U>X  ^^*7^ 

3ia)  ^ Xa)  . . W 

•♦’  g~  ( U)  ) Cwi  <<?f  ( Cc»  ti>X  ) I 


(5*1«12) 


8«B  Solatlsaw  CTslag  tog  Fpiqagaay  Boimd 


at  X 


Iho  bouadary  oaodltioa  of  ocastaat  Maoh  saMi#as*  of  «^s  flea 
at  tha  ecaibustlo&  ahsabar  axlt  is  asod  as  m approiitetle«i  l»  RafatHmoa  d* 


4 


f 


4W 


i$-  • 


Jhi#  bpnndiury  o««adltiett  has  been  shew  to  b»  ©orreot  om3.s»‘  f®?  Iw 
freqxmoy  oeelllAtlcme  (Befemio#  6),  end  the  en^osr  ©f  uiAag  tfeii  (Rpgijp®3d*« 
aate  boundary  opudltlon  iia  deteneli^Uag  the  stability  boundary  i$  ©%p#©t©d 
to  ittorease  idtb.  iaoreaslng  frequeaeies  ©f  the  different  Modes  of  ^resshfe 
osoillatloas*  For  ^9  purpose  of  ooeiparisea,  w shall  first  Intestiigite  the 
result  by  using  Idds  lew  frequency  botiadary  oandItloQfi. 

Ihe  eeodltiw  of  eenstent  Maeh  nvoiber  otn  be  writtim  as 

o*'  id  ^ k . ^ -o 

^ f f " 

In  the  ease  of  iseatropio  flour,  we  huwe  ^ S • u x 

By  eUaiiuiting  in  equation  (2*l«f),  ee  get  a'relatisn  SMong  V . ^ 

«t«d  Jhus  the  feemdary  oonditi®&  for  the  funetiiaii  («>. 

at  x»\  is  tbtaiaed  uailer  idie  apprcadmatlat  tA*  <<■  \ as 

foUenS 

X «.  I-  ~ 

Equate  the  two  expretaions  for  vp  O)  froet  equations  (Bnl«lt)  and  (S»B«1) 
and  separate  this  relation  into  the  real  and  the  leaglnary  parts* 


4 _ 

f a ^ r lii  * ^ 


3.  “f  lin  iO  I 1 . 


udiiere 


* IS^ 

_ 'X  Fi  M V , — 

f ~ 


(S.2.B) 


(8.2.S) 


n -I 

-1 


f . 


I 

I 

i 


i 


ir 


fvm  «quAtlons  (3«£»S)  vtt  mn  S»t«min«  th«  arititftl 

'mXatn  ef  T «s<&  th*  required  ®f  'U  o&rx^sgptmding  to 

(frtqufaeloc  &t  aoutml  os«ill9Ltiis3<s  111®  romlt  @f  iis  s^ma  ia 

;ilg^s  S tad  6»  fQT  a dis«u88l@a  of  Idi®  it  if 

ffliffioitntlgr  mooajnit®  to  a«»  th«  simplif i«d  fosansXas 


-i  Wi  U}'t 


W 3 M 


H9  ■•«  obviously  that  for  real  solutioas  of  T and 


(9.2«i) 
tho  miu« 


of  B Bust  bo  biggor  than  sons  lossr  13xlt  abioh  Is  gitoa  approxlBatsly  as 


■ ^ ^ T' V'  M 

(5.8.6) 

For  thoss  values  of  a aot  xuoh  bljg;or  thaa  the  aiaiwm  value  glvaa  ia 
oquatioB  (S.8.6),  as  Bust  have  -V«nM3a-  o(m)  or  <a)=- V<tt  — 

Thus  as  oaa  oonoluAs  that  the  frequeaoies  of  the  oastabls  oseillatloas,  if 


say,  aould  be  elese  to  the  aatural  frequeaoy  V<  'tr 


of  the  gas  systea. 


The  fuadsneatal  uastable  node  has  its  frequenoy  arouad  "W  • Cones* 
queatly,  the  nlainwn  value  of  a aould  bo  of  the  order  of  15  or  80  ahl^ 
sesns  too  large  as  a pressure  esqpaneat.  If  tho  value  of  a of  the  sysiioigi 
should  happen  to  be  this  big«  then  the  uastable  range  of  the  values  of 
sad  't  are  given  approx±eately  by 

tO , c uJ 

OTT-  (*5^  -mM  1 ^ ^ ^ ^ 

“ • - ' ''■■•■  ■ -n  I n II  -T  I --  ^ MMiprI  I l.il  «...  mi 


Ifil^ 


' \ 


iAi«v«  ri  » 0,1,2  • • o aueetfisift  hi^es*  imstabl*  }re»g«a  ot 


tiio  of  -f  » and  W<*  .1,2,3 


• a » 


tho  suaeoiaiina  hi^ar 


»odts  of  th«  oaolllations*  Th«  fundoMUbal  mod*  If  ropraffntad  by  V<  > 1 


fha  boustdaiy  ooaditioQs  at  Urn  aoBdiustioa  obf^ar  mXt  at  all 
froquanoidt  of  oiolllatloas  ar«  tha  ooatiauity  at  aoy  iastaat  of 
T#  loolty,  tho  dasueity,  soKd  tha  praggura  partuzbatioaa  the  ossijnstioa 
ohanber  and  those  la  the  dalaval  nossle*  Zhas  the  hi|^  frequenoy  boaadai^ 


nostle*  In  Hefeaimoe  7,  the  pastioalar  ease  of  unsteady  floir  In  a deLaral 
aossle  adth  linecnp  steady  state  aeloeity  distribution  in  the  subsoaiio  part 
of  the  noiiiile  -is.  snalyiod  ‘srith  th@  bo^dary  os^itisn  of  & nm<= 

•ingular  sonlo  throat*  The  fraetioxtal  density  perturbatioex  for  arbitrary 
frequineiea  is  obtained  by  sol'ting  a hypsrgeei^lsriS'  difforsnti&l 
The  ealoulated  reculte  are  gi^en  in  graphioal  fosa*  The  ratio  of 
is  ealoulated  and  plotted  against  the  redueed  frequenoy 


/9=  = 


eith  Z 


eurvos  as  reproduoed  in  Figures  2 and  3 oaa  be  used  to  detensdne  the 
stability  boundary  of  the  high  frequenoy  esoillatiesu  in  the  ocobustion 
ehaaiber*  Sinoe  the  botmdary  ralue  of  the  oseillation  is  sueh  a otaplioated 
function  of  the  frequenoy  as  gi*ren  by  these  grains,  the  solution  of  the 
oseillation  problem  oan  be  oarried  out  mly  numerioally  for  speoifie  eaamplee 
The  linearised  snail  perturbation  equations  in  te  as  of  the  fraoticnal 
amplitudes  of  the  dietuibanoos  in  the  oopbustion  ohember  are 
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^(1)*  (l-y 
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f SfS'^px 


By  •liolim'liizig 
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Xe 
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Afc  atftticsi  X ■ Ij  ■®e  sh&ll  agts  a®bses*ig'fe  ’F  t©  iadlsals®  bb®  fl«y 
proparbUc  oa  th®  nossl®  aid®  wad  subs@s>ipl;  _ to  ix^©&t®  tii®  flop  ]^ep®rtl®8 
on  th®  ohMib®r  aid®*  thus  th®  bousdazy  aojaAitics?^  at  s * 1 are 


fj 


Qtr»d 


Suhttitatitts  ®%uatiw  (8«1«8)  iat®  ®<ittatlott  (8«8«t)  and  divlllac  thsTMig^ 
, <m®  ebtaiiu 


VV(U  S Y)  “ ) 


(3.8,8) 


idi®r® 

X-iY=  (l-i^)  -U-^)(||-)* 

•nd  F and  a ar®  d®fin®d  by  ®qQatl8e«  (8*2  ^S)  a®  Amotione  of  ^ 71.  wsd 
"C  • is  siTW  by  aquation  (S«I«12)  sad  eaa  be  easily  wittaa  in 

terms  of  &>,  P and  G*  X and  T ar®  funotlea®  ot  P , U ei^  vritH 
eipsa  in  Figaros  2 and  3.  For  gitan  value  of  H the  valu®  of  ^ aorres- 
ponding  to  a dafinit®  angular  fr®qu«noy  u3  inorsasas  pith  dsereasing 
Tolooity  gradiact  , !•©•  with  inoraasing  length  of  th* 

subsoaio  part  of  th®  deJn^l  nossl®  • If  the  subsonio  part  of  th®  nossl®  is 
vary  shorty  th®  ralu®  of  for  the  fuadanontal  mod®  of  ©soillation  is 
vary  soall  and  th®  boundary  oondition  of  ocnstant  Uaoh  nu9d»«r  at  the  ®xlt 


«.  f 


T 


I I 


m 20 


of  ihf  eoatbuKtim  ohcstboi^  imld  b«  & oI»t®  Apps>(!Kisa>ti,@i&«  Wi.m  X6<4>. 

th«  botmda^  ot»idlti«»i  far  aRain&ti<m«  of  sl-^  ftngulw  im^ 
q^ano^  ehiHs@  tomMo  the  bomtdery  oonditicms  far  osaiHatioiaa  of  hi^im 
.redaoed  frequwe^^  ^e  laagth  of  the  subsoala  pert  of  the  deXeml  aoxsl® 
of  roakete  ere  acuelly  l/S  or  l/4  of  the  a6n3>ttetiooi  ohesiber 
fore  la  the  folloelag  oeleuletlea  hs  ehelle  as  a l^ioal.  eaoeaple, 

Uk«  laatuml  freqiueaoieR  oorreapaad  to  lategral.  mlaee 

of  jS  asd  the  fuadeaumtal  orgea  ^|ie  aode  la  gi^rea  by  P » le  ^e  wlaea 
of  X and  Y ae  ealoalated  froa  equatioa  (3*S#4}  are  plotted  ia  Figore  4« 

2 aiad  7 Am  of  the  order  of  uoity  or  aBaller* 

Shut  eepamtiag  the  real  ead  the  imeglsary  giesrts  ia  eq.wti^ 
(9*9«i$}g  ws  Obtain,  tba  folleeiag  two  equati^^  trm  vM.&h  -m  ®5gn  S±^  . 

F fdsd  O.fer  gleea  ealue  of  tha  freqaoaoy  of  the  aeatral  eeeiUatieia* 
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(9*9,8) 


U) 
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^ F U ~ (F  - Cr  >J 


If  a ie  aueh  a xagoitude  that  "‘~Zu  order  of  ualty 

thea  both  F aad  G avet  be  of  the  order  of  uaity^  Eeaoe  ic 

of  the  order  of  M or 

u>  = U-n*  ~ o< M ) 

Shut  the  frequenolei  of  the  aactable  preeeare  oaelXlatloai,  if  aay«  are 
eloae  to  the  aatoral  frequeaoiac  of  the  gas  ayatscso 
The  Taluea  of  a aad  x are  glrea  by 


7nM 

00 


F*'+ 


(9,M) 
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( 


iih*  ef  <jO  € ttsOzes.  ia  th«  quadsmt  oonisi«tfat  witlbt 

CoS  * t-  G-/(Ig£fl) 

Oaiauilfttion  is  «arri«d  oat  f@r  th«  @M«e  with  1 « 0«G@,O«lOi^0eS0 
for  h «*  0 fti«  roKuito  skm  plottM  ta  Flirts  6 @«  For 

th«  ronaltti  ohtf^lMS  by  uiiag  oqmtiem  with  tho  l«sw  fr®%vtto«y  beu^ry 

ooBiiitiw.  at  th©  oMteitiea  ©h«b@r  o^.t  &m  ale©  pl@li®S  im  tteM 


Thosf  oar?o»  mr«  not  intondod  to  b©  ao@amt@  quantitati'«®iy  ospooialiy  for 
iMTgtr  waMos  of  ^ or  1^  itdiioh  ar@  not  n©giigiblj  .Kmll  owparwd  to  un&ty^ 
but  to  show  th«  brand  of  variation  of  diff@ra&t  paramtora  and  'teir  ©rdor 
of  ssaipsitndoo  fheso  oarrsa  ar®  oloooly  amlagons  to  the  oorrespmMag 
ottrves  for  ayatfoa  with  ooaomtratod  e^^stioa.  at  gimi  axial  io«tion« 
(Roformoo  S)*  She  laajor  dlfforonoe  is  that  the  values  of  n for  th®  ®m^ 
of  uniforaly  distrlbutsd  oonbuotlim  is  of  the  order  of  nagnltude  of  iq- 
"^'wfeii.l«  it  i|  of  order  of'  unity  'for  the  oaso  of  osa^tratod  oosdiu@ti€s«* 

Tho  ourrs  of  a wi«  |3  with  tho  hi|^  froquonoy  boundary  oonditloa  art  shlftod 
toward  largor  taluos  of  a and  lower  vaiuss  of  ^ (orbi)  frm  oerros^/ondiag 
ourre  wiidi  tho  low  f roquonoy  boundary  oonditicno  iho  ourv»s  of  • T v«« 
aro  also  shifted  toward  lower  values  of  • 4s  a result^  if  the 
low  froquonoy  boundary  oondltioa  is  used  as  an  apprc^imtioBg,  the  lestabl® 

■4 

rsnge  of  the  values  of  "T  for  a given  value  of  n is  overostisatod  as  shown 
in  Fl^re  7 and  the  wlni,wi  values  of  n ocapatible  with  unstable  oeellla* 
tioos  of  frequonoy  a)  are  uadorestisated  as  shown  in  Figure  9» 

OOHSlUTBlOiB 

The  following  oonolusi<ma  with  respeet  to  tho  stability  b^vior 
of  the  high  frequonoy  osoillatloas  in  a liquid  propellant  roohot  with 
uttifomly  distributed  oombustiou  are  obtained  t 

The  frequenoies  of  the  unstable  pressure  osoillatione.  If 

, ^ t f 

* ♦ 


•xsyg  atr»  elom«  to,  mi  in.  .gonornl  than,  tha^  natura.3.  sfrfquonoloi 

of  tha  t^a®« 

of  a gi-ran  m&&®  in  auoh  a 

ayttom  oan.  beo«^  matahla  i^han  tha  pros  sura  exponont  n ia  blggar  tten 
aartain  Maisaum  mlua  aM  at  tha  sassa  tiM  tha  prassura  aansitiw  ttes  lag 
. -X  of  the  syetm  lias  oas*taisi  Siserot®  ranges  of  mlua@«  fhe  isiinlsim 
f«lu«  of  n sottpatible  with  tha  unatahle  fuadazsaatal  mod®  of  oa^llation  ii? 
of  tha  order  of  15  or  ^0  seems  too  big  t@  b®  % pr@s@ur® 

•Shar®fer®s  it  is  not  -mr^  likely  tSiat  there  eould  b®  tmatabl®-  Mgh  fr®qu«noy 
osoillations  in  sueh  a ooahustion  gyst«ai» 

Sa  fh®  unstabla  ranges  of  the  mluas  of  t of  u glT»n  mod®  of 
oseillfttion  ®r®  fimotlons  ot  n and  th©  m^xitud®  of  aeoh  rang®  inor®aa®a 
a Inorfas® a from  th®  ainljssi  aiilu®  of  n* 

th 

4«  Th®  ninimtoa  tnlu®  of  n of  th®  k mod®  of  osolllatloa  in« 

cr®&s@8  faster  than,  linearly  ■with  iaeroasing  k » Th®  ftmdossKSSs^l  aaj^® 
is  th®  first  to  boooM®  tmatabla  If  T is  in  the  prop«r  range*  Th®  mt®  of 
inofses®  of  this  lainisasa  'mlu®  of  s Is  muoh  largsr  wh®a  the  hi|^  frequanoy 
bou»d«ry  oondltion  is  us®d  than  wh®ii  the  low  fr®qu®noy  boundary  oonditlon 
is  USOdo 

So  The  miaiaum  walue  of  n deoreases  ^rsetioislly  linearly  when 
th®  Toloeity  gradient  M in  steady  state  inoreasesj  that  is,  ^^®n  th®  looal 
burning  rat®  of  the  uniformly  distxdbuted  eraabustioffl  iner®ai$s« 

' d*  The  approstsiation  of  using  the  low  fr«qa«®ey  beun&ary  om® 

dition  in  detemitting  th®  stability  bouadery  results  in  an  ©wsr^estifflatioaii 
of  th®  unstabl®  regions  by  giving  too  lev  minism  valuer  of  n and  too 
largo  unstabl®  ranges  of  the  values  of  T « Th®  error  inoreases  with  the 
higher  modes  of  eseillation*  The  result  of  using  the  low  frequonoy  boundary 


eos^itlea  fer  th®  ftmAanmtftl  »od«  of  osoillfttioa®  l.e,  hoortrwirj,  ® good 
qwilitfttiTe  g«dd®« 

It  sfeoald  t®'  r«mrk®d  fimUy  tliat  th®  f@r@goi^ig  omolusism® 
p»rt®&n  oalj  to  th»  ol®!  modos  ot  ogqillaticns  m tte  smlysla  md® 

on  the  b»*i«  of  o»»  diaonalca®!  flow#  fh®  j»o««lbiIity  of  h®vi.\g  ti^anwml 


nMos  it  not  oocuildorod* 
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